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I.  Research  Summary 

This  final  report  summarizes  research  activities  performed  by  the  Center 
fcr  Accelerator  Science  and  Technology  of  the  University  of  Texas  at  Arlington 
di  ring  the  period  of  July  1,  1987  and  June  30,  1990  under  Grant  No.  AFOSR- 
8',-0280r:>The  main  goals  of  the  successful  research  program  are  to  provide 
theoretical  basis  to  analyze  the  wake-fields  in  the  elliptically  shaped  wake-field 
cavity,  and  to  examine  the  feasibility  of  applying  the  new  idea  of  wake-field 
acceleration  to  beams  of  lasers.  The  theorectical  research  is  very  important 
due  to  the  world-wide  interest  in  finding  new  acceleration  principles  with 
higher  accelerating  gradients  suitable  for  future  particle  accelerators.  In  par¬ 
ticular,  very  exciting  and  new  experimental  result^  from  Argonne  National 
Laboratory  seems  lo*  indicate  that  wake  field  acceleration  is  indeed  a  very  vi¬ 
able  alternative  for  this  endeavor.  \ 

Although  our  own  experimental  work  has  not  been  completed  in  the 
period  due  to  primarily  an  inadequate  experimental  support  for  the  20  MeV 
electron  linac  as  noted  on  the  technical  report  of  July,  1990,  a  great  deal  was 
accomplished  in  the  theoretical  analysis.  We  hereby  report  the  research  ac¬ 
tivities  undertaken  during  this  period.  The  research  performed  under  this  con¬ 
tract' is  divided  into  six  categories .  and  main  achievements  in  each  category  are 
as  follows. 

(1)  Investigation  of  wake-field  acceleration  using  the  elliptical  wake-field 
cavity  _  /  , 

Our  theoretical  works  on  this  subject  were  very  successful  and  major 
progress  has  been  made.  An  analytical  method  was  developed  to  calculate  the 
wake-fields  in  a  meta’lir  elliptical  pillbox  cavity.  By  solving  sets  of  the 
Maxwell’s  equations  in  the  transformed  plane,  we  were  able  to  express  the 
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wake-fields  analytically.  This  analytic  method  simplifies  the  wake-fields  calcula¬ 
tion  significantly,  and  we  can  readily  optimize  the  cavity  design  without  much 
computational  efforts.  Our  analysis  was  shown  to  agree  extremely  well  with 
data  of  existing  3-D  computer  ct  de  as  well  as  with  experimental  data 
presented  recently  by  the  Argonne  ^oup.  This  work  was  well  recognized  and 
published  in  the  refereed  journal  (Particle  Accelerators).  Modal  analysis  has 
been  extended  to  include  the  effects  of  arbitrary  sized  beam  aperture  in  a 
cavity.  Excellent  agreement  between  the  experimental  data  and  our  analysis  was 
demonstrated.  The  second  work  was  also  published  in  refereed  journal  (J.  of 
Applied  Physics,  in  progress). 

-L  (2)  Theoretical  study  of  dielectric  wake-field  cavity 

It  is  found  that  the  dielectric-lined  elliptical  wake-field  cavity  will  yield 
an  acceleration  gradient  of  the  same  order  of  magnitude  that  can  be  obtain 
with  a  metallic  cavity,  and  the  deflection  modes  are  supressed  if  the  energy 
of  the  driving  beam  is  sufficiently  high.  This  new  acceleration  scheme  provides 
a  simple  solution  to  the  transverse  instability  problem  in  the  wake-field  ac¬ 
celeration. 

(3)  Computer  code  development  for  a  wakeTield  cavity  analysis^ 

Based  on  the  theoretical  analyses  developed  here,  very  efficient  and  time 

saving  computer  codes  are  developed. 

(4)  Development  of  wake-field  measurement  system 

A  novel  electrostatic  probe  to  measure  the  wake-fields  directly  was 
developed.  It  was  shown  that  we  can  measure  with  reasonable  accuracy  the 


4 


strength  of  the  wake  potentials  in  a  cavity  without  using  an  expensive  energy 
analyzing  magnet  and  beam  manipulating  equipments. 

(5)  Application  of  Wake  field  cavities:  Study  of  a  soft  x-ray  laser  scheme 

This  work  doals  with  an  application  of  the  elliptical  wake-field  cavity.  We 

have  performed  pioneering  analysis  for  applying  the  wake-  fields  generated  in 
an  elliptical  cavity  as  an  electric  wiggler  in  free  electron  lasers.  Theoretical 
study  of  a  soft  x-ray  EFEL  scheme  using  two-beam  elliptical  pillbox  wake-field 
cavity  has  been  conducted.  It  is  found  that  the  scheme  provides  sufficient  gain 
as  a  coherent  radiation  source  down  to  the  x-ray  regime.  Our  work  was 
presented  in  a  number  of  published  and  refereed  articles. 

(6)  Development  of  a  laser  photocathode 

The  design  and  construction  of  a  laser  photocathode  was  conducted  for 
a  use  as  an  intense,  short  electron  beam  source  for  the  wake-field  acceleration 
study.  Numerical  study  of  optimal  cathode  geometry  was  conducted  and  beam 
currents  were  measured.  It  was  demonstrated  that  this  type  of  photoemitter 
can  be  used  as  an  efficient  short,  intense  electron  beam  source  for  a  compact 
linac  design. 
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II.  Research  Performed  by  Category 

In  the  following,  we  describe  in  detail  the  six  research  tasks  we  ac¬ 
complished.  First  we  list  the  publications  in  which  the  results  were  presented 
in  detail.  A  summany  of  the  results  is  then  given. 

Task  1: 

Investigation  of  wake-field  acceleration  using  the  ellipticalwake-fielcLcavitv- 
Purpose:  to  develop  a  rigorous  and  theoretical  basis  to  calculate  the  wake-fields 
in  an  elliptical  cavity. 

Publications  to-date: 

1.  J.  S.  Yang  and  K.  W.  Chen,  An  analytical  solution  of  wake-fields  in 
an  elliptical  pillbox  cavity,  will  appear  on  Particle  Accelerators,  Vol  23,  1990. 

2.  S.  H.  Kim,  K.  W.  Chen  and  J.  S.  Yang,  Madal  analysis  of  wake-  fields 
and  its  application  to  elliptical  pillbox  cavity  with  finite  aperture,  will  appear 
on  Nov.  1,  1990  issue  of  J.  of  Applied  Physics 

3.  J.  S.  Yang  and  K.  W.  Chen,  Wake  potential  in  a  semi-elliptic  pillbox 
cavity,  on  Advanced  Accelerator  Concepts,  AIP  conference  proceedings  193, 
Lake  Arrowhead,  CA  1989. 

4.  K,  W.  Chen  and  S.  H.  Kim,  Wake-field  acceleration  and  compact  ac¬ 
celerator  considerations,  SPIE  Vol.  875  (SPIE,  Bellingham,  WA  1988),  pp.  223- 
233. 

Results 

Theoretical  study  on  the  wake-field  generation  using  the  elliptical  cavity 
has  been  conducted.  Major  progrers  has  been  made  in  this  study:  development 
of  an  analytical  method,  and  analysis  of  cavity  with  arbitrary  sized  aperture.  It 
was  known  from  theoretical  and  experimental  works  that  the  impedance,  which 
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is  the  Fourier  transformation  of  the  wake-fields,  of  a  closed  cavity  and  that 
of  a  cavity  with  finite  aperture  are  practically  the  same  at  the  low  frequencies. 
Therefore,  the  long  range  wake-fields  of  a  short  bunch  or  the  wake-fields  of 
a  very  long  bunch  are  almost  the  same  for  these  two  different  structures.  Since 
the  wake-field  acceleration  scheme  requires  a  very  short,  intense  driving  bunch 
to  excite  acceleration  fields  of  the  order  of  100  MeV/m  and,  furthermore, 
delay  distance  between  the  driving  and  trailing  beam  is  very  large  compared 
to  the  bunch  length,  we  can  estimate  with  good  accuracy  the  wake  potentials 
of  the  most  practical  cavities  by  neglecting  aperture  effects.  Under  this  as¬ 
sumption  and  using  the  coordinates  transformation,  we  were  able  to  express 
the  wake-fields  and  wake  potentials  of  an  elliptical  cavity  analytically.  With  this 
analytical  method,  we  can  easily  estimate  maximum  energy  gain  of  the  ac¬ 
celerated  particles,  and  simplify  the  optimization  of  the  cavity  design  without 
time  consuming  3-dimensional  computation.  It  is  demonstrated  that  our  analysis 
agrees  extremely  well  with  that  of  existing  computer  code  WELL. 

Although  aperture  effects  are  not  serious  for  the  long  range  wake-  fields 
as  long  as  aperture  is  small  compared  to  the  cavity  cross  section,  we  have  to 
consider  its  effects  to  calculate  the  short  range  wake-fields  or  to  estimate  the 
energy  loss  of  the  driving  beam  and  to  consider  the  beam  dynamics  calculation. 
For  this  purpose,  the  theoretical  analysis  was  expanded  significantly  to  extend 
the  modal  analysis  methods  to  include  the  effects  of  arbitrary  apertures  in  the 
cavity.  A  complete  modal  analysis  was  carried  out  by  using  the  Floquet’s 
theorem  and  an  obvious  requirement  that  the  energy  gain  over  all  acceleration 
cavity  of  repeating  structure  must  be  proportional  to  the  number  of  pillboxes. 
In  addition,  by  assigning  proper  boundary  conditions  on  the  apertures  and  con¬ 
sidering  the  parity  and  continuity  of  the  fields  at  the  aperture,  the  formalism 
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to  determine  the  Fourier  coefficients  for  the  repeating  stricture  consisting  of 
many  pillbox  cavities  with  arbitrary  aperture  was  derived.  The  results  of  this 
analysis  were  compared  with  experiments  done  by  the  Argo  me  Group.  An 
excellent  agreement  between  the  experiment  data  and  o>  r  analysis  was 
demonstrated  as  in  publication  2. 

Task  2 

Study  of  wake-Fields  excited  in  the  dielectric  structure 

Purpose:  to  investigate  the  effects  of  dielectric  material  in  the 

elliptical  cavity  for  a  use  as  a  wake-field  device. 

Publications: 

1.  J.  S.  Yang  and  K.  W.  Chen,  Wake-fields  in  a  dielectric-loaded  elliptic 
waveguide,  submitted  to  Particle  Accelerators  for  publication,  July,  1990. 

Results 

The  wake-field  acceleration  scheme  using  dielectric-lined  elliptical  cavity 
was  conducted.  It  is  well  known  that  the  electron  beam  loses  energy  and  excite 
slow  EM  waves  (Cerenkov  radiation)  behind  when  it  passes  through  a 
dielectric-lined  structures.  The  dielectric  structure  can  be  used  as  a  wake-field 
device  as  well  as  other  applications  such  as  storage  cell,  energy  modulator, 
etc.  It  was  shown  experimentally  by  Argonne  group  that  the  dielectric-lined 
circular  cavity  has  great  advantages  over  other  wake-field  schemes;  acceleration 
gradient  is  of  the  same  order  of  magnitude  that  can  be  achieved  with  metallic 
cavity  and  the  transverse  wake  potential  can  be  made  quite  small  for  the  ultra- 
relativistic  driving  beam.  The  proof-of-principle  of  acceleration  by  the  wake- 
fields  generated  in  the  circular  dielectric  structure  has  been  experimentally 
demonstrated  by  the  Argonne  group.  However,  it  was  not  clear  whether  the 
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transverse  modes  are  also  supressed  in  the  elliptical  dielectric  cavity  as  in  the 
circular  cavity.  Extraction  and  re -acceleration  of  driving  beam  can  be  simplified 
significantly  if  the  driving  beam  and  accelerated  beam  pass  through  different 
paths.  Here,  we  extend  the  idea  of  acceleration  using  dielectric-lined  elliptical 
cavity  to  supress  the  transverse  deflection  modes  and  still  provide  staging 
capability.  The  investigation  of  the  wake-fie’ds  excited  in  the  elliptical  dielectric 
cavity  was  carried  out.  Sets  of  the  Maxwell’s  eqs.  are  solved  analytically  in  the 
transformed  plane.  It  was  shown  that  only  fundamental  modes  are  excited  in 
the  dielectric  structure,  and  all  the  higher  modes  and  transverse  deflecting 
modes  are  inversely  proportional  to  the  square  of  the  driving  beam  energy. 
This  property  of  lack  of  the  transverse  modes  is  a  unique  feature  of  a 
dielectric  structures.  The  model  calculation  showed  that  this  cavity  can  be  used 
as  a  wake-field  acceleration  device  which  has  a  large  acceleration  gradient  of 
the  order  of  100  MeV/m,  no  appreciable  transverse  forces,  and  still  provides 
the  easy  staging  capability.  Excellent  agreement  with  the  Argonne  wake-field 
experiment  was  demonstrated. 

Task  3 

Computer  Codes  Development  for  an  Elliptical  Pillbox  Cavity 
Publications:  J.  S.  Yang,  Results  of  a  PC-based  wake-field  code,  CAST  internal 
report  CAST-89-0115 

Results 

As  a  result  of  the  theoretical  works  developed  here,  a  simple  and  effi¬ 
cient  numerical  code  was  developed  to  calculate  all  the  resonant  modes  and 
loss  parameters  of  an  elliptical  pillbox  cavity.  Excellent  agreements  with  the 
3-dimensional  numerical  code  WELL  as  well  as  with  experimental  data  were 
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Uwmonstrated  in  publication  1. 


Task  4 

Development  of  a  Wake-Field  Measurement  System 

Purpose:  to  develope  a  method  to  measure  the  intensity  of  the  wake-  fie  Js 
directly  in  an  elliptical  cavity. 

Publications: 

1.  S.  H.  Kim  and  K.  W.  Chen,  An  electrostatic  probe  system  for  the 
measurement  of  wake-field  in  a  two-beam  elliptical  pillbox  cavity,  submitted  to 
Review  of  Scientific  Instrumentation 

Results 

A  unique  system  to  measure  the  wake-fields  generated  in  an  elliptical 
cavity  was  developed.  Utilizing  the  facts  that  wake-fields  are  predominated  by 
a  few  lower  order  modes  and  the  wavelength  of  the  fundamental  mode  is  of 
the  order  of  the  cavity  size,  an  electrostatic  probe  with  two  conducting 
spheres  can  measure  the  longitudinal  and  transverse  wake-fields  of  an  elliptical 
cavity.  The  sailent  advantage  of  this  probe  system  is  deduced  from  the  ex¬ 
perimental  result  that  the  voltage  amplitude  at  the  oscilloscope  is  proportional 
to  the  wake-field  amplitude.  The  wake-field  is  about  5  KV  for  3  nC  of  driving 
charge  for  a  cavity  used  in  the  Argonne  group  experiment.  Since  we  can  easily 
measure  a  50  V  pulse,  which  is  the  expected  signal  level  from  the  probe  for 
0.03  nC  driving  charge  of  our  linac,  easily  by  an  oscilloscope,  this  alleviates 
the  need  that  a  powerful  buncher  be  installed  to  a  system  using  a  low-current 
injector  or  the  need  of  a  system  using  a  witness  bunch.  This  measurement 
procedure  will  simplify  greatly  the  measurement  and  test  system  of  the  wake- 
field  experiment. 
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Task  5 

Theoretical  Study  of  a  Soft  X-Ray  Free  Electron  Laser  Scheme 

Purpose:  to  develop  a  soft  x-ray  laser 

Publications: 

1.  S.  H.  Kim  and  K.  W.  Chen,  Soft  x-ray  free  elec  »on  laser  using  a  two 
beam  elliptical  pill-box  wake-field  cavity,  SPIE  Vol.  875  (1988) 

2.  S.  II.  Kim,  Free  electron  lasing  in  a  longitudinal  electric  wave,  Phys. 
Lett.  135 A,  39  (1989) 

3.  S.  H.  Kim,  Stimulated  bremsstrahlung  of  soft  x-ray  in  a  two  beam 
wake-field  cavity,  Proc.  of  International  Conf.  on  Phenomena  in  Ionized  Gases 
(ICPIG)  XIX,  Yugoslavia,  10-14  July  (1989) 

4.  S.  H.  Kim,  Free  electron  lasing  in  transverse  undulating  magnetic  field, 
Phys.  Lett.  135 A,  44  (1989) 

Results 

The  lasing  by  stimulated  bremsstruhlung  of  a  relativistic  dilute  electron 
beam  passing  through  a  spatially  periodic  longitudinal  electrostatic  field  (static 
electric  wiggler)  or  a  traveling  undulating  longitudinal  electric  field  (traveling 
electric  wiggler)  was  investigated.  We  found  that  the  elliptical  pillbox  cavity  is 
suitable  to  generate  the  traveling  wiggler  fields.  If  we  inject  the  lasing  electrons 
in  the  opposite  direction  of  the  driving  beams  which  excite  the  wake-  fields 
along  the  foci  of  the  elliptic  cavity,  the  lasing  electrons  are  forced  to  emit 
the  laser  radiation  by  both  the  longitudinal  and  transverse  wake-fields.  This 
scheme  is  based  on  the  fact  that  we  can  produce  the  longitudinal  and 
transverse  wake-fields  whose  amplitudes  are  more  than  100  MeV/m  in  the  el¬ 
liptical  pillbox  cavity.  Until  the  strength  of  the  laser  wave  arrives  that  of  the 
transverse  wake-field  excited  by  the  driving  beam,  the  transverse  motion  of  the 
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lasing  electrons  is  confined  by  the  transverse  wake-field,  and  both  stimulated 
and  unstimulated  bremsstrahlung  take  place  similar  to  the  radiation  in  the  con¬ 
ventional  magnetic  wiggler.  It  was  shown  that  a  laser  using  a  traveling  electric 
wiggler  provides  sufficient  gain  to  be  used  as  a  coherent  radiation  source  down 
to  the  soft  x-ray  regime.  The  gain  increases  with  the  inverse  of  the  laser 
wavelength,  while  that  of  a  conventional  free  electron  laser  using  a  transverse 
undulating  magnetic  field  operating  in  the  Compton  regime  decreases  with  the 
laser  wavelength.  This  lasing  scheme  is  entirely  new  and  is  quite  important  for 
workers  in  this  field. 

Task  6 

Development  of  a  laser  photocathode  for  the  generation  of  high  current,  short 
electron  bunch 

Purpose:  To  use  in  the  wake-field  acceleration  and  compact  accelerator  design 
study 

Publications: 

1.  K.  W.  Chen,  Y.  C.  Chae  and  J.  Choi,  Development  of  a  laser 
photocathode  for  use  in  wake-field  acceleration  and  compact  accelerator  design 
studies,  SPIE  Conf.  on  Innovative  Science  and  Technology,  Los  Angeles,  CA, 
SPIE  Vol.  875  (SPIE,  Bellingham,  WA  1988) 

Results 

Experimental  study  of  a  back-illuminating  short-pulse  photocathode  for  a 
use  as  a  short,  intense  electron  source  has  been  conducted.  Thin  films  of  CssSb 
was  deposited  on  the  cathode  using  laser  ablation  method.  The  performance 
test  of  the  constructed  cathode  was  made  with  the  Q-  switched  Nd-YAG  laser 
pulse  (532  nm)  directed  from  the  rear  of  the  cathode  surface.  It  was  observed 
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that  a  current  of  2.5  mA/cm2  is  obtained  at  6  KV/cm  potential  in  a  modest 
vacuum  of  10'7  torr.  With  subsequent  improvements  and  refinements  in  design, 
it  was  demonstrated  that  use  of  this  type  can  be  extended  to  other  compact 
linac  designs  for  commertial  applications. 
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Participants  of  the  wake-field  research  program  include:  Dr.  K.  W.  Chen, 
PI,  and  Dr.  S.  H.  Kim,  senior  scientist  and  graduate  students,  Mr.  J.  S.  Yang 
Y.  C.  Chae  and  J.  Choi.  Other  research  students  involved  in  this  program  are 
T.  Pham  and  N.  Nguyen.  Mr.  Yang  completed  his  Ph.D.  degree  in  Mechanical 
Engineering  in  August  1990  with  the  support  from  this  grant.  Mr.  J.  Choi 
obtained  a  M.S.  degree  in  Mechanical  Engineering  in  1988  and  subsequently 
moved  to  the  University  of  Wisconsin  to  pursue  a  Ph.D.  degree.  Mr.  Chae 
has  subsequently  received  a  M.S.  degree  and  has  since  moved  to  the  University 
of  Houston.  He  is  currently  pursuing  his  Ph.D.  degree  at  the  Argonne  National 
Laboratory. 
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V.  Summary  of  Publications 


l.  Journal  Publications 

1.  J.  S.  Yang  and  K.  W.  Chen,  An  analytical  solution  of  wake-fields  in  an 
elliptical  pillbox  cavity,  will  appear  on  Particle  Accelerators,  Vol  23,  1990. 

2.  S.  H.  Kim,  K.  W.  Chen  and  J.  S.  Yang,  Madal  analysis  of  wake-fields  and 
its  application  to  elliptical  pillbox  cavity  with  finite  aperture,  will  appear  on 
Nov.  1,  1990  issue  of  J.  of  Applied  Physics. 

3.  S.  H.  Kim,  Free  electron  lasing  in  a  longitudinal  electric  wave,  Phys.  Lett. 
135 A,  39  (1989) 

4.  S.  H.  Kim,  Free  electron  lasing  in  .transverse  undulating  magnetic  field,  Phys. 
Lett.  135 A,  44  (1989) 

5.  S.  H.  Kim  and  K.  W.  Chen,  An  electrostatic  probe  system  for  the 
measurement  of  wake-field  in  a  two-beam  elliptic  pillbox  cavity,  submitted  to 
Review  of  Scientific  Instrumentation 

2.  Conference  Publications 

1.  S.  H.  Kim,  Stimulated  bremsstrahlung  of  soft  x-ray  in  a  two  beam  wake-field 
cavity,  Proc.  of  International  Conf.  on  Phenomena  in  Ionized  Gases  (ICPIG) 
XIX,  Yugoslavia,  10-14  July  (1989) 

2.  J.  S.  Yang  and  K  W.  Chen,  Wake  potential  in  a  semi-elliptic  pillbox  cavity, 
on  Advanced  Accelerator  Concepts,  AIP  conference  proceedings  193,  Lake 
Arrowhead,  CA  1989. 

3.  K.  W.  Chen  and  S.  H.  Kim,  Wake-field  acceleration  and  compact  accelerator 
considerations,  SPIE  Vol.  875  (SPIE,  Bellingham,  WA  1988),  pp.  223-233 

4.  S.  H.  Kim  and  K  W.  Chen,  Soft  x-ray  free  electron  laser  using  a  two 
beam  elliptical  pill-box  wake-field  cavity,  SPIE  Vol.  875  (1988) 
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5.  K.  W.  Chen,  Y.  C.  Chae  and  J.  Choi,  Development  of  a  laser  photocathode 
for  use  in  wake-field  acceleration  and  compact  accelerator  design  studies,  SPIE 
Conf.  on  Innovative  Science  and  Technology,  Los  Angeles,  CA,  SPIE  Vol.  875 
(SPIE,  Bellingham,  WA  1988) 

3.  Internal  Reports  and  Publications 

1.  CAST-88-0013,  Status  report  on  UTA  Linac  installation  and  testing 

2.  CAST-89-0810,  Wake-fields  in  a  dielectric  loaded  elliptic  waveguide 

3.  CAST-88-0304,  Radiation  shielding  structure  design  and  radiation  safety 
adequacy  considerations 

4.  CAST-89-0115,  Results  of  a  PC-based  wake-field  code 

4.  Thesis  Completed 

1.  J.  Choi,  M.S.  Thesis,  "Design  and  construction  of  a  high  current  electron 
injector  utilizing  laser  back-illuminated  photoemission",  June  1988,  Mechanical 
Engineering  Department. 

2.  J.  S.  Yang,  Ph.D.  Dissertation,  "Wake-field  Acceleration  of  Charged  Particles 
With  an  Elliptical  Cavity",  August  1990,  Mechanical  Engineering  Department. 

V.  Conclusion 

The  research  on  the  wake-field  acceleration  with  an  emphasis  on  the  use 
of  elliptical  cavities  has  been  conducted.  Major  progress  has  been  made  in  the 
analysis  of  the  wake-fields  in  elliptical  structures  and  its  application  for  genera¬ 
tion  of  laser  beams.  It  was  demonstrated  that  the  use  of  elliptically  shaped 
wake-field  cavity  offers  promising  way  for  reaching  acceleration  gradients  of  a 
few  hundreds  MeV/m,  which  is  suitable  for  applications  such  as  compact,  high 
brightness  linacs.  The  elliptic  geometry  allows  us  to  simplify  multi-stage  ac¬ 
celerator  design  and  to  overcome  the  low  transformer  ratio  and  transverse 
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wake-fields  problem.  We  can  achieve  much  higher  acceleration  gradients  by 
selecting  optimal  cavity  geometry  and  beam  distribution.  With  our  analytical 
method,  we  can  simplify  the  optimization  of  the  cavity  design  without  time 
consuming  computational  work.  To  verify  our  theoretical  analyses,  we  compared 
the  results  with  experimental  data  in  circular  wake-field  cavities  obtained  by 
the  Argonne  group,  and  excellent  agreement  are  demonstrated.  Furthermore, 
it  was  shown  that  we  can  obtain  the  wake-  fields  in  a  circular  cavity  exactly 
from  our  analytical  formula  for  the  wake-fields  in  an  elliptical  cavity.  Also,  the 
results  of  our  analytical  method  were  compared  with  existing  numerical  code 
and  very  good  agreement  between  these  two  methods  are  domenstrated. 
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THE  ANALYTICAL  SOLUTION  OF  WAKE-FIELDS  IN  AN 
ELLIPTICAL  PILLBOX  CAVITY 


J.  S.  YANG  and  K.  W.  CHEN 

Center  for  Accelerator  Science  and  Technology ,  The  University  of  Texas  at  Arlington , 

Arlington,  TX  76019 

( Received  June  26.  1990) 

The  wake  potential  of  a  bunch  of  relativistic  charged  particles  traversing  an  elliptical  pillbox  cavity  is 
derived  analytically  in  the  limit  of  vanishing  aperture.  It  is  found  that  the  resonant  modes  of  an  elliptical 
cavity  can  be  expressed  in  terms  of  Mathieu  functions.  Calculation  results  are  presented  and  compared 
with  numerical  ones. 


1.  INTRODUCTION 

The  wake-field  produced  by  a  bunch  of  relativistic  charged  particles  in  an  elliptical 
pillbox  cavity  is  important  not  only  for  future  high  gradient  electron  accelerators1, 
but  also  for  a  use  as  an  electric  wiggler  for  some  proposed  free  electron  laser  schemes2. 
The  suitability  is  based  on  the  estimation  that  the  acceleration  gradient  will  exceed 
100  MeV/m  per  gC  of  driving  bunch  charge,  and  on  wavelengths  of  the  order  of  a 
few  centimeters. 

The  principle  of  acceleration  by  the  wake-field  generated  in  a  metallic  cavity  has 
been  experimentally  verified3.  Also,  other  wake-field  acceleration  schemes  using  a 
plasma  medium4  or  a  dielectric-loaded  cavity5  have  been  experimentally  investigated. 
However,  these  wake-field  schemes  have  not  been  demonstrated  to  have  significantly 
larger  acceleration  gradient  than  that  in  a  metallic  cavity.  Elliptical  cavities  have 
been  investigated  by  several  authors 1,2,6>9.  It  was  shown6,9  that  in  elliptical  cavities 
the  transformer  ratios  are  rather  limited,  and  that  a  strong  transverse  wake-field, 
which  might  lead  to  beam  instabilities,  is  also  excited.  However,  we  can  overcome 
the  low  transformer  ratio  and  the  transverse  deflection  problem  by  using  multi-stage 
schemes.  Although  the  transformer  ratio  would  not  be  as  high,  a  wake-field  acceler¬ 
ator  not  based  on  the  impedance  transformation  principle  could  be  achieved  by  using 
multi-stage  schemes  with  short  stages;  the  driving  beam  is  replaced  with  a  new  one 
or  replenished  in  energy  after  each  stage.  By  rotating  subsequent  groups  of  cavities, 
the  overall  transverse  deflection  of  the  accelerated  beam  can  be  minimized. 

The  wake  potential  in  an  elliptical  cavity  can  be  obtained  either  by  modal  analysis 
or  by  numerically  solving  Maxwell’s  equations  in  the  time  domain.  In  the  previous 
modal  analysis7,8,  the  wake-field  is  expressed  in  a  Fourier  series  based  on  the  vector 
eigenfunctions  of  the  unit  pillbox.  For  infinitely  repeating  structures,  the  problem  for 


the  entire  acceleration  cavity  is  reduced  to  that  in  a  pillbox  cavity  by  using  Floquet’s 
theorem  on  the  periodicity.  The  modal  analysis  can  be  generalized  to  a  cavity  of 
arbitrary  shape  when  we  car  calculate  the  resonant  modes.  For  a  cavity  with  finite 
apertures  we  should  use  numerical  methods  to  find  resonant  modes  accurately. 
Previously  the  wake-fields  n  an  elliptic  pillbox  cavity  with  finite  aperture  were 
calculated  using  the  numerical  code  well9,  which  directly  solves  Maxwell’s  equa¬ 
tions  in  the  time  domain. 

In  this  article,  we  do  not  take  into  account  the  aperture  effect.  The  analytic  solution 
which  will  be  formulated  in  this  article  exists  and  is  of  interest,  even  though  it  is  an 
approximation  for  the  cavity  with  beam  holes,  since  we  can  readily  estimate  the 
maximum  energy  gain  of  the  accelerated  particle.  In  view  of  these  considerations,  we 
try  to  obtain  an  analytic  expression  for  the  resonant  modes  in  an  elliptical  pillbox 
cavity  in  the  limit  of  vanishing  aperture.  Using  the  mode  analysis,  the  longitudinal 
and  transverse  wake  potentials  are  derived  in  terms  of  Mathieu  functions.  It  is  also 
shown  that  we  can  derive  exactly  the  same  expressions  for  the  wake  potentials  in  a 
circular  cavity  as  in  Ref.  7  when  the  ellipse  tends  to  a  circle. 


2.  ELLIPTIC  PILLBOX  CAVITY 

2.1.  Solution  of  Homogeneous  Helmnoltz  Equation  in  an  Elliptic  Cavity 

Consider  an  elliptic  pillbox  cavity  as  shown  in  Fig.  1.  For  a  cavity  of  elliptic  cross 
section,  the  eigenfunctions  can  be  found  in  terms  of  known  functions  by  transforming 
the  Cartesian  coordinates  to  the  confocal  elliptic  coordinates  as  shown  in  Fig.  2.  In 
these  coordinates,  boundary  conditions  on  the  elliptic  cavity  are  readily  satisfied.  The 


FIGURE  1  Elliptic  pillbox  cavity. 
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elliptic  coordinate  variables  (£,  tj,  z)  are  defined  by10 

X  =  h  cosh  £  cos  t], 

Y  —  h  sinh  ^  sin  rj,  (1) 

Z  =  z. 


where  h  is  the  semi-interfocal  distance. 

The  Helmholtz  equation  is  then  transformed  to 


+  +  —  +  «■-»- a 

A2(sinh2  <J  +  sin2  rj)  \  8 <j;2  8 rj2 )  dz2 


(2) 


Following  the  method  of  separation  of  variables,  we  seek  to  find  solutions  of  the  form 

«i,  z)  =  fmnMz).  (3) 

Substituting  Eq.  (3)  into  (2)  and  dividing  by  'f',  Eq.  (2)  is  split  into  three  ordinary 
differential  equations: 

/"(£)  ~{a  —  2q  cosh  2£)/(£)  =  0,  (4) 

g"(rj)  +  (a  -  2q  cos  2fj)gfa)  =  0,  (5) 

w"(z)  +  /?2w(z)  =  0,  (6) 

where  2 q  -  y2h2/2,  y2  =  »c2  —  /?2,  and  a  is  an  arbitrary  separation  constant.  The  above 
Eqs.  (4)  and  (5)  are  called  the  Mathieu  equations.  The  solutions  of  these  equations 
are  the  Mathieu  functions10. 


2.2.  Resonant  Modes  of  a  Cavity 

For  a  closed  elliptic  pillbox  cavity,  we  consider  the  solutions  of 

(Vi  +  y2m,  n)  =  o  (7) 

in  the  region  0  <  q  ^  £0  and  0  <  rj  <2n,  wh  re  1 1/  =  Ez  for  a  TM  wave  and  Hz  for 
a  TE  wave.  Any  combination  of  the  product  of  the  solution  of  Eqs.  (4)  and  (5)  is  also 
the  solution  of  Eq.  (7).  In  addition  to  th-.  boundary  conditions,  the  following 
conditions  must  be  satisfied: 

(i)  continuity  of  \j/  on  the  interfocal  line, 

«A(0,  rj)  =  t/t(0,  —  q),  (8) 

(ii)  continuity  of  gradient  of  i/t  on  the  interfocal  line, 

~  m,  n))(- 0  =  -  4  -n)k- o-  (9) 

oc,  oq 

Among  the  possible  combinations,  the  only  permissible  form  of  the  solution  which 
satisfies  above  two  conditions  is 

•Z'  =  X  CmCeJZ,  q]cem(tj,  q)  +  X  SmSejq,  q)sem(q,  q),  (10) 

m»0  m  = l 

with  the  factor  exp (j(cot  —  /Jr))  being  omitted.  Here,  Cm  and  Sm  are  arbitrary 
constants.  The  functions  cem(q,  q)  and  sejrj,  q)  are  respectively  the  even  and  odd  type 
Mathieu  func  ons  of  the  first  kinds  of  integral  order,  and  Cem(£,  q)  and  Sem(£,  q)  are 
the  modified  Mathieu  functions  of  the  first  kinds  of  integral  order.  These  functions  are 
given  by10 


00 


ce2n{rj,  <?)  =  X  -4</r")  cos  2r7, 

f-0 

(ID 

ce2n  +  l{ri,q)  =  X  cos  (2 r  +  l)q, 

r  *  0 

(12) 

se2a  +  i(n<q)  =  x  B2,0*2'  sin  (2r  +  2 

r.*0 

(13) 

s«2„  +  ifa.  <?)  =  x  «2Vln  sin  (2f  +  l)»f. 

(14) 

r-0 


for  n  =  0, 1, 2, ... ,  and  the  coefficients  A  and  B  are  function  of  q.  Modified  Mathieu 
functions  Cem(£,q)  and  Sem{£,  q)  have  the  same  forms  as  in  Eqs.  <  1 1)— ( 14)  except  q, 
sin,  and  cos  are  replaced  by  <J,  sinh,  and  cosh  respectively. 

We  first  notice  that  symmetry  of  i^(4,  q)  is  determined  by  cem(q,  q)  and  sem(q,  q). 
From  Eqs.  (1 1)— (14)  we  see  that  the  first  term  in  Eq.  (10)  is  an  even  function  and  the 
second  term  is  an  odd  function  with  respect  to  rj.  Therefore,  the  second  term  in  Eq. 
(10)  is  always  zero  on  the  median  plane  (y  =  0  plane)  where  wake  potentials  are  to 
be  evaluated,  and  does  not  contribute  to  the  calculation  of  the  longitudinal  wake 


potential.  For  this  reason,  we  only  consider  the  even-type  modes  of  a  TM  wave, 
ipm  =  Ccm(i ;,  q)cem(t j,  q),  for  the  wake  potential  calculation. 

The  boundary  corditions  are: 


h)  =  0. 

(15) 

There  remains  the  condition  that  no  tangential  component  of  electric  field  exists  at 
the  end-plate  walls  at  z  =  0  and  d,  which  is  satisfied  if  we  choose 

P  =  0,1,2,... 
a 

(16) 

where  d  is  the  gap  distance  of  the  pillbox  cavity  shown  in  Fig.  1.  From  Eq.  (15)  we  have 

Cetfo.  <?)  =  0. 

(17) 

Let  qmn  be  the  nth  root  of  Eq.  (17)  for  mode  m.  Then  we  can  calculate  the  resonant 
frequencies  from  the  root  qmn.  Combining  2 q  =  y2h2/2  and  y2  —  k2  —  /}2,  the  resonant 
frequencies  are  given  by 

v-MiHt)?1 

(18) 

The  wavelength  of  the  dominant  mode  is  then  given  by 

.  7th  nxbee 

^010  f -  f - •> 

v^oi  v  oi 

(19) 

where  ez  and  xb  are,  respectively,  the  eccentricity  and  semi-major  axis  of  the  boundary 
ellipse.  The  ratio  x010/xb  is  plotted  against  ec  in  Fig.  3.  In  this  figure,  we  see  that  as 
ec  -» 0,  i.e.,  as  an  ellipse  tends  to  a  circle,  the  ratio  approaches  2.61,  which  is  the  ratio 
of  wavelength  to  the  radius,  k/rQ  =  27r/y0l,  for  the  TM010  mode  of  a  circular  pillbox 
cavity  of  radius  r0,  where  x01  is  the  first  zero  of  Bessel  function  J0. 

The  field  components  are  given  by 

pn 

35  n)cem(tj,  qmn)  cos  z, 

(20) 

r>  P^  f-  '  i  Z  \  /  \ 

rj  A  Qmn) 

* 

(21) 

p7l  pit 

(22) 

-  n  .  Cem(f,  qmnjcem(rj,  qmn)  cos  z, 
®mnd  d 

(23) 

=  n  J  qm„)cem(rj,  qmn)  cos  z, 

0«.nd  d 

(24) 

=  0, 


(25) 


where 


Dm.  = 


4qnin  (cosh  2<J  -  cos  2>i] ,/: 


(26) 


The  resonant  mode  patterns  for  some  of  the  lower-order  modes  are  shown  in 
Fig.  4. 


3.  WAKE  POTENTIAL  CALCULATION 
3.1.  Longitudinal  Wake  Potential 

Consider  the  test  charge  which  is  traveling  through  the  focus  axis  of  an  elhptic  pillbox 
cavity  (£  =  0,  =  ?r)  and  trailing  the  driving  charge  Q,  which  is  traveling  through 

the  other  focus  axis  {£  =  rj  =  0),  by  a  fixed  distance  s  in  the  r-direction.  The  delta 
function  longitudinal  wake  potential  Wz  is  defined  as  the  energy  gained  by  unit  test 
charge.  Bane  et  al.7  obtained  the  wake  potential  as  an  infinite  sum 

^(s)  =  2  £  fc,  cos 

r  Vt(r')Vx(r)  fa>iS\ 

=  ?  C°S  \~ ) 


for  s  >  0 


(27) 


(C)  (d) 


FIGURE  4  Configuration  of  resonant  modes  in  an  elliptic  pillbox  cavity,  ((a)  m  *  0,  (b)  m  *  1,  (c)  m  =  2, 
(d)  m  «  3). 


where  kx  is  the  loss  parameter  and  Vx  is  the  voltage  induced  by  a  point  charge  Q. 
The  stored  energy  uA  is  given  by 

y  a*  a  =  (28) 

where  aA  is  the  vector  eigenfunction.  The  vectors  r'  and  r  represent,  respectively,  the 
transverse  position  of  the  path  of  the  driving  charge  and  that  of  the  test  charge.  Using 
the  field  components  in  Eqs.  (20H22),  the  voltage  Vx  becomes 

Ki(r)  =  j*  dz  exp  ^"^.(0,  0,  z) 

=  Cem(0,  qmn)cem(0,  qma)  j  dz  exp  (^j  cos  ~  z,  (29) 

K(r)  =  |  dz  exp  n,  z) 

f d  ( j(oz\  pn 

=  Cem(£,  qmn)cem(ii,  qmn) J  dz  exp  cos  —  r, 


(30) 


and  further 


2 Cmn|  - 


vrv,  = 


1  -  (  —  l)p  cos 
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(31) 


where  Cmn  =  Cejq,  qmn)  Cem( 0,  <?mn)  cejq,  qmn)  cej 0,  <?mn). 
The  stored  energy  u1  becomes 
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^mnp  2 


*x*xdV 
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->  m2 

*"  (mn  * 


qmn)cel(>1,  O  cos 
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s«/K 


(32) 


h2 


After  integration  using  dV  =  —  (cosh  2<f  -  cos  2/j)  dc  dq  dz,  it  is  further  reduced  to 


U  = 
mnp  g 


£07lh2d  V  C  ) 


»|2 

•  mn  J 


Ce id  qma)  x  (cosh  2q  -  0  J  dq. 


(33) 


where 


I  f2* 

©n.  =  “  ce2m(q,  qmn)  cos  2 q  dq 
71  Jo 

+  f form =0,2,4,... 


r  *0 


(34) 


l«’)2  +  X  1 4*  3,  for  m  -  1, 3, 5, . . . 


r  *0 


and  A\m)  are  the  coefficients  for  the  series  representation  of  the  Mathieu  function 
cejq,  <?)  g»ven  in  Eqs.  (1  i)  and  tl2). 

Substituting  Eq.  (31)  and  (33)  into  (27),  the  delta  function  longitudinal  wake 
potential  on  the  accelerated  beam  path  (q  =  0,  q  =  rt)  becomes 


J  <30  30  <30 

III 

£0Mm*0»*l  p  ■  0 


SpC^l  -(-l)pcos?^) 
- L  cos 


qmn  N  ff 


where  e.  =  \  for  p  =  0  and  1  for  p  #  0.  Nmn  and  Cm„  in  Eq.  (35)  are  given  by 


Nmn  = 


Ce2JZ,  qmn)  *  (cosh  2<J  -  0  J  dq. 


(35) 


(36) 


Cm»  =  (CeJ  0,  ^mn  ))2«m(0,  q 

[(f^r)  form  =  0,:, 4,... 

form  =  1,3,5 — 

For  a  circular  pillbox  cavity,  Wz  is  expressed  analytically7'1*  in  the  form  of 


A  an  --u 

K(s)  =  — ;  X  X 

n  »  !  p  =  ( 


—  ( —  l)p  cos  — j 


Zon*f  l(Xon) 


for  m  =  0 


where  yon  is  the  nth  zero  of  the  Bessel  function  JQ. 

When  the  boundary  ellipse  tends  to  a  circle  of  radius  r0,  the  confocal  hyperbolae 
in  Fig.  2  become  radii  of  the  circle  r,  and  the  confocal  ellipses  become  concentric 
circles  of  that  radius.  In  this  case,  Eq.  (35)  is  reduced  exactly  to  Eq.  (38).  We  can 
easily  show  this  by  using  the  limiting  properties  of  Mathieu  functions.  The  Mathieu 
functions  cem(q,q)  and  Cem(£,  q)  degenerate  into  the  following  forms'  as  h-»0 
and  c  -*  30,  while  keeping  the  product  h  cosh  -*•  r: 

q  -*  0  as  h  -*  0,  (39) 

(Ji,  form  =  0 

cem(n,  9)  “*  |CQS m(p  form/0 

Cem&  q)  —  pmJm(yr),  (41) 

A[m)  -♦  0  (except  /4(™'  -+  1  for  m  #  0,  and  A{°'  -*  %/!),  (42) 

where  pm  is  a  constant  multiplier,  and  Jm(x)  is  the  Bessel  function  of  the  first  kind. 
Then,  Eq.  (37)  becomes 

Cmn  =  (Cem(0,  <?m„))2cem(0,  q„jcejn,  qmn) 

-*•  (pm3m(0)  cos  mtp)2  =0,  for  m  #  0 
-  (PoJom2/2  -  PS/2,  for  m  —  0.  (43) 

It  is  apparent  from  Eq.  (43)  that  contributions  from  m  ^  0  modes  become  zero 
as  expected  as  an  ellipse  tends  to  a  circle.  The  denominator  in  Eq.  (35)  for  m  *  0 
becomes 


0 

,  CeKt 

Jo 


q„n)  cosh  2S,  dt, 


£j> 


Jl( y&  dr 


■  jXoVfrU. 


where  y  =  yoJr0.  Substituting  Eqs.  (43)  and  (44)  into  (35),  we  can  get  exactly  the 
same  expression  for  the  wake  potential  in  the  circular  pillbox  cavity  as  in  Eq.  (38). 


No  closed  expression  is  known  for  the  infinite  sum  in  Eqs.  (35)  and  (38),  and  they 
must  be  evaluated  numerically. 

If  the  driving  bunch  has  a  Gaussian  charge  distribution 


/.(z)  =  -7=—  exp 
s/2na 

then  the  bunch  wake  potential  Uz  becomes 
Ut(s)  =  I  2(z)Wz(s  -  z)  dz 


v'  2/n 
£0ti  da 


111 


£pCmnl  1  -(-lJPcOS^* 


4mn  ‘^mn 


(45) 


x 


«w(s  ~  Z)  , 
c 


(46) 


When  s  <r,  the  bunch  wake  potential  becomes 


<30  00  CO 


U'i5)XJVdZ  1  Zexp(-(^)^/2 

So71"  m  — 0  n  =  1  p-0 


(Drr 


£p^mn|  1  ~  (  —  l)PCOS 

9m*  ^mn 


C 

^mnp^ 


COS 


Wmnpf 


(47) 


From  this  equation,  we  can  see  that  contributions  from  the  modes  whose  resonant 
wavelengths  are  much  shorter  than  the  bunch  length  2 a  become  negligible.  For  the 
dominant  mode.  Eq.  (47)  can  be  conveniently  written  as 


U 


(i  -  gc2)1/2;.2c0l 

E0n2  del SN0l 


e 


-  II2(2mo/1)2 


(' 


(48) 


where  S  =  rcxhyb  is  the  cross  sectional  area  of  the  cavity,  A  is  the  wavelength  of 
dominant  mode,  and  ec  is  the  eccentricity  of  the  boundary  ellipse.  From  Eq.  (48),  we 
see  that  the  wake  potential  scales  as  cu*2  and  5“ '. 


3.2.  Transverse  Wake  Potential 


From  the  Panofsky-Wenzel  theorem,  the  transverse  wake  potential  is  related  to  the 
longitudinal  wake  potential  by 


3Wx(5) 

ds 


=  V,  Wz{s). 


(49) 


From  this  relation  we  can  write  the  transverse  delta  function  wake  potential  as7 


WJs)~£c 


^(QVx  Vk(t) 
2ulwk 


(50) 


On  the  accelerating  beam  path  (£  =  0,  q  =  n),  we  have 

=  CeJ 0,  qi)ceJ0,  qjcem(7 1,  qJV.CeJO,  qx) 


(51) 


Since  the  driving  charge  and  test  charge  are  assumed  to  pass  through  ea:h  focus 
axis,  only  £ -component  of  the  transverse  wake  potential  exists  at  the  foci  of  the 
elliptic  cavity.  Thus, 


V.CeJO,  qx)  =  \  lim 


1 


4—o  h  sinh  ?  dq 


z  CeJ  c,  qx) 


(liw. 

\n  r«0 


for  m  =  0, 2, 4, . . . 
for  m  =  1, 2, 5, . . . 


(52) 


Therefore,  the  transverse  delta  function  wake  potential  is  given  by 

<axd\ 


a 


5pc/l  (  1  )p  cos 
- - - sin 


a>iqiNx 


(?■)■ 


(53) 


where 


Cx  =  CeJO,  qma)cem( 0,  ?mn  )cem(tt,  ^mn  )V±Cem( 0,  <?mn) 

( I  ^S?)3( i(2r)2A<2”j  for m  =  0,2,4,... 

u  ( ‘J(  ^0(2r  +  I,2x*+  •)  for  m = 3’ 5’  •  • •  • • 


(54) 


For  the  driving  bunch  of  a  Gaussian  charge  distribution,  the  transverse  bunch 
wake  potential  becomes 


UJs)=  j  Mz)WJs  -  z)  dz 
J  —  00 


40  3D  00 


j-  j-  j 

t^naa  m«0  n*  1  p*o 


^mnp^7  mn  Nmn 


L"p( 


^  ^  •  ^mnp 


2*V 


]  sin - 2  (s  —  z)  dz. 

c 


(55) 


4.  NUMERICAL  EXAMPLES  AND  DISCUSSIONS 


The  longitudinal  and  transverse  wake  potentials  in  an  elliptic  pillbox  cavity  are 
calculated  by  using  Eqs.  (46)  and  (55).  We  choose  the  same  cavity  dimensions  and 
bunch  length  as  in  Ref.  9  to  compare  the  results.  The  cavity  dimensions  and  beam 
parameter  are: 

major  axis  2xb  =  10  cm, 
minor  axis  2 y„  =  6  cm, 
gap  distance  d  =  2  cm, 
a  =  5  mm. 

The  wavelength,  /.mnp,  and  the  loss  parameters  of  a  Gaussian  bunch,  k^(a),  for 
some  of  lower-order  resonant  modes  are  tabulated  in  Table  I,  while  in  Table  II  the 


TABLE  I 

Wavelength  and  Loss  Parameter  kmB^a)  of  an  Elliptical  Pillbox  Cavity  (xb  =  5  cm.  yb  =  3  cm, 

d  »  2  cm,  a  =  5  mm) 


m 

<Iml 

^rnlO 

(cm) 

'L.u 

(cm) 

(cm) 

(V/pC) 

n(ff) 

(V/pC) 

kmil(a) 

(V/pC) 

0 

1.7353 

9.5394 

3.6883 

1.9574 

3.2345  x  10“ 

1.6366  x  10“ 

6.8638  x  10“ 

1 

3.3522 

6.8634 

3.4559 

1.9201 

1.2629  x  10-1 

1.2993  x  10* 2 

5.7848  x  10“ 

2 

5.6530 

5.2853 

3.1895 

12821  x  10“ 

4.3203  x  10“ 

10914  x  10“ 

3 

8.6577 

4.2708 

19194 

1.8112 

15776  x  10“ 

8.4739  x  to*2 

4.5668  x  10- 3 

4 

12.3689 

3.5731 

16648 

1.7452 

10007  x  10" 1 

1.1262  x  10- 1 

6.9154  x  10“ 

5 

16.7792 

3.0678 

2.4343 

1.6754 

1.1060  x  10-' 

1.0854  x  10- 1 

7.7733  x  10- 3 

m 

1ml 

'•mlO 

(cm) 

'■mil 

(cm) 

'‘•mil 

(cm) 

km20(<r) 

(V/pC) 

kml  i(o) 

(V/pC) 

*»«(<») 

(V/pC) 

0 

11.3563 

3.7289 

17276 

1.7625 

18997  x  10'* 

1.4188  x  10-* 

8.4089  x  10' * 

I 

14.6278 

3.2856 

15389 

1.7084 

1.9755  x  10' 5 

1.4828  x  10“ 

9.8517  x  10“ 

2 

18.4878 

2.9225 

13598 

1.6505 

6.1309  x  10“ 

7.4392  x  lO'5 

5.6547  x  10“ 

3 

22.9665 

2.6221 

11930 

1.5902 

1.0984  x  10“ 

13849  x  lO"* 

11224  x  lO'3 

4 

28.0957 

2.3708 

10395 

1.5287 

1.1261  x  10“ 

5.4398  x  10“ 

5.8271  x  lO' 3 

5 

33.9196 

11577 

1.8990 

1.4668 

4.5322  x  t0'4 

9.1853  x  10“ 

1.2279  x  1C 5 

TABLE  II 


Wavelength  and  Loss  Parameter  k^Jia)  of  a  Circular  Pillbox  Cavity  (r0  ”  3.873  cm,  d  »  2  cm,  o  —  5 

mm) 


n 

Xm, 

(cm) 

/om 

(cm) 

''•O.J 

(cm) 

komdia) 

(V/pC) 

*0.l(<7) 

(V/pC) 

*o.i(<7) 

(V/pC) 

i 

1405 

10.118 

3.720 

1.962 

7.235  x  10* 1 

3.212  x  10“ 

1.338  x  10“ 

2 

5.520 

4.408 

1962 

1.821 

6.031  x  10* 1 

1.818  x  lO'’ 

9.614  x  10“ 

3 

8.654 

1812 

2.300 

1.630 

1.162  x  10“ 

1.700  x  lO" 1 

1.361  x  10“ 

4 

11.792 

1064 

1.834 

1.436 

4.419  x  10“ 

4.693  x  lO' 2 

7.128  x  10“ 

5 

14.931 

1.630 

1.509 

1.263 

3.944  x  lO" 3 

1710  x  10“ 

1.458  x  10“ 

same  information  is  given  for  the  circular  pillbox  cavity  that  has  the  same  cross 
sectional  area  (r0  =  v^y„). 

Figures  5  and  6  are,  respectively,  the  curves  of  the  longitudinal  wake  potential  on 
the  accelerated  beam  path  and  that  on  the  driving  beam  path,  in  which  different 
number  of  modes  are  included  (solid  lines  for  24  modes  and  broken  lines  for  12 
modes).  From  these  figures,  we  see  that  the  mode  summation  converges  rapidly, 
indicating  clearly  that  the  wake  potential  is  dominated  by  a  few  lower  modes.  It  was 
pointed  out  that  the  wake  potential  inside  the  bunch  is  difficult  to  calculate  because 
of  the  slow  convergence  of  mode  summation8.  However,  it  is  not  clear  in  these  figures 
whether  the  series  converges  rather  slowly  for  positions  inside  the  driving  bunch. 

On  the  accelerated  beam  path,  we  obtained  about  125  MeV/m/^C  acceleration 
gradient,  while  about  110  MeV/m /fxC  was  obtained  by  Y.  Chin9  (see  Fig.  7).  The 
elliptical  cavity  in  this  example  calculation  does  not  represent  the  maximum  accelera¬ 
tion  gradient  that  can  be  achieved.  As  discussed  earlier,  the  longitudinal  wake 
potential  is  proportional  to  the  number  of  particles  in  the  driving  bunch  and  inversely 
proportional  to  the  cross  sectional  area  of  the  cavity.  Also,  it  is  dependent  on  the 
distribution  of  charges  within  a  bunch  and  the  eccentricity  of  a  cavity.  By  choosing 
appropriate  parameters,  one  can  achieve  an  even-higher  acceleration  gradient. 

The  transverse  wake  potentials  on  the  accelerated  beam  path  and  driving  beam 
path  are  shown  in  Figs.  8  and  9  by  the  broken  lines.  The  corresponding  longitudinal 
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FIGURE  10  Plot  ol  normalized  wake  potentials  on  the  driving  beam  path  (results  of  Ref.  9). 

wake  potentials  are  plotted  together  in  order  to  see  if  we  can  find  the  positions,  such 
as  the  point  A  in  Fig.  8,  at  which  the  accelerating  potential  is  large  while  the  transverse 
potential  is  small.  The  charge  to  be  accelerated  should  be  positioned  at  such  a  point 
in  order  to  avoid  large  transverse  deflection  during  acceleration.  Figures  7  and  10 
are,  respectively,  the  curves  for  the  longitudinal  and  transverse  wake  potentials  on 
the  accelerated  beam  path  and  driving  beam  path  calculated  by  the  numerical  code 
WELL9.  In  these  calculations,  the  effect  of  beam  apertures  of  0.5  cm  radius  was 
considered.  Comparing  these  analytical  results  (Figures  8  and  9)  with  numerical 
method  (Figures  7  and  10),  we  found  very  good  agreement  in  both  magnitudes  and 
frequencies.  However,  magnitude  of  the  longitudinal  and  transverse  wake  potentials 
are  higher  for  this  analytical  method  which  does  not  include  the  beam  aperture  effects. 
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The  potential  of  the  wake-field  produced  by  a  bunch  of  relativistic  charged  particles  passing 
through  a  pill-box  cavity  is  expressed  by  using  Floquet’s  theorem,  and  an  obvious  requirement 
that  the  energy  gain  over  ill  acceleration  cavity  of  many  pill  boxes  must  be  proportional  to  the 
number  of  pill  boxes,  based  on  the  previous  modal  approach  ( BWW  theory).  It  is  found  that 
the  wake-field  is  consisted  of  two  classes  of  modes:  the  longitudinal  modes  which  are 
independent  of  the  aperture  and  the  pill-box  gap,  the  hybrid  (pill-box)  modes  which  are 
dependent  of  the  pill-box  gap.  The  wake  field  is  predominated  by  the  fundamental  longitudinal 
mode  whose  wavelength  is  on  the  order  of  the  effective  diameter  of  the  cavity,  and  its 
magnitude  is  inversely  proportional  to  the  cross  sectional  area  of  the  cavity  for  practical 
cavities  with  small  apertures.  Both  longitudinal  and  transverse  wake-fields  due  to  the 
longitudinal  modes  in  an  elliptical-pill  box  cavity  are  expressed  analytically  in  a  closed  series 
form  by  solving  exactly  the  longitudinal  eigenmode/equation  in  the  elliptical  cylindrical 
coordinates  in  terms  of  Mathieu  functions.  It  is  found  that  both  longitudinal  and  transverse 
wake-fields  whose  amplitudes  per  driving  charge  are  greater  than  100  MV/m/pC  can  be 
generated  in  an  elliptical  cavity. 


I.  INTRODUCTION  Fourier  series  based  on  the  scalar  and  vector  eigenmode 

functions  of  the  entire  cavity.  By  the  law  of  causality.  Bane  et 
i  Even  though  a  charged  particle  in  uniform  motion  on  a  al.1  expressed  the  part  of  the  Fourier  series  expanded  on  the 
j  straight  line  in  free  space  of  a  pill-box  cavity  does  not  radiate  scalar  eigenmode  functions  in  terms  of  the  part  on  the  vector 

!  (or  does  not  produce  “acceleration  fields”),  it  can  produce  eigenmode  functions.  Thus  the  wake  field  can  be  expressed 

j  “velocity  field,”  which  are  independent  of  acceleration.  in  a  Fourier  series  based  on  only  the  vector  eigenmode  func- 

These  velocity  fields  induce  currents  at  the  metallic  surface  tions  so  that  the  general  formulation  of  the  wake/neld  be- 

of  cylinder  and  disk  plates,  which  in  turn  produce  the  elec-  come  much  tractable.  The  synchronism  between  particle 

trie  and  magnetic  field,  commonly  called  the  wake  field,  be-  and  mode  was  used  to  formulate  implicitly  the  Fourier  coef- 

hind  the  trajectory  of  the  charged  particle.  ficients  in  a  Fourier  series  expansion  of  the  wake  field.7  So 

The  wake  fields  produced  by  a  bunch  of  relativistic  far,  from  the  synchronism,  only  the  acceleration  gradient  in 

charged  particles  in  pill-box  cavities  are  important  not  only  a  cavity  in  the  limit  of  vanishing  aperture  has  been  explicitly 

for  future  high-current  compact  electron  accelerators,1  but  calculated,  and  it  was  shown  explicitly  to  be  proportional  to 

also  for  a  use  as  an  electric  wiggler  for  a  proposed  free  elec-  the  inverse  square  of  the  pill-box  gap.8  This  inverse  square 

tron  laser  scheme.2  The  suitability  is  based  on  an  estimate  law  is  appropriate  when  the  gap  distance  p  is  much  longer 

that  the  wake  electric  field  exceeds  100  MV/m  per  pC  driv-  than  the  thickness  of  the  disks.  When  p— 0,  the  pill-box  cav¬ 
ing  bunch  charge,  and  the  wavelength  is  on  the  order  of  a  few  ity  is  same  as  a  plain  cylinder  cavity  with  inner  radius  equal 

centimeters.  The  principle  of  acceleration  by  the  wake  field  ,  to  the  aperture  radius  (i.e.,  the  inner  radius  of  the  disks), 

in  metallic  cavities  was  experimentally  verified.3  Also,  other  .  Thus,  the  wake  field  gradient  approaches  to  the  zero  asp— 0. 

wake  fields  in  either  plasma  medium4,5  or  a  dielectric  cavity0  j  In.  the  above  modal  analysis,  the  disks  are  assumed  to  be 

have  been  experimentally  investigated.  The  energy  gain  in  infinitely  thin  so  that  the  inverse  square  law  holds  for  any  p. 

the  recent  plasma  wake  experiment5  is  about  4  MeV  whereas  The  principle  of  the  modal  analysis  for  a  physical  quan- 

those  in  the  previous  wake  field  experiments3,4  6  are  on  the  '  tity  is  to  determine  the  Fourier  coefficients  based  on  the 
order  of  100  keV.  complete  basis  functions.  Since  the  basis  functions  are  com- 

In  the  rf  acceleration  cavity,  the  current  source  is  out-  plete,  the  Fourier  coefficients  are  uniquely  determined, 
side  the  ca .  ity,  and  the  phase  velocity  of  an  electromagnetic  Therefore,  once  the  Fourier  coefficients  are  determined  to 
wave  is  determined  entirely  by  the  cavity  structure.  In  con-  satisfy  some  proper  conditions  pertinent  to  the  physical 
trast,  the  wake  field  is  produced  by  the  driving  charge,  j  quantity,  there  are  no  other  proper  Fourier  coefficients  than 
Therefore,  the  phase  velocities  of  the  wake  field  modes  are.  1  these  coefficients.  In  this  article,  using  an  obvious  attribute 
determined  by  both  the  velocity  of  the  driving  charge  which  \  j  of  the  acceleration  gradient  in  the  periodic  cavity,  we  show 
produce  the  wake  field,  and  the  cavity  structure.  \  I  the  formulation  to  d  termine  the  Fourier  coefficients  for  any 

The  experiments3  showed  the  predominance  by  the  fun-  M  pill-box  cavity  of  arbitrary  aperture.  Using  the  formulation, 
damental  mode,  indicating  that  the  modal  analysis  i^appro-  |  the  wake-potential  in  a  cavity  of  vanishing  aperture  is  ana- 
priate  tool  for  the  evaluation  of  the^wake  field.  In  TO  pre-  [  lyrically  expressed,  and  the  result  is  exactly  the  same  as  the 
I  vious  modal  analysis,7  the  wake/field  is  expanded  in  a  result  from  the  synchronism. 


To  apply  the  modal  analysis,  we  first  find  the  basis  func¬ 
tions  (or  the  mode  functions).  By  using  Floquet’s  theorem 
on  the  periodicity,  the  problem  all  over  the  acceleration  cav-  1 
ity  can  be  reduced  to  that  in  a  unit  pill  box.  Since  the  mode  i 
functions  satisfy  the  elliptic  partial  differential  equations,  we 
should  specify  boundary  conditions  all  over  the  enclosing 
surface  to  determine  mode  functions  uniquely.  On  the  metal-  ! 
lie  surface,  the  infinite  conductivity  condition  can  be  as-  j 
signea  a:  the  pertinent  boundary  condition  as  usual.  How¬ 
ever,  so  fur  the  problem  as  to  what  is  the  boundary  condition 
on  apert  ire  has  not  been  investigated;  the  previous  modal 
|  analysis  s  applied  to  the  cavity  in  the  limit  of  vanishing  aper- 
I  ture,  wh  „h  does  not  have  such  problem,  and  the  numerical 
methods  which  solve  directly  the  Maxwell  equations  and  | 
boundaiy  conditions  all  over  the  entire  cavity  (not  the  unit 
pill-box  i  in  the  time-  and  space-mesh  domain,  e.g.,  the  TBCI  ' 
and  T3  code,  have  been  used  for  the  cavities  with  nonzero 
aperture.7,9  We  show  here  how  should  the  boundary  condi-  I 
tion  on  aperture  be  assigned  by  considering  the  parity  and  i 
the  continuity  at  the  aperture.  Even  though  the  modal  analy-  [ 
sis,  in  principle,  cannot  take  into  account  the  effect  of  large  1 
aperture,  the  specification  of  the  boundary  condition  is  still 
of  academic  interest.  For  large  apertures,  we  should  use  oth¬ 
er  methods. 10 

As  is  mentioned,  one  important  application  of  the  wake 
field  in  a  pill-box  cavity  other  than  the  acceleration  is  to  use 
the  wake  field  as  a  wiggler  field  for  a  free  electron  laser.  For 
this  application,  the  cavity  should  be  a  two-beam  configura¬ 
tion  such  as  an  elliptic  cavity2  or  an  annular  circular  cavity7 
so  that  the  lasing  electrons  can  be  injected  in  the  opposite 
direction  of  the  driving  charge  to  be  wiggled  by  the  wake- 
field  produced  by  the  driving  charge.  Being  an  axi-symmet- 
ric  configuration,  the  annular  cavity  cannot  have  the  electric 
field  in  the  transverse  direction  (i.e.,  the  direction  perpen¬ 
dicular  to  the  cavity  axis)  along  the  cavity  axis  (“transverse 
wake  field”).  Therefore,  as  far  as  free  electron  lasing  is  con¬ 
cerned,  the  annular  pill-box  cavity  is  similar  to  the  usual  rf 
acceleration  cavity  in  which  the  electric  field  undulating 
along  the  cavity  axis  can  be  routinely  made  more  than  10 
MeV/m.  The  rf  cavity  cannot  support  both  a  TM  mode  of ) 
the  same  wavelength'which  has  an  electric  field  in  the  trans-  • 
verse  direction.  Since  the  wake-field  is  produced  by  the  driv¬ 
ing  charge  moving  inside  the  cavity  as  is  already  mentioned, 
the  elliptic  wake-field  cavity  have  both  the  longitudinal  elec¬ 
tric  field  and  the  transverse  electric  field  of  the  same  wave¬ 
length.  If  we  inject  the  lasing  electrons  in  the  opposite  direc¬ 
tion  along  the  path  passing  through  the  foci  of  the  elliptic 
disks  on  the  other  side  of  the  foci  through  which  the  driving 
charge  is  passing,  the  lasing  electrons  are  forced  to  emit  the 
laser  radiation  by  both  the  longitudinal  wake-field  and  the 
transverse  wake-field.  Until  the  field-strength  of  laser  wave 
arrives  that  of  the  transverse  wake-field  (“weak  laser  re¬ 
gime”),  the  transverse  motions  of  the  electrons  are  still  con¬ 
fined  by  the  transverse  wake-field  field.  Similar  to  the  pres¬ 
ent  free  electron  lasing  in  the  magnetic  wiggler,  both 
stimulated  bremss,rahlung  and  unstimulated  bremsstrah-  j 
lung  by  the  transversely  bound  electrons  requires  the  trans-  ! 
verse  current  produced  by  the  transverse  wake-field  along, 
and  take  place  at  the  same  wavelengths.2  A  transverse  wake 
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field  whose  amplitude  and  wavelength  are  on  the  order  of 
100  MV/m  and  5  cm,  respectively,  is  equivalent  to  a  magnet¬ 
ic  wiggler  on  the  order  of  3  kG  and  5  cm,  respectively,  which 
are  about  the  standard  values  for  the  magnetic  wigglers  used 
in  the  present  free  electron  lasers.  Therefore,  even  if  we  just 
consider  the  transverse  wake-field,  the  elliptic  wake-field 
cavity  is  superior  to  the  conventional  magnetic  wiggler  if  we 
can  produce  the  transverse  wake  fi<  .Js  whose  amplitudes  are 
more  than  100  MV/m. 

If  we  inject  an  electron  beam  in  the  direction  of  the  lon¬ 
gitudinal  electric  wiggler,  the  transverse  motions  of  the  elec¬ 
trons  are  not  bounded  by  the  electric  field,  which  contrasts 
with  the  transverse  motions  of  the  electrons  in  the  free  elec¬ 
tron  laser  using  a  magnetic  wiggler.  Also,  the  longitudinal 
motions  of  the  electrons  are  not  bounded  by  the  longitudinal 
electric  wiggler  if  the  initial  kinetic  energies  are  sufficiently 
large  which  is  the  usual  case  in  any  free  electron  laser.  Free 
electron  lasing  by  the  electrons  whose  transverse  and  longi¬ 
tudinal  motions  are  not  bound  by  the  wiggler  does  not  re¬ 
quires  a  transverse  current  produced  by  the  wiggler  along.' 
The  transverse  motion  in  the  longitudinal  electric  wiggler  is 
primarily  governed  by  the  incident  electromagnetic  wave, 
and  the  incident  laser  wave  is  the  major  field  to  drive  the 
source  current  with  the  help  from  the  longitudinal  electric 
wave.  Unlike  transversely-bound  FELs,  stimulated  emission 
in  the  longitudinal  electric  wiggler  occurs  as  a  two-quantum 
Stark  emission1 1,12  in  which  the  electron  make  the  transi¬ 
tions  from  a  free  statd'viata'' vifiuai  free  state"  "The  two-quan¬ 
tum  Stark  emission  by  the  bound  electrons  are  used  as  a 
diagnostic  means  in  plasma  spectroscopy.  By  this  two-quan¬ 
tum  Stark  emission,  a  light  is  emitted  in  the  same  direction  as 
the  applied  electric  field  by  either  stimulated  emission  or 
spontaneous  emission.  It  had  been  theoretically  shown  that 
the  shorter  the  wavelength  of  the  free  electron  laser  which 
we  want  to  create,  the  more  advantageous  is  the  longitudinal 
electric  wiggler  compared  with  the  transverse  magnetic 
wiggler.2 

The  emission  in  the  so-called  Smith-Purcell  laser13  is  a 
spontaneous  emission  in  both  transverse  and  longitudinal 
wake  fields  generated  by  non-relativistic  electron  bunches 
when  they  pass  over  a  ruled  grating,  which  acts  similarly  as 
the  pill  box.  Since  the  electrons  are  non-relativistic  in  this 
laser,  they  can  experience  the  ripple  of  the  wake  field  which 
they  produce,  We  can  guess  that  the  electrons  experience 
•extraordinar^/large  electric  fields  since  the  gap  distance  is 
extremely  small  in  the  ruje|grating.  ■  ^  c  WhU  ><  t«  |'y  I  / 
Previously,  the  wake  fields  in  fin  elliptic  pill-box  cavity 
was  calculated  by  solving  directly  the  Maxwell  equations 
and  boundary  conditions  over  the  entire  cavity  in  the  time- 
and  space-mesh  points.9  The  advantage  of  the  direct  solution 
of  the  Maxwell  equations  is  that  the  aperture  effect  can  be 
taken  into  account.  However,  the  direct  solution  requires  a 
very  large  number  of  numerical  processing.  In  contrast  to 
the  direct  solution,  any  modal  analysis,  in  principle,  cannot 
take  into  account  properly  the  aperture  effect  when  the  aper¬ 
ture  is  large.  However,  the  experiment  showed  that  the  pre¬ 
dominance  of  the  wake  field  by  the  fundamental  mode  even 
for  a  cavity  of  considerably  large  aperture,3  indicating  that  a 
modal  analysis  can  render  a  reasonably  good  approximate 
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result®  unless  high  accuracy  is  demanded.  Both  longitudinal 
and  transverse  wake  fields  due  to  the  longitudinal  modes  can 
be  calculated  to  any  degree  of  precision  since  we  can  express 
these  fields  analytically  in  closed  series  forms.  This  expres¬ 
sion  is  possible  by  solving  analytically  the  longitudinal  eigen- 
mode^equation  in  the  elliptical  cylindrical  coordinates  in 
terms  of  Mathieu  functions14  as  will  be  shown  in  this  article. 


II.  THE  LONGITUDINAL  WAKE  POTENTIAL 


c*ew 


For  clarity,  let  us  consider  a  wake-field  consisting  of  a 
very  large  number  of  elliptic  pill  boxes.  A  point  charge  Q 
passes  through  om^fthe  foci  of  the  elliptic  disks  jyith-a- 
speed  c.  Let  the  traftwerie" coordinates~oTThese  centers  be 
rx  =  (xo,0)  and  (  —  xo,0).  The  delta-function  longitudinal 
wake  potential  Wt  is  defined  as  the  energy  gained  by  a  unit 
negative  point  charge  with  a  velocity  c  traveling  in  the  path 
passing  through  the  other  elliptic  foci,  at  a  distances  in  the  z 
direction  behind  the  accelerating  point  charge  Q. 

The  pill-box  cavity  is  a  periodic  structure,  the  wake  field 
in  this  cavity  should  be  very  amenable  to  a  modal  analysis. 
However,  since  the  wake  fie)  1  is  consisted  of  both  electro¬ 
magnetic  (vecior)  and  electrostatic  (scalar)  fields,  the  com1 
plete  basis  functions  for  the  Fourier  series  expansion  should 
include  both  scalar  and  vector  eigenmode  functions.  A  trac¬ 
table  analysis  of  the  Fourier  analysis  (or  modal  analysis) 
based  on  these  two  different  classes  of  basis  functions  had 
been  nearly  impossible.  This  difficulty  is  resolved  by  Bane  et 
alP  They  found  the  relationship  between  the  scalar  wake 
field  and  the  vector  wake  field,  which  is  arisen  by  the  causal¬ 
ity.  Accordingly,  imposing  the  causality  and  using  Max¬ 
well’s  equations,  they7  (BWW  theory)  first  obtained 


uai  (z  +  np)  =  exp(m/?a  )aai  (z),  (5) 

where /?„  is  an  arbitrary  constant.  By  means  of  Eq.  (5),  we 
have 

I  rL  >  rp  n„ 

aaI{.z)cxpUo>az/c)dz  =  aai(z)cxp(icoaz/c)dz  V 
I  Jo  Jo  „~“o 

Xexp[/(/Ja  +  coap/c)n],  (6) 

r  I 

where  n0  =  L/p  —  1. 

Even  though  any  value  of 0a  satisfies  the  geometric  con- 
)  dition  that  the  cavity  is  a  periodic  structure,  only  a  particular 
'  value  satisfies  other  physical  requirements.  An  obvious  re¬ 
quirement  is  that  the  energy  gain  over  the  entire  cavity, 

|  Wz  (s^s  proportional  to£  (or  the  number  of  the  pill  boxes). 

|  To  satisfy  this  requirement,  we  find  from  Eq.  (6)  that  /?„ 
must  be  given  by 

Pa=  -COaP/C.  (7) 

Substituting  Eq.  (7)  into  Eq.  (6)  and  combining  the 
j  resulting  equation  with  Eq.  ( 1 ),  we  derive  the  longitudinal 
wake  potential  (or  the  acceleration  gradient)  as 

W  ArrO 

‘.-J.  =  ZJLZ  (Jco,0./la )  Va  (  -  xo,0,Aa  )cos(s/A„ ) .  ( 8 ) 

I  L  P~  a 


where 


y a  ( x,yj.a )  =  f  <,(  x,yj )  exp  ( iz/Aa )  dz. 
Jo 


Here,  all  distances  such  as  z  and  Aa  are  measured  in  units  of 
the  pill-box  gap p,  and  a£  =  ca  a„  where  ca  is  a  constant  such 
that  !,  ' 


^(i)  =  4ir£?2  |J  0aiU)exp(/<ua2/c)<fc|  ^  (1) 

where  aa  ii  the  normalized  vector  eigenfunctions  given  by 
the  following  equation  and  boundary  conditions: 


V2aa  -1 — r  aa  =  0 
<r 


Vaa  =0  (3) 

everywhere,  and 

a„  X  n  =  0  (4) 

on  the  surface  of  the  metallic  enclosure  of  the  accelerator 
cavity  where  h  is  the  normal  vector  of  the  surface. 

From  this  point,  the  wake  potentials  is  formulated  by 
using  an  obvious  attribute  of  the  periodic  structure  which 
the  previous  BWW  theory7  has  not  used.  Also,  the  math¬ 
ematical  treatment  is  different  from  that  in  the  BWW  theo- 

From  this  point,  the  wake  potentials  is  formulated  by 
using  an  obvious  attribute  of  the  periodic  structure  which 
the  previous  BWW  theory7  has  not  used.  Also,  the  math¬ 
ematical  treatment  is  different  from  that  in  the  BWW  theo- 


Since  the  acceleration  cavity  is  a  repetitive  structure 
with  period  p,  au  (z)  is  required  by  Floquet’s  periodicity 
theorem  to  satisfy 


J  Jk  J <-*?dir=SaJi, 


where  Vp  denotes  the  space  in  the  pill-box.  Equation  ( 8 )  is  a 
similar  form  of  the  Coulomb  law  in  that  the  acceleration  is 
inversely  proportional  to  the  square  of  a  distance,  p. 

If  the  driving  inch  has  a  Gaussian  distribution 

g(z)  =  (2tt)  "  ,/2exp(  -  r/2cr)/cr, 
the  wake  potential  for  s$><7  is  given  by 

■— s  yz  )  vz  ( -  xt„oAa ) 

L  p-  a 

Xexp(  —  cr/2A  2)cos(5/Aa ).  (11) 

From  Eq.  ( 1 1 ),  we  find  that  the  modes  whose  wavelengths 
are  much  shorter  than  the  length  of  the  Gaussian  driving 
bunch,  2cr,  do  not  contribute  significantly  to  the  wake  field, 
and  the  shorter  the  wavelength  of  a  mode,  the  smaller  is  the 
contribution  to  the  wake-field  from  this  mode. 


\  III.  CLASSIFICATION  OF  WAKE  MODES  AND 
BOUNDARY  CONDITIONS 

1  For  brevity,  we  drop  the  superscript  p  which  denotes  the 
I  pill-box  eigenfunction  unless  there  is  a  chance  of  confusion, 
/  and  choose  the  origin  of  the  pill-box  coordinates  at  the  center 
of  the  pill-box.  Then,  the  pill-box  vector  eigenfunctions 
a  =  a,  +  azz  satisfy  the  following  equations  and  boundary 
conditions: 


+  at/A 2  =  0, 

(12) 

V2a,  +  a,/ A 1  =  0, 

(13) 

%+V.x-o 

dz 

(14) 

in  the  pill  box  j 

ax  =  0 

(15) 

on  the  disk  plates  at  2  ==  p/2, 

i 

aA  Xn  =  0, 

(16) 

ft 

H 

II 

O 

(17) 

/(rt)=0, 

on  the  periphery  of  the  cross  section,  and 


aSib-o, 


jik 


)dx  dy  =  1, 


on  the  periphery  of  the  cn  ss  section  (the  ride  surface  of  the 
cylinder). 

From  Eqs.  ( 14)  and  '  15),  we  should  impose 

da, 

—  =  0  (18) 

dz 

on  the  disk  plates  at  z  =  ±  p/2. 

We  must  impose  a  condition  that  a  and  its  derivatives 
are  continuous  functions  of  z  as  z  crosses  an  aperture  to  the 
next  adjacent  pill  box. 

Since  a  should  satisfy  Eqs.  (12)  and  (13),  and  the  pill 
box  is  symmetric  with  respect  to  the  reflection  about  the 
center  plane  (2  =  0)  of  the  pill  box,  a,  must  be  either  sym¬ 
metric  (even  parity  with  respect  to  2)  or  anti-symineiric 
(odd  parity  with  respect  to  2).  However,  because  of  Eq. 
(14),  a!  and  a ,  must  have  different  parities  (the  pair  of  ei¬ 
ther  a,  and  da^/dz  or  aL  and  da,/dz  should  have  the  same 
parity). 

We  need  only  to  prescribe  the  mode  functions  in  the 
region  0<2<p/2»»The  modes  can  be  classified  into  following 
classes.  k. 

(i)  The  odd  hybrid  modes:  These  modes  are  dependent 
on  the  pill-box  gap  and  aperture,  and  their  a,  have  the  odd 
parity  (consequently,,  Vfli/Ya  must  have  the  odd  parity). 
Accordingly,  \  \ 

d  \  v 


V,Jtr1)+^-  =  0, 
inside  in  the  cross  section 


where  S  is  the  cross  section  of  the  pill  box.  These  longitudinal 
modes  are  independent  on  the  apertures,  and  their  wave¬ 
lengths  are  determined  only  by  the  cross  section. 


•V.  WAKE  POTENTIAL  IN  A  CIRCULAR  PILL-BOX 
CAVITY 

The  modal  analysis  presumes  the  periodicity  of  the  ac¬ 
celeration  cavity,  and  thus  neglects  the  velocity  field  directly 
arrived  at  the  accelerated  charge  compared  to  the  electric 
and  magnetic  field  produced  by  the  currents,  which  are  in¬ 
duced  at  the  surface  of  the  disk  plates  and  the  cylinder  sur¬ 
face  by  the  velocity  field.  This  assumption  can  be  approxi¬ 
mately  met  when  the  aperture  is  much  smaller  compared  to 
the  cross  section  of  the  pill  box  for  practical  lengths  of  the 
cavity.  Therefore,  the  modal  analysis,  in  principle,  renders 
more  accurate  values  as  the  aperture  becomes  smaller,  and 
should  not  be  applied  to  the  limit  of  no  disks.  The  experimen¬ 
tal  results3  showed  that  the  wake-field  is  predominated  by 
the  fundamental  longitudinal  mode  even  when  the  ratio  of 
the  aperture  radius  (r0)  to  the  cavity  radius  (/?„)  is  consider¬ 
ably  large  ( rjR ~ ]  for  cavity  1  in  the  experiment).  There¬ 
fore,  we  can  take  the  wake-field  in  the  limit  of  vanishing 
aperture  as  a  good  approximation  of  the  wake  field  for  most 
practical  cavities  whose  aperture  radii  are  not  greater  than  \ 
of  the  cavity  radius. 

For  a  circular  pill  box  with  finite  aperture  radius, 
at  =  a‘r  and  Eq.  ( 16)  is  automatically  satisfied  for  this  pill 
box. 

Since  in  the  cylindrical  coordinate  system,  Eqs.  (12) 
and  ( 1 3 )  are  written  as 

d2a  1  dar  a  d2a  ar 

+ - - - _+  +  =0,  U4) 

dr  r  dr  r  dr  A' 


at  z  =  0  and  on  the  aperture  at  2  =  p/2. 

( ii )  The  even  hybrid  modes:  These  modes  are  dependent 
on  the  pill-box  gap,  but  independent  of  the  aperture.  Their  at 
and  da,/dz  have  the  odd  parity.  Equations  (15),  (18),  and 
the  continuity  of  the  odd  functions  at  the  aperture  are  com¬ 
bined  to 


at  7  =  0  and  2  =  p/2. 

(iii)  The  pure  longitudinal  (ax  =0)  modes:  The  wave¬ 
lengths  of  these  modes  are  completely  independent  of  both 
the  pill-box  gap  and  the  aperture.  From  the  fact  that  a!  =  0 
satisfies  Eq.  (12)  for  any  A,  Eq.  (15),  and  Eq.  (16), 
•>«,*  —A*i  )*  are  eigenmode  if 


V^i+-L:?+¥i+TT=0’  (25> 

dr  r  dr  dzr  A1 

and  a ,  and  a,  are  finite  -u  r  =  0,  we  have 

,  =  a,  =  0,  ^  =  0  (26) 

■  dr  dr 

at  r  =  0. 

Further,  Eqs.  (24)  and  ( 14)  are  combined  a  single  sec¬ 
ond-order  differential  equation  with  respect  to  2: 

£31  +  ^  =  ^..  (21) 

dr  A  ~  dz  dr 

For  convenience,  we  will  omit  the  superscripts  “even”  and 
"odd”,  and  the  subscript  “hyb”,  where  there  will  be  no  possi¬ 
ble  change  of  confusion. 


A.  Wake-field  gradient  from  all  even  hybrid  modes 

Solving  Eq.  (25)  with  boundary  conditions 
a,(Rotz)  =  0  [Eq.  ( 17)  ],  da,  (0j)/dr  =  0  [Eq.  (26)], and 


(34) 


<?a. (<•,())  _  da:(r,p/2)  _ 
dz  dz 

[Eq.  (20)],  we  obtain 

az(rj)  =AJ0(vrr)cos(v2z),  (28) 

where  J„  (g)  is  Bessel  function  of  order  n,  vr  —  £„/R0  with 
♦'eing  the  nth  root  of  J0{£)  =0,  vt  —  2 nir/p  with 
i.  =  1,2,3,...,  and/)  is  the  normalization  constant  to  be  deter¬ 
mined  later. 

Substituting  Eq.  (28)  into  Eq.  (27),  and  solving  the 
resulting  equation  with  boundary  condition  ar(r,Q)  =0, 
which  is  derived  from  a,  being  odd  for  z,  a,(Q>z)  =0, 
<9a,(0,z)/<9r  =  0  [Eq.  (26)  ],  and  a,  (r, p/2)  =  0  [Eq.  (20)], 
we  obtain 


a,(rj)  =  —  ^(v,/v,)/,(v,r)  sin(vJz).  (29) 


From  Eqs.  (28),  (29),  and  ( 10),  we  find  the  normaliza¬ 
tion  constant  as 


,*=  /I .  *■  -Ti- 

\  ft  •^O^l(Sn)  vn.l m 


(30) 


where  =  [  (£„/R0)-  +  (mir/p)2]w2. 

Combining  Eqs.  (8),  (9),  (28),  and  (30),  we  find  that 
the  wake  field  gradient  from  all  even  hybrid  modes  is  given 
by 


i 

y  n  —  \  m  =*  even 


m  >  0 


[  1  -  (  -  D^cosfv^lJcostv^j) 

JULHt 


(31) 


|  B.  Wake-field  gradient  from  all  longitudinal  modes 


1 


t 


The  calculation  of  the  longitudinal  modes  from  Eqs. 
(21 )  to  (23)  is  elementary.  We  find  that  the  wake  field  gra¬ 
dient  from  all  longitudinal  modes  is  given  by 


W, 


J.lont  (£)  _  8 Q  [ 

n* 


COS(v„oP)]CQS(v„0S) 

JUDSl 

(32) 


C.  Wake-field  gradient  from  all  odd  hybrid  modes 

For  ro^0,  [f^n^must  satisfy  Eq.  (25)  and  the  bound¬ 
ary  conditions 

a,(r,0)=0,  a,(Ruj)=  0, 
a,(r,p/l)  *0  for  0<r<r„, 
da,(r,p/2) 


dz 

da  A  CU) 

dr 


■■  0  for  ru<r<7?u 


=  0, 


(33) 


and  a,  should  be  determined  by  Eq.  (27)  and  the  boundary 
conditions 


daAr.O)  daAr,p/2) 

- - - =  0,  - - - =  0,  for  0 <r<r0, 

dz  dz 

aAr,p/ 2)=0  for  r0<r<R0, 

da  A  0^) 

a,(0,z)=,  — ^ - =  0, 

l  dr 


We  can  solve  these  equations  only  by  the  numerical  meth¬ 
ods.15  The  result  from  the  numerical  methods  is  not  impor¬ 
tant.  The  reason  is  as  follows.  First,  the  modal  analysis  does 
i  not  take  .nto  account  properly  the  aperture  effect  as  men- 
i  tioned  before.  Therefore,  the  wake  field  gradient  from  all 
odd  hybrid  modes  by  means  of  an  modal  analysis  is,  in  prin- 
i  ciple,  appropriate  only  for  small  apertures  where  an  analyti¬ 
cal  good  estimate  is  readily  available.  Second,  we  can  readily 
find  that  the  wake  field  gradient  from  all  odd  hybrid  modes 
from  the  above  equations  is  anyhow  negligible  compared  to 
that  from  all  longitudinal  modes  for  most  practical  cavities. 
In  the  limit  of  vanishing  aperture  (i.e.,  ru— 0),  the  boundary 
condition  at  z  =  p/2  is  now 


daAr,p/2) 

dz 


=  0,  aAr,p/2)  =  0. 


(35) 


Then,  the  analytical  solution  for  the  wake-potential  from  the 
odd  hybrid  modes  is  given  by 

_  16Q  A 

f  n~ 

*-•  P  n  —  1  m  —  odd 

m  >0 

x  [!  -  (1  —  )mcos(vnwp)]cos(y„m5) 

(36) 

In  the  limit  of  no  disks  (i.e.,  r0— /?0),  the  boundary  con¬ 
dition  at  z  =  p/2  is  written  as 


Mv/H-o.  *3^1 -0. 

dz 


(37) 


Then,  the  wake-field  gradient  from  all  odd  hybrid  modes  is 
given  by 


_  16Q  ^  ” 

r  n2  2*  2* 

"  y  n  =  1  m  mt  even 

m  >  O 

x  [!  -  (  -  l)mcos(v„t,„/>)]cos(v„„,s) 

JULHl 

(38) 

This  wake-field  gradient  in  the  limit  of  no  disks  has  no  ap¬ 
propriate  physical  meaning.  Only  the  magnitude  of  this  field 
can  be  used  only  a  crude  estimate  of  the  component  of  the 
wak^field  gradient  which  depends  on  the  aperture. 

Sincere  eigenmodes  of  the  odd  hybrid  modes  is  in  the 
range  between  the  eigenmode  in  the  limit  vanishing  aperture 
and  that  in  the  limit  of  no-disks,  we  find  the  largest  eigenval¬ 
ue  (/l)  of  the  odd  hybrid  modes  is  about  p/tr. 

The  wake/potential  from  all  modes  in  the  limit  of  van¬ 
ishing  aperture  is  given,  by  adding  Eqs.  (31),  (32),  and 
(36),  as 


W,(S)  _  8g  y  [1  -  (  -  1  )  mCOS(  VnimP)  3 

1  ^  P"  -« 

Xcos  (v„ms).  (39) 

!  This  equation  is  the  exactly  same  as  Eq.  ( 8 )  of  Ref.  8  which 
I  is  derived  from  the  synchronism  between  the  accelerated 
|  particle  and  the  wake-field  modes. 

:  From  the  above  equations,  we  find  that  the  largest 

;  eigenvalue  of  the  longitudinal  mode  is  2irRu/gt~2.6R0, 

|  while  the  largest  eigenvalue  of  the  hybrid  mcdes  is  about 
j  p/tr.  This  guarantees  that  the  wake-field  from  all  hybrid 
;  modes  are  small  compared  to  those  from  the  longitudinal 
modes  for  most  practical  cases.  Further,  if'Ro./<7>0.26,  we 
have 

exp(  -  2.89<r/RS)cos(2.4s/R0),  (40) 

L  R  o 

|  where  Q,  R0,  and  s  are  in  Gaussian  units. 

To  demonstrate  the  above  features  of  the  modal  analy¬ 
sis,  the  wake  potentials  in  a  cavity  of  R0  =  10,  r0  —  3,  a  =  4 
|  are  calculated.  In  Fig.  1,  both  the  wake  potential  from  the 
nine  most  dominant  longitudinal  modes  and  the  one  from 
;  the  nine  most  dominant  even  hybrid  modes,  expressed  by 
Eqs.  ( 25 )  and  ( 26)  are  drawn  as  the  functions  ofs.  The  wake 
potential  from  the  nine  odd  hybrid  modes^omputed  numeri¬ 
cally  by  solving  Eqs.  (13),  (24),  and  (25)  by  means  of  the 
finite  difference  method.  The  largest  eigenvalues  of  the  lon- 
j  gitudinal  modes,  even  hybrid  modes,  and  odd  hybrid  modes 
I  are  4.158,  0.159,  and  0.215,  respectively.  Therefore,  for  the 
i  driving-charge  length  of  a  =  4,  the  wake  potential  from  all 
j  hybrid  modes  are  practically  zero  compared  to  that  from  the 
I  longitudinal  modes  as  shown  in  Fig.  1.  The  wake  potential 


from  the  longitudinal  modes  clearly  shows  the  predomi- 
;  nance  by  the  fundamental  longitudinal  mode. 

In  Fig.  2,  the  theoretical  energy  increments  over  the  to- 
|  tal  acceleration  length  of  a  cavity  with  the  same  cavity  geom- 
I  etry  of  cavity  1  used  in  the  experiment3  are  drawn  for  both  a 
!  driving  charge  having  the  length  which  makes  the  theoreti- 
|  cal  curve  fit  best  the  experimental  one  (curve  a),  and  an- 
!  other  driving  charge  having  the  same  length  as  the  experi¬ 
mental  one  (curve  b).  The  contribution  from  all  hybrid 
i  modes  is  not  considered  since  it  is  theoretically  negligible, 
j  Even  though  we  do  not  see  so  remarkable  predominance  of 
j  the  fundamental  mode  in  curve  b  as  in  the  experimental  re- 
i  suit.  Therefore,  we  can  conjecture  that  either  the  bunch 
|  length  was  not  measured  accurately  in  the  experiment,  or  the 
I  aperture  effect,  which  will  be  considerably  large  in  the  ex- 
j  perimental  case  where  r^/Ru=.\,  suppresses  more  strongly 
|  high-modes  than  the  fundamental  mode  or  the  finite  con¬ 
ductivity.  The  latter  reason  is  more  probable  since  the  result 
J  from  the  TBCI  code,  which  can  calculate  accurately  the 
wake  field  for  any  aperture,  is  in  excellent  agreement  with 
j  the  experimental  data. 

i 

!  V.  WAKE  FIELDS  IN  AN  ELLIPTIC  PILL-BOX  CAVITY 
|  A.  The  formula  of  the  transverse  wake  potential 

Consider  again  the  exciting  charge  Q  traveling  the  one 
i  foci  of  the  elliptic  cavity  at  v  =  c.  The  delta-function  trans- 
j  verse  wake  potential  Ux  is  defined  as  the  transverse  momen¬ 
tum  kick  experienced  by  a  test  charge  following  at  a  distance 
;  s  in  the  z  direction  on  the  other  foci  and  also  at  u  =  c.  Bane  et 
j  al.1  derived 


FIG.  1 .  Wake  potentials  in  unit  of  'iQ  /p'  as 
the  functions  of  the  distance  s  between  the 
driving  charge  and  the  accelerated  charge  in 
a  cavity  R„  =  10,  r„  =  3,  a  =  4.  Here,  p,  R„ 
and  2<r  are  the  cavity  radius,  the  aperture 
radius,  and  the  length  of  the  driving  charge 
in  unit  of  the  pill-box  gap  p,  and  Q  is  the 
driving  charge,  (a)  The  wake  potential  from 
the  fifteen  longest  longitudinal  modes,  (b) 
The  wake  potential  from  the  odd  nine  long¬ 
est  modes,  (c)  The  wake  potential  from  the 
even  nine  (m  =  2,4,6;  n  =  1,2,3)  longest 
hybrid  modes. 


WAKE  POTENTIAL 


FIG.  2.  Energy  increment  (in  unit  ofKeV)  over 
the  t  Hal  acceleration  length  as  a  function  of  the 
dista  ice  s  between  the  accelerated  charge  and 
the  i  riving  charge  in  unit  of  picosecond.  The 
cavi'  '  has  the  following  parameters:  cavity  radi¬ 
us  R,  =  19.5  mm,  aperture  radius  ru  =  6.3,  pill¬ 
box  ,ap  p  =  6.3  mm,  driving  charge  Q=  3.1 
nC,  number  of  pill  boxes  =  60,  and  rms  length 
of  tiie  driving  charge  <7  =  8, 4  mm.  The  contri¬ 
butions  from  the  fifteen  most  dominant  longitu- 
din  il  modes  are  considered.  ( a )  <7  =  8  mm.  ( b ) 
o  =  4  mm. 


=  <*o.CU„  )V,Fa(-  Xo,(Ua  )Xa 

Lt  pa 

Xsin(s/Aa),  (41) 

Accordingly,  the  transverse  wake  potential  driven  by  the 
driving  charge  of  an  arbitrary  charge  distribution  is  defined 
by 

U1(s)=J  Ui(s  +  s')gU')<fr\  (42) 

where  g(s)  is  the  charge  distribution  of  the  driving  bunch. 

B.  The  wake  potentials  in  an  elliptical  cylindrical  cavity 

From  the  previous  theoretical  result  ( Fig.  1 )  and  an 
experiment,2  we  found  that  the  aperture  of  practical  size 
does  not  affect  significance  wake  field,  and  the  wake  field 
is  predominated  by  the  longitudinal  modes.  Therefore,  we 
neglect  the  hybrid  modes,  and  consider  only  the  wake  fields 
contributed  from  the  longitudinal  modes. 

To  obtain  the  wake  fields  from  the  longitudinal  modes  in 
an  elliptical  cavity  whose  periphery  of  the  cross  section  is 
given  by 


I  we  must  solve  Eq.  (21 )  with  the  boundary  condition  such 
|  that /=  0  on  this  periphery.  We  first  note  that  the  solution 
j  must  be  an  even  or  odd  function  for y  since  both  the  equation 
[  Eq.  (21)]  and  the  boundary  condition  are  invariant  when.y 
i  is  replaced  by  —  y.  However,  the  odd  functions  do  not  con- 
1  tribute  to  either  longitudinal  or  transverse  wake-potential  at 
|  the  foci  whose  y  values  are  zero,  and  accordingly  should  be 
!  excepted.  Also,  we  note  that  /and  V  J  must  be  finite  at  the 


foci  since  both  the  longitudinal  and  transverse  wake-poten¬ 
tials  are  finite. 

The  solution  can  be  found  in  terms  of  known  functions 
by  the  transformation  of  the  coordinates  system  to  the  ellip¬ 
tical  cylindrical  coordinates.14 

The  elliptical  cylindrical  coordinates  (g,rj)  are  related 
with  the  cartesian  coordinates  by  the  following  equations: 

x  =  h  cosh  £  cos  77, 

y  =  h  sinh  £  sin  77,  (44) 

where  2 h  is  the  distance  between  two  foci,  i.e.,  h  =  y]x2h  —y\. 
Then,  the  problem  is  transformed  to  solve 

~~r  +  +  2k  2(cosh  —  cos  2ij)f=  0.  (45) 

<?£-  dy- 

in  the  region  0<£<£o,  0<t7<2n-with  the  boundary  condition 

:  Mo>V)  =  0,  (46) 

where  k  —  h  /IX  and  =  cosh  ~  1  (xb/^x2b  —yi ) . 

By  means  of  the  separation  of  variables,  we  have 
A£>V)  =^(^1^(47)  where  </>(.£)  and  ^(77)  satisfy  the  fol¬ 
lowing  Mathieu  equations: 

(a-2*2cosh2£)V>  =  0,  (47) 

—r^r  +  (a  —  2k2  cos 277)^  =  0,  (48) 

a 77* 

where  a  is  the  separation  constant. 

Since  we  consider  only  the  even  functions  of y,  we  have 
i  that  ^(77)  =  t  —  77).  From  this  parity  consideration,/ 
being  finite  at  the  foci  and  Eqs.  (46)— (48 ),  the  solutions  are 
expressed  as 


Ai.v)  =  '£'2  N"«,Cem  (£,q)cem  (7),q) 

in  tf 

+  1.'2N  (§,q)cem  ( y,q ),  (49) 

m  n 

where  q  —  k  2,  Nm<l  and  N  'mq  are  the  normalization  constant, 
and 

Csm(i>Q  >q)  —  Se,„(£u>q)  =0,  (50) 

for  any  value  of  q. 

Since  df/dx  =  VJ-xat  the  foci,  i.e.,  (£17)  =  (0,0)  and 
(0,tt),  is  given  by 

dx  h  i-u  \sinh£  dc,  J 

m?i=-±  ta(_!_  (51) 

dx  h  s-o\sinh£  dg  ) 

the  condition  that  the  transverse  wake-field  is  finite  leads  to 
the  discard  of  all  Sem  (£,$)cem  (77,$)  terms  in  Eq.  (49). 
Further,  we  are  concerned  with  large  eigenvalue  A,  which 
corresponds  to  small  values  of  q  under  the  condition  Eq. 
(50).  Accordingly,  we  can  ignore  the  terms  with  m>l. 
Therefore,  we  obtain 

Ai.V)  =  X  NiCe» (i>q)ceu(.T],q),  (52) 

where  q' s  are  the  real  roots  of  the  following  equation: 
0=1—}$  cosh  2|„  +  ^q2  cosh  4£0 

-  tis  9* (A  cosh  65o  -  7  cosh  2£u) 

+  Tjhu  94(cosh  8|o  -  320  cosh  4£0)  +  0($3), 

(53) 

and  TV,  is  the  normalization  constant  given  by 


-v=-  r  *  r 

N\  2  Jo  Jo 


Xd77(cosh  2£  -  cos  2 17)  [Ce0(£,$)ce0(77,$)  ]2.  (54) 

If  $’s'  and  TV,’ s  are  determined,  then  the  longitudinal 
wake  potentials  on  the  same  path  as  that  of  the  driving  bunch 
by  a  gaussian  driving  charge  of  length  a  is,  for  s><7, 

H'(O.O)  4 irQh'-^Nl  4,  %  P/tx 

- - - =  — ~y—a\q)sm-(ylq/h) 

L  p-  u  q  .  ,17-1-, 


Xexp(  —  2$trV/t2)cos(2y$s//i),  (55) 


VA0.0)  2 irQh'-^Nl  „  x  .  „  r,l% 

- ; - = - —  X  — r  a  (<?)0($)sm-(V$//») 

L  P-  *  Jq 


X  exp  (  —  2qcr/h 2 )  sin  ( 2 \fqs/h ) ,  ( 57 ) 


UAOfi) 


and  the  transverse  wake  potential  on  the  opposite  path  is 

Ux(0,ir)  2nQh2^N]  x0/  s  .  ,,r#l4 

- - - = - Y  —aHq)P(q)sin-Uq/h) 

L  p-  u  .fa 


Xexp(  —  2qcr/h  2)sin(2\[qs/h),  (59) 


UA0,ir) 


fi(q)  =  —  lim  (-  dCe°  f.^) 
q  s-o  \  smht  dt  J 


=  -2  +  }$  +  —  + 

From  hqs.  (43)  to  (48),  we  find  that  the  wake  potential 
is  symmetrical  with  respect  to  x  =  0  plane  in  the  lowest 
(m  =  0)  modal  analysis. 

Figure  3  shows  the  longitudinal  wake  potentials  com¬ 
puted  by  the  analytical  formulation  given  by  Eq.  (55)  for 
cavities  of  Q  —  1  nC,p  =  4mm,  and  xb/p  =  10  (i.e.,xfc  =  4 
cm),  a  =  8  mm,  and yb/xb  =  0.9, 0.5, 0.3  as  the  eccentricity 
parameter.  From  the  curves,  we  find  that  the  dominant 
wavelength  increases  as  yb/xb  increases,  that  is,  the  effective 
radius  given  by  R0  =  Jxbyb  increases.  We  also  find  that  the 
amplitude  of  the  longitudinal  wake  potential  does  not  be¬ 
come  maximum  as yb/xb  -»  1 ,  indicating  that  an  elliptic  cross 
section  has  a  larger  acceleration  gradient  than  the  circular 
cross  section.  The  acceleration  gradient  can  be  greater  than 
100  MeV/m  for  a  cavity  having  total  driving  charge  of  1  /iC, 
pill-box  gap  of  4  mm,  and  an  elliptical  cross  section  of  xb 
=  10  cm  and  yb  —  5  cm. 

Figure  4  shows  the  transverse  wake  potentials  for  the 
same  cavities  as  in  Fig.  3.  From  the  curves,  we  find  that  we 
can  make  the  undulating  electric  field  in  the  transverse  di¬ 
rection  whose  amplitude  and  wavelength  are  greater  than 
100  MeV/m  and  on  the  order  of  10  cm,  respectively,  with  a 
cayity  having  total  driving  charge  of  1  fiC,  pill-box  gap  of  4 
mm,  and  an  elliptical  cross  section  of  xb  =  10  cm  and  yb  =  5 


and  the  longitudinal  wake  potential  on  the  opposite  path  is 

r,(0,tr)  4 irQh2^N\  4,  %  .  „  P/Lx 

- - - =  — y  a  ($)sin*(V9/A) 

L  p-  T  a 


Xexp(  —  IqcP/h  2)co9(2\[qs/h). 


where 


a(q)  =  CV„(0,$)  =  ceo(0,$)  =  ce„(jr,$) 

Similarly,  the  transverse  wake-potential  on  the  same 
path  as  that  of  the  driving  charge  is 


|  Vi.  CONCLUSION  AND  DISCUSSIONS 

The  potential  of  the  wake-field  produced  by  a  bunch  of 
relativistic  charged  particles  passing  through  a  pill-box  cav- 
j  ity  is  expressed  as  a  Fourier  series  based  on  the  vector  eigen- 
|  mode-functions  of  the  unit  pill  box.  The  Fourier  coefficients 
are  uniquely  determined  by  Floquet’s  theorem,  and  an  ob- 
i  vious  requirement  that  the  energy  gain  over  all  acceleration 
cavity  of  many  pill  boxes  must  be  proportional  to  the  number 
!  of  pill  boxes.  The  eigenmode  functions  must  satisfy  the  ellip- 
tic^type;  partial  differential  equations  derived  from  Max¬ 
well’s  equations.  The  complete  boundary  conditions  on  the 


FIG.  3.  Longitudinal  w;  <e  potentials  Tor 
three  values  of  yb/xb  in  u  t  elliptical  cavity 
having  Q=  1  nC.p  =  4  mir,  andxb/p  =  10 
(i.e.,  xk=4  cm),  c  -8  mm.  (a) 
y>/xb=  0.9.  (b)  yjx„  =0.5.  (c)  yjxb 
=  0.3. 


whole  surface  of  the  unit  pill  box,  which  are  needed  to  deter¬ 
mine  uniquely  the  eigenmode  functions,  are  assigned  by  tak¬ 
ing  into  account  the  symmetricity  of  the  equation  and  ! 
boundary  conditions,  the  continuity  on  the  aperture,  and  the 
zero  divergence  of  the  eigenmode  function,  and  the  infinite 
conductivity  on  the  conducting  surface. 

From  the  above  formulation,  we  derive  the  wake  field  in 
a  similar  form  of  the  Coulomb  law.  In  particular,  the  wake 
potential  of  the  circular  cavity  in  the  limit  of  vanishing  aper¬ 
ture  derived  by  our  formulation  is  identical  with  that  derived 
by  the  usual  method  using  the  synchronism  between  particle 


and  wave  mode.  It  is  found  that  the  wake  field  consists  of 
three  classes  of  modes:  the  longitudinal  modes  which  are 
!  independent  of  the  aperture  and  the  pill-box  gap,  the  even 
hybrid  modes  which  are  independent  of  the  aperture,  but 
dependent  of  the  pill-box  gap,  and  the  odd  hybrid  modes 
which  are  dependent  on  both  the  aperture  and  pill-box  gap. 
The  contributions  from  both  the  hybrid  modes  can  be  ne¬ 
glected  compared^those  from  the  longitudinal  modes.  We 
have  found  that  the  wake/fbtential  in  most  practical  cavities 
)  is  predominated  by  the  contribution  from  the  fundamental 
|  longitudinal  mode  whose  wavelength  is  2.6  times  of  the  di- 


FIG.  4.  Transverse  wake  potentials  for  three 
values  of  yj xb  in  the  same  cavity  as  in  Fig.  3. 
(a)  y„/xb  =0.9.  (b)  y„/xb  =0.5.  (c)  yb/xb 
=  0.3. 
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ameter  of  the  pill  box,  and  is  inversely  proportional  to  the 
cross  sectional  area  of  the  cavity.  The  result  of  this  modal 
analysis  is  in  excellent  agreement  with  an  experiment  except 
that  the  length  of  the  driving  charge  differs  by  a  factor  2.  We 
can  attribute  the  difference  to  the  aperture  effect,  which  can¬ 
not  be  taken  into  account  properly  by  the  modal  analysis 
when  the  aperture  is  large  as  in  the  experiment. 

Unlike  the  circular  cavity,  we  can  generate  the  trans¬ 
verse  wake  field  in  the  elliptical  cavity.  Botlithese  fields  can 
serve  as  a  catalyzing  field  for  stimulated  bremmstrahlung, 
and  thus  are  important  in  the  application  for  a  free  electron 
laser.  Both  the  wake  fields  from  the  longitudinal  modes  have 
been  analytically  obtained  by-thw  nseof  .the  Mathiej#  func¬ 
tions.  The  analytical  result  shows  ma?  we* can  generate  both 
longitudinal  and  transverse  wake  fields  having  amplitude 
greater  than  100  MV/m  and  wavelength  on  the  order  of  10 
cm  if  we  can  produce  a  driving  £harg£  bunch  qfy  the  order  of 
1  fiC  in  an  elliptical  cavity.  ^ 
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ABSTRACT 


WAKE-FIELD  ACCELERATION  OF  CHARGED  PARTICLES 
WITH  AN  ELLIPTICAL  CAVITY 

Publication  No. _ 

Jin-Seok  Yang,  Ph.D. 

The  University  of  Texas  at  Arlington,  1990 
Supervising  professor  :  K.  Wendell  Chen 

Acceleration  of  charged  particles  by  the  wake  fields  excited  in  the  el¬ 
liptical  cavity  by  an  intense,  short  electron  bunch  is  investigated.  Two 
cavities  are  considered  as  a  wake-field  device:  disk-loaded  metallic  cavity  and 
dielectric-loaded  cavity.  The  analytical  methods  are  developed  to  evaluate  the 
wake  fields  and  wake  potentials  in  these  cavities.  Using  modal  analysis 
method,  it  is  shown  that  the  longitudinal  and  transverse  wake  potentials  in 
the  disk-loaded  elliptical  cavity  can  be  expressed  analytically  in  terms  of  the 
Mathieu  functions.  The  aperture  effects  are  not  considered  in  this  calcula¬ 
tion.  The  longitudinal  wake  potential  in  a  circular  cavity  are  derived  from 
that  in  an  elliptical  cavity,  and  it  is  shown  to  be  exactly  the  same  formula 
that  was  derived  earlier.  These  analytical  results  are  compared  with  numeri¬ 
cal  ones,  and  excellent  agreements  in  both  frequencies  and  amplitudes  are 
observed  between  these  two  methods.  Also,  the  wake  fields  in  the  dielectric- 
loaded  structures  are  derived  analytically  using  the  Fourier  transform  method. 
It  is  shown  that  the  dielectric-loaded  structure  can  support  strong  accelera- 


v 


tion  field,  and  has  a  useful  property  that  the  transverse  wake  fields  can  be 
made  quite  small  for  the  ultra-relativistic  driving  beam.  It  is  demonstrated 
that  acceleration  gradient  greater  than  a  few  hundred  MeV/m/uC  can  be 
achieved  in  the  disk-loaded  elliptical  cavity  or  in  the  dielectric-loaded  ellip¬ 
tical  cavity. 
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CHAPTER  I 


INTRODUCTION 

The  invention  and  continued  development  of  particle  accelerators  and 
the  associated  technologies  have  had  a  profound  impact  on  many  fields  of 
pure  and  applied  science  and  on  variety  of  branches  of  modern  technology. 
In  the  fields  of  high  energy  and  nuclear  physics,  a  marvelous  complement 
of  accelerator  facilities  exists  around  the  world.  In  other  area  of  science, 
synchrotron  light  sources  and  accelerator-driven  pulsed  neutron  sources  have 
opened  up  revolutionary  new  research  opportunities  in  materials,  chemistry, 
and  biological  research.  In  industry  and  medicine  there  are  literally 
thousands  of  accelerators  in  use  in  health  care  treatment,  radiation  steriliza¬ 
tion,  radiation  processing,  ion  implantation,  microchip  production,  etc. 

During  the  past  few  decades,  accelerator  energies  have  been  increasing 
from  KeV  to  TeV  at  the  rate  of  an  order  of  magnitude  in  about  every 
seven  years  due  mainly  to  imaginative  advances  in  accelerator  technology 
and  design.  As  the  demand  of  higher  accelerator  energy  is  increased,  it  be¬ 
comes  more  important  to  produce  higher  acceleration  gradient  to  decrease 
the  size,  cost  and  complexity  of  accelerators,  especially  for  the  future  linear 
colliders.  This  has  motivated  the  development  of  new  ideas  for  charged  par¬ 
ticle  acceleration,  which  are  expected  to  produce  acceleration  gradient  much 
higher  than  that  conventional  accelerators  can  produce.  Typical  acceleration 
gradient  of  conventional  accelerators  is  about  20  MeV/m.  The  new  accelera¬ 
tion  schemes  under  active  theoretical  and  experimental  research  are:  (1) 
Laser-Plasma  Acceleration,1'5  (2)  Plasma  Wake-field  Acceleration,6*9  (3)  Wake- 
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field  Acceleration,10"16  etc. 

Over  the  last  few  years,  substantial  progress  both  theoretically  and  ex¬ 
perimentally  has  been  made  on  these  new  acceleration  schemes.  The  plas¬ 
ma-based  acceleration  schemes  utilize  the  extremely  large  high  acceleration 
gradient  associated  with  plasma  waves  excited  by  a  number  of  ways,  for  ex¬ 
ample,  by  using  high  power  laser  beam1,2  (Laser  Wake-field  Acceleration), 
beating  of  two  relatively  lower  power  laser  beams3  (Plasma  Beat  Wave  Ac¬ 
celeration),  or  using  relativistic  electron  beams6  (Plasma  Wake-field  Accelera¬ 
tion).  For  those  schemes  involving  plasmas,  acceleration  gradients  of  the 
order  of  1  GeV/m  are  theoretically  possible.  But  plasma  devices  are  very 
complicated  and  their  practical  applications  have  a  number  of  unresolved  is¬ 
sues  which  have  not  been  fully  explored.  These  are:  fine  tuning  of  plasma 
density  and  laser  frequency,  production  of  driving  beam  with  a  slow  rise 
time  and  very  rapid  fall  time,  laser-plasma  instability,  focusing  of  laser 
beam,  etc  The  plasma  accelerators  really  exhibit  advantages  over  other  con¬ 
cepts  only  at  frequencies  greater  than  about  100  GHz,  and  the  problem  of 
producing  a  suitable  driver  beam  which  is  capable  of  propagating  a  suffi¬ 
ciently  long  distance  within  the  plasma  becomes  more  difficult  as  the  fre¬ 
quency  is  increased. 

In  the  moderate  acceleration  gradients  (several  hundred  MeV/m),  a 
class  of  wake-field  acceleration  schemes  based  on  metallic  cavity10"14  or 
dielectric-loaded  structure15,16  (or  combination  of  both)  are  promising  new 
concepts  in  view  of  their  structual  simplicity  and  ease  of  phase  matching. 
In  these  wake-field  schemes,  an  intense,  low-energy  beam  excites  the 
electromagnetic  fields  (so  called  "wake  fields")  in  a  metallic  cavity  or  in  a 
medium,  which  are  then  used  to  accelerate  a  second,  less  intence  beam  to 
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high  energies.  The  periodic  structure  plays  the  necessary  role  to  produce 
the  wake  fields  and  to  keep  constant  phase  between  an  accelerated  beam 
and  the  wake  fields  in  a  metallic  cavity,  while  in  the  dielectric-loaded  struc¬ 
ture  a  driving  beam  excites  the  slow  electrom  rgnetic  waves  through  the 
Cerenkov  radiation  mechanism. 

Let  us  briefly  discuss  some  features  of  early  works  on  these  wake-field 
acceleration  schemes.  Recently,  proof-of-principle  tests  for  the  wake-field  ac¬ 
celeration  schemes14,16  have  been  demonstrated  at  the  Argonne  National 
Laboratory  (ANL).  It  was  shown  that  a  few  hundred  MeV/m  acceleration 
gradient  can  be  achieved  in  the  disk-loaded  metallic  cavity  and  in  the 
dielectric-loaded  cavity  provided  the  driving  beam  is  sufficiently  intense  (of 
the  order  of  1  micro-Coulomb).  Wake-field  cavities  tested  have  a  common 
geometry,  i.e.,  both  the  driving  and  accelerated  bunches  pass  through  the 
same  path  (so  called  "co-linear  geometry").  For  the  co-linear  structures,  there 
is  a  fundamental  theorem  tv  *t  the  transformer  ratio,  defined  as  the  ratio 
of  the  maximum  energy  gain  per  particle  of  the  accelerated  bunch  to  the 
maximum  energy  loss  per  particle  of  the  driving  bunch,  can  not  exceed  two 
when  the  driving  bunch  has  a  symmetric  charge  distribution.  This  theorem 
complicates  the  practical  applications  of  wake-field  acceleration  schemes, 
especially  for  high  energy  linear  colliders. 

However,  we  can  overcome  this  problem  by  either  increasing  the  trans¬ 
former  ratio  or  staging  the  cavities.  The  transformer  ratio  can  be  increased 
by  a  number  of  ways,  for  example,  by  mixing,17  injecting  an  asymmetrical 
driver  bunch,18  injecting  multiple  bunches,13  non-linear  effects  in  dielectric, 
using  cavity  of  annular  geometry,10  etc.  But  practical  realization  of  such 
methods  meets  with  some  difficulties  and  limitations,  and  has  not  been 
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developed  yet.  Other  method  to  overcome  the  low  transformer  ratio  is  to 
use  multi-stage  scheme,  i.e.,  acceleration  in  short  distance.  For  accelerators 
of  low  transformer  ratio,  it  is  inevitable  to  replenish  the  energy  of  the  driv¬ 
ing  bunch  or  to  replace  used  driving  bunch  by  a  new  one  frequently.  For 
the  co-linear  cavities,  re-acceleration  of  the  decelerated  driver  or  replacement 
of  the  driving  charge  appears  difficult  and  technically  unattractive,  as  it  re¬ 
quires  disposing  of  the  spent  driver  beam  and  introducing  new  one  without 
disturbing  the  accelerated  beam.  Multi-staging  can  be  readily  realized  with 
non-concentric  cavities  which  have  separate  paths  for  the  driving  and  ac¬ 
celerated  beams. 

Among  the  co-linear  wake-field  devices,  the  cavity  of  an  annular 
geometry  had  been  once  rigorously  investigated.10  Here,  we  briefly  discuss 
the  features  of  this  well-known  geometiy.  The  largest  transformer  (impedance 
transformation)  ratio  is  expected  in  a  concentric  arrangement  as  shown  in 
figure  1.1.  In  this  scheme,  high  current  ring-shaped  driver  bunch  is  used  to 
accelerate  a  low  current  beam  following  the  central  axis.  For  a  symmetric 
driving  bunch,  transformer  ratio  in  excess  of  20  can  be  realized.  Utilizing 
asymmetrical  shaping  of  the  driving  beam,  it  has  been  shown  numerically 
that  an  optimal  transformer  ratio  as  high  as  100  is  possible  in  this  annular 
wake-field  device.  However,  the  practical  realization  of  this  device  or  the 
asymmetrical  shaping  of  the  driver  beam  meets  with  rather  stringent  ex¬ 
perimental  limitations.  While  the  production  of  short,  intense  annular 
electron  ring  beam  is  already  difficult,  but  feasible  (e.g.  workers  at  DESY, 
West  Germany,  once  tried  such  an  experiment),  the  additional  requirement 
that  the  electron  ring  beams  must  possess  no  rotational  movement  while 
traversing  the  series  of  annular  slots  is  exceedingly  difficult  because  the 
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Figure  1.1.  Sketch  of  an  annular  wake-field  cavity 


inner  disks  must  be  supported  in  place  with  wires  or  rods.  The  ring  beam 
will  be  lost  if  rotation  is  present.  In  addition,  even  if  the  driver  beam  does 
not  rotate  appreciably,  the  supports  destroy  the  apparent  rotational  symmetry 
by  creating  holes  in  the  ring,  thus  make  the  ring  beam  much  vulnerable 
to  possible  transverse  beam  instabilities. 

The  past  experimental  program  at  DBSY  is  not  adversely  affected  by 
the  these  difficulties  because  of  its  very  short  accelerating  structures  of  40 
cm.  For  a  long  accelerator  structure,  say  even  for  a  1  GeV  machine,  these 
problems  could  become  quite  troublesome.  Staging  could  help,  but  not  so 
attractive  due  to  its  concentric  geometry.  Replacing  the  annular  beam  by 
cylindrical  beamlets  symmetrically  located  around  the  annular  ring  has  been 
discussed.  Here  again  the  possible  transverse  instability  inherent  in  this  ar¬ 
rangement  is  difficult  to  eliminate  and  also  required  phase  synchronism  is 
difficult  to  achieve. 

A  different  approach  was  employed  in  the  UTA  Wake-field  and  FEL 
program.  Multi-stage  acceleration  scheme  was  adapted  to  overcome  the  low 
transformer  ratio  and  still  ptovide  high  acceleration  gradient.  For  the  staging 
purpose,  we  choose  an  elliptical  geometry  because  of  its  inherent  staging 
capability.  Unlike  the  co-linear  geometries,  replenishing  the  energy-depleted 
driver  beams  is  practical  for  the  elliptical  structures  since  two  beam  aper¬ 
tures  of  unequal  size  around  the  foci  of  an  ellipse  provide  two  separate 
paths  for  the  driving  and  accelerated  bunches.  Although  the  transformer 
ratio  is  not  as  high,  adequate  acceleration  gradients  in  such  a  structure  are 
possible  when  the  driving  beam  passes  through  a  high-impedance  aperture 
but  is  replenished  periodically  as  shown  schematically  in  figure  1.2.  Since 
the  driving  beam  is  to  deliver  energy  in  a  single  stage,  it  is  possible  to 
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Figure  1.2.  Multi-stage  wake-field  accelerator  using  elliptical  cavities 
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do  so  at  a  high  rate.  Because  of  the  staging  capability,  it  is  now  possible 
to  obtain  the  correct  synchronism  to  overcome  the  usual  "phase  slippage" 
between  the  faster  accelerated  beam  and  the  slower  driver  beam.  Another 
unique  application  of  the  elliptical  wake-field  cavity  which  cannot  be  realized 
with  the  co-linear  geometry  is  that  it  can  be  used  as  electric  wiggler 
in  a  Free  Electron  Laser19  (FEL),  as  shown  schematically  in  figure  1.3.  Two 
separate  paths  along  the  focus  axis  of  the  elliptical  cavity  allow  the  lasing 
electron  beam  to  be  injected  in  the  opposite  direction  of  the  driving 
bunches.  The  lasing  electron  bunches  wiggles  in  the  rippled  wake  electric 
fields  produced  by  the  driving  bunches,  analogous  to  the  case  of  the  mag¬ 
netic  wigglers  in  conventional  FELs.  This  new  lasing  scheme  promises  suf¬ 
ficient  gains  in  the  soft  x-ray  region  so  that  it  might  be  feasible  to  operate 
it  as  a  soft  x-ray  laser  without  requiring  reflecting  mirrors. 

Based  on  these  ideas,  investigations  of  wake-field  acceleration  and  its 
applications  with  the  emphasis  on  the  elliptical  geometry  have  been  going 
on13,19'21,23  at  the  Center  for  Accelerator  Science  and  Technology  (CAST), 
UTA,  in  the  last  few  years,  part  of  which  is  the  subject  of  this  dissertation. 

The  purpose  of  this  dissertation  is  to  set  up  appropriate  theoretical 
means  to  analyze  the  wake  fields  excited  in  the  elliptical  cavity  which  is 
loaded  with  metallic  disks  or  a  medium.  The  wake  fields  and  wake  poten¬ 
tials  in  the  disk-loaded  cavity  can  be  evaluated  by  either  modal  analysis 
method  or  numerically  solving  the  Maxwell  equations  directly  in  the  time 
domain.  Because  analytical  calculations  are  only  Known  for  closed  cylindrical 
cavities,  the  wake  fields  in  an  elliptical  cavity  have  been  roughly  estimated 
from  those  in  the  equivalent  circular  cavity.  Furthermore,  no  method  has 
been  developed  to  calculate  wake  fields  in  a  dielectric-loaded  elliptical 


Figure  1.3.  Elliptical  cavity  as  an  electric  FEL  wiggler 
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cavity.  The  numerical  code  WELL,23  which  directly  solves  the  Maxwell  equa¬ 
tions  in  the  time  domain,  has  been  recently  developed  to  study  the  wake- 
field  effects  in  the  elliptical  pill-box  cavity.  This  particular  code  calc  dates 
the  wake  fields  and  wake  potentials  in  an  elliptical  pill-box  cavity  by  using 
the  finite  difference  method,  and  the  beam  aperture  effects  are  considered 
in  the  calculations. 

In  this  investigation,  new  methods  are  developed  to  calculate  the  wake 
fields  in  the  elliptical  structure  analytically  in  the  limit  of  vanishing  aperture. 
We  do  not  take  into  account  the  aperture  effects  based  on  the  recent  ex¬ 
perimental  work,14  which  showed  that  the  wake  fields  are  predominated  by 
a  few  lower-order  longitudinal  modes  even  when  the  ratio  of  the  aperture 
radius  to  the  cavity  radius  is  considerably  large.  At  lower  frequencies,  it 
was  known  that  the  Impedance  of  a  cavity  with  finite  aperture  and  that  of 
a  closed  cavity  are  practically  same.37  Therefore,  the  wake  fields  of  bunches 
which  is  long  compared  to  the  aperture,  or  the  long  range  wake  fields  of 
short  bunches  are  almost  the  same  for  the  two  different  structures.  Since 

wake-field  acceleration  schemes  require  very  short,  intense  bunches  to  excite 
strong  accelerating  field,  and  the  distance  between  the  driving  bunch  and 
accelerated  bunch  is  usually  large  compared  to  the  driving  bunch  length,  the 
wake  fields  in  the  limit  of  vanishing  aperture  is  a  good  approximation  for 
the  most  practical  cavities  as  long  as  apertures  are  not  too  large.  Even 

though  it  is  a  limiting  case  for  the  cavities  with  finite  apertures,  it  is  of 
particular  interest  since  the  analytical  solution  exist,  which  will  be  formulated 
in  this  investigation.  Therefore,  we  can  readily  estimate  maximum  energy 

gain  of  the  accelerated  particles  and  the  parametric  dependence  of  the  wake 

potentials  on  the  geometry  and  beam  parameters.  Furthermore,  it  becomes 
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possible  to  include  the  effects  of  dielectric  material  inside  the  cavity  with 
this  analytical  approach. 

A  short  review  of  the  concepts  of  the  wake  fields,  the  wake  potential, 
and  detailed  modal  analysis  method  are  presented  in  Chapter  2.  We  present 
analytical  formulations  of  the  wake  fields  and  wake  potentials  in  an  elliptical 
pill-box  cavity  in  Chapter  3.  The  wake-field  effects  in  the  dielectric-loaded 
elliptical  cavity  are  described  in  Chapter  4,  while  in  Chapter  5  discussions 
and  conclusions  are  drawn. 


;HAPTER  II 


NORMAL  MODE  EXPANSION  OF  THE  WAKE  POTENTIALS 
2.  1.  Introduction 

The  wake  fields  are  the  electromagnetic  fields  induced  by  a  bunch  of 
relativistic  charged  particles  passing  through  the  structure  of  varying  shape, 
and  their  concepts  are  important  not  only  for  the  applications  to  future  high 
current  compact  electron  accelerators,  but  also  for  the  calculation  of  energy 
loss  and  beam  stability  in  the  high  energy  particle  accelerators.  We  begin 
with  brief  introduction  of  the  concepts  of  the  wake  fields  and  wake  poten¬ 
tials. 

Consider  a  point  charge  Q  in  free  space.  When  a  charge  is  not 
moving,  the  electric  field  lines  are  directed  outward  in  all  radial  directions 
as  shown  in  figure  2.1(a).  But  the  electromagnetic  fields  carried  by  a 
relativistic  point  charge  Q  in  free  space24  is  Lorentz  contracted  into  a  thin 
disk  as  shown  in  figure  2.1(c)  with  angular  spread  of  the  order  of  1/y, 
where  y  being  the  relativistic  factor.  The  disk  actually  shrink  into  a  delta 
function  thickness  in  the  ultra-relativistic  limit  of  the  particle  velocity 
v  =  c,  as  shown  in  figure  2.1(d).  In  case  the  charge  moves  along  the  axis 
of  a  cylindrically  symmetric  pipe  with  perfectly  conducting  wall,  the  fields 
in  the  pipe  are  identical  to  the  free  space  fields,  because  the  sole  function 
of  the  pipe  wall  is  to  truncate  the  field  lines  by  terminating  them  onto 
the  image  charges  on  the  wall  as  shown  in  figure  2.2(a).  If  the  charge 
moves  off-axis,  or  if  the  cross  section  of  the  pipe  is  not  a  circle,  the  fields 
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trailing  charge 


Figure  2.2.  Fields  in  a  perfectly  conducting  smooth  pipe  (a)  and  in  a  structure 

with  perturbing  walls  (b) 
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in  the  pipe  are  perturbed  from  their  free  space  values.  But,  no  fields  are 
left  behind  the  moving  charge  since  the  fields  are  Lorentz  contracted  into 
a  delta  functioi  in  its  longitudinal  distribution. 

Therefore,  in  free  space  or  in  a  perfectly  conducting  smooth  pipe,  a 
particle  does  not  see  the  fields  carried  by  other  particle  unless  the  two 
particles  move  side  by  side  with  exactly  same  longitudinal  position,  in  which 
case  they  see  each  other’s  fields  but  do  not  experience  any  force  because 
the  electric  force  and  the  magnetic  force  cancel  exactly  in  the  limit  of  par¬ 
ticle  velocity  v  =  c. 

When  there  are  irregularities  or  discontinuities  in  the  structure,  the 
trailing  charge  still  will  not  see  the  direct  fields  in  the  wavefront  moving 
with  the  driving  charge.  However,  this  wavefront  are  scattered  from  the 
boundaiy  discontinuities,  and  this  scattered  radiation  can  reach  the  trailing 
charge  and  exert  forces  parallel  and  perpendicular  to  its  direction  of  motion. 
These  scattered  waves  are  called  the  wake  fields,  and  the  force  integrated 
over  the  total  passage  time  is  usually  called  the  wake  potential  (more 
precisely  delta  function  wake  potential). 

If  the  wall  is  smooth  but  not  perfectly  conducting,25”6  a  charge  will 
generate  behind  it  the  wake  fields  due  to  the  small  longitudinal  electric 
field  Ei  at  the  wall,  which  is  related  to  the  dissipation  in  the  wall.  How¬ 
ever,  the  effects  of  resistive  wall  are  much  smaller  than  other  factors.  Most 
of  the  wake  fields  come  from  effects  associated  with  discontinuities  and 
obstacles  in  the  structures.  Therefore,  in  the  following  section  and  thereafter, 
finite  conductivity  effects  will  not  considered  in  the  wake  fields  calculations. 
It  is  also  assumed  in  this  chapter  that  the  driving  charge  and  the  test 
charge  move  with  velocity  of  light  c,  and  their  paths  through  a  cavity  or 
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structure  are  not  significantly  changed  from  what  would  have  been  in  the 
absence  of  the  induced  wake  fields. 

2.  2.  Method  of  Calculation 

The  wake  fields  and  wake  potentials  in  a  cavity  can  be  evaluated  by 
either  modal  analysis  method  or  numerically  solving  the  Maxwell  equations 
in  the  time  domain.  In  the  modal  analysis  method,  the  eigenmodes  of  an 
empty  cavity  are  used  to  calculate  the  time  development  of  the  fields. 
Analytical  calculations  of  the  eigenmodes  are  only  known  for  cavities  with 
simple  geometries,  usually  for  closed  cylidrically  symmetric  cavities.  For  this 
reason,  early  studies  on  the  wake-field  cavities  were  mainly  concentrated  on 
the  circular  geometry.  For  non-axisymmetrical  structures,  one  has  to  use 
numerical  methods. 

The  straightforward  way  to  calculate  wake  fields  and  wake  potentials 
in  the  structure  is  to  solve  the  Maxwell  equations  numerically.  The  numeri¬ 
cal  code  TBCI27,28  has  been  extensively  used  for  the  cylidrically  symmetric 
structures,  and  MAFIA29  code  T3,  a  three  dimensional  version  of  the  code 
TBCI,  was  recently  developed.  Previously,  the  particular  code  WELL,22  the 
finite  difference  solution  of  the  Maxwell  equations  in  the  time  domain,  has 
been  developed  for  the  calculation  of  the  wake  fields  in  an  elliptical  cavity. 
In  this  code,  the  aperture  or  beam  hole  effects  are  included  in  the  cal¬ 
culation.  However,  we  need  not  take  into  account  the  aperture  effects 
seriously  for  the  most  practical  cavities  whose  apertures  are  usually  very 
small  compare  to  the  cross  sectional  area.  This  is  clearly  supported  by  the 
experiment,14  which  showed  that  the  wake  fields  are  predominated  by  a  few 
lower-order  longitudinal  modes  even  when  the  ratio  of  the  aperture  radius 
to  the  cavity  radius  is  considerably  large  (**1/3  for  the  cavity  1  in  that  ex- 
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periment).  Since  the  impedance  of  a  cavity  with  finite  aperture  and  that  of 
a  closed  cavity  are  practically  same  at  lower  frequencies,  the  wake  fields 
of  a  bunch  whose  length  is  long  compared  to  the  aperture  or  the  long 
range  wake  fields  of  a  shot  bunch  are  almost  the  same  for  the  two  dif¬ 
ferent  structures.  Since  we  are  more  interested  in  evaluating  the  long  range 
wake  fields  induced  by  a  short  driving  bunch,  a  closed  pill-box  cavity  is  a 
good  approximation.  Even  though  the  final  result  is  the  limiting  case  for 
the  cavity  with  beam  holes,  it  is  of  particular  interest  since  we  can  formu¬ 
late  an  analytical  solution  under  the  assumtion  of  zero  aperture.  The  lon¬ 
gitudinal  and  transverse  wake  potentials  can  be  calculated  to  any  degree  of 
precision  since  we  can  express  these  fields  analytically  in  closed  series  form. 
Also,  we  can  readily  estimate  the  maximum  energy  gain  of  the  trailing  par¬ 
ticles,  frequency  contents,  and  the  dependence  on  the  cavity  geometry  and 
beam  parameters. 

2.  3.  Normal  Mode  Analysis  of  the  Wake  Fields 

Consider  a  closed,  empty  cavity  with  perfectly  conducting  walls.  The 
electric  fields  and  the  magnetic  fields  induced  by  a  bunch  of  charged  par¬ 
ticles  traversing  the  cavity  can  be  derived  from  a  scalar  potential  £>(x,r)  and 
a  vector  potential  A(x,t)  together  with  proper  boundary  conditions.  First,  let 
us  start  from  the  Maxwell  equations  (in  MKS  units) 


V  x  E  = 

V  ft  ’ 

(2.1) 

VXB=„,J  +  ^, 

(2.2) 

v-e=£. 

(2.3) 

V  •  B  =  0, 

(2.4) 
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where  £o  and  /r<j  are  respectively  the  permittivity  and  permeability  of 
vacuum.  By  introducing  a  vector  potential  A  (x,  r)  and  a  scalar  potential 
<p(x,t)  which  are  given  by  the  relations 


B  -  V  x  A, 


the  inhomogeneous  Maxwell  eqs.  (2.2)  and  (2.3)  become 


(2.5) 

(2.6) 


’  c2  at 


C2  dif 


t*0  J> 


vv= 


(2.7) 

(2.8) 


where  the  Coulomb  gage,  V-A  =  0,  is  used  in  deriving  eqs.  (2.7)  and  (2.8). 
The  vector  potential  A  (x,r)  and  a  scalar  potential  <p(x,t)  due  to  current 
and  charge  sources  can  be  expanded  in  terms  of  the  othogonal 
eigenfunctions  su(x)  and  <px(x)  as 


A(x,/)  =  £  qx(t)  ai(x), 
<P(*,  0  =  £  n(t)  <px{x), 


(2.9) 

(2.10) 


where  a*  and  <px  are  the  sets  of  vector  and  scalar  eigenfunctions  of  an 
empty  cavity  respectively,  which  satisfy  the  following  homogeneous  Helmholtz 
equations  and  boundary  conditions: 


V2ai  +  (— )2a;  =  0, 
c 

(2-11) 

V-aa  =  0  in  V, 

(2.12) 

n  x  a*  =  0  on  the  boundary, 

(2.13) 
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and 


=  o. 

(2.14) 

<px  =  0  on  the  boundary, 

(2.15) 

where  n  is  the  unit  vector  normal  co  the  wall.  In  general,  (ox  *  Q*. 

Substituting  eqs.  (2.9)  and  (2.10)  into  (2.7),  dotting  this  result  with  ar,  and 

integrating  over  the  cavity  volume  V  gives 

(qx  +  (ox2  qx )  fax  •  ax-  dV  =  poc2  /J  •  a;  dV.  (2.16) 

If  the  eigenfunctions  ax  are  normalized  by 

y  fax- ax- dV=  ux  dxx-,  (2.17) 

where  ux  is  the  normalizing  factor  and  dxx-  is  the  Kronecker  delta,  then 
eq.  (2.16)  becomes  simply 

q\  +  m2qx  =  ~fJ*xdV.  (2.18) 

Similarly,  beginning  with  the  Poisson’s  equation  given  in  (2.8),  we  can  get 

the  expansion  coefficients  for  the  scalar  potential,  which  is  given  by 

rx(t)  =  YfJp<pxdV,  (2.19) 

where  <px  are  normalized  by 

y  fV<px  *  dV  =  Tx  dxx-.  (2.20) 

The  expansion  coefficients  qx(t)  and  rx(t)  can  be  determined  from  eqs. 
(2.18)  and  (2.19)  provided  the  source  terms  J,  p  and  the  eigenfunctions 
ax  and  <px ,  which  will  be  derived  in  the  following  chapter,  are  known.  Ac- 
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cordingly,  we  can  rewrite  the  electric  field  E  and  the  magnetic  induction 
B  in  terms  of  these  quantities  as 


E(*  0  =  -  5  (^(0  a*(x) +  V0i(x)) , 

(2.21) 

B(x,  l)-J  qx{t)  V  x  ai(x). 

(2.22) 

These  are  the  wake  fields  left  by  the  sources.  In  order  to  solve  eq.  (2.18) 
and  (2.19),  we  have  to  know  the  source  terms. 

2.  4.  Wake  Potentials  in  a  Cavity 

2.  4.  1.  Longitudinal  Wake  Potential 

Consider  the  cavity  shown  in  figure  2.3,  where  a  driving  charge  Q  is 
assumed  to  enter  the  cavity  at  z  =  0,  t  =  0  and  to  exit  at  z  -  L. 
The  driving  charge  Q  is  passing  through  the  cavity  with  velocity  v  -  c,  fol¬ 
lowed  by  a  test  charge  by  a  fixed  distance  zo  in  the  z-direction.  Let  r  and 
ro  be  respectively  the  transverse  coordinate  of  the  test  charge  and  that  of 
the  driving  charge.  The  delta  function  longitudinal  wake  potential  Wz  is 
defined  as  the  energy  gained  by  the  unit  negative  point  charge  following 
at  a  distance  zo  behind  the  driving  charge,  divided  by  Q,  or 

W^(r,  zo)  =  -q  ft dz  Ez( r,  z,  Ot-c+ioyc.  (2.23) 

The  source  terms  due  to  the  exciting  charge  are 

P( x.  0  =  2  <*(r  ~  *>)  <Kz  ~  v/),  (2.24) 

J(x,  t)  =  zpv  =  zpc(x,  t),  (2.25) 

where  z  is  the  unit  vector  in  the  z-direction  and  particle  velocity  of 
v  =  c  is  assumed.  After  substituting  eq.  (2.25)  into  (2.18),  equation  for  qx 


Figure  2.3.  The  driving  and  accelerated  beams  in  a  cavity 
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becomes 


qx  +  o)x2  qx  = 


Qc_ 

2ui 


0,  t<  0 

axt  (it),  ct ),  0  <  t  <  L/c 
0,  t  >  L/c 


(2.26) 


where  ait  is  the  2-component  of  t>e  vector  eigenfunction  at. 

With  the  initial  conditions  <?i(0)  =  qx(0)  =  0,  i.e.,  no  field  exists  before 
the  driving  charge  Q  enters  the  cavity,  the  solution  for  qx  is 


qx (r)  =  ~  /J’"’0,  Wc)  dt'  sin£u^(r  -  /')  ct)  for  t  >  0. 


(2.27) 


Similarly,  substituting  source  term,  eq.  (2.24),  into  (2.19)  gives 


rx(t)  = 


_Q_ 

2Tx 


0,  t  <  0 

•  <px  (ib,  ct),  0  <  t  <  i/c 

0,  /  >  t/c 


(2.28) 


Substituting  these  expressions  of  the  expansion  coefficients  into  eqs.  (2.21) 
and  (2.22),  we  can  get  the  electric  field  E  and  the  magnetic  induction  B 
at  any  position  in  a  cavity,  which  are  then  used  to  calculate  the  longitudinal 
and  transverse  wake  potentials  by  way  of  eqs.  (2.23)  and  (2.32).  By  using 
the  electric  field  Ez  and  imposing  causality,  i.e.,  no  wake  field  will  be 
produced  ahead  of  the  driving  charge,  Bane  et  al.12  obtained  the  following 
equation  for  the  delta  function  longitudinal  wake  potential: 


wz{r,  n»,  zo)  =  ^  £  ;0L  dz  (qxtr^r)  <«,(«•,  A  + 


kV(tb)  Vx(r)  ,0), iN 

■L  “*<7* 


(2.29) 


where  V\  is  a  complex  conjugate  of  Vx.  Notice  that  the  longitudinal  wake 
potential  is  solely  expressed  in  terms  of  the  z-component  of  the  vector 


eigenfunction,  and  is,  in  general,  functions  of  the  transverse  positions  of  the 
driving  charge  and  a  test  charge.  In  eq.  (2.29),  Vx  is  defined  as 

W  =  So  dz  exp(^~ )  <wx(r,  2).  (2.30) 

Once  the  delta  function  wake  potential  given  in  eq.  (2.29)  has  been 
calculated,  it  can  be  used  as  a  Green’s  function  to  calculate  th  j  bunch 
wake  potential  Uz  in  and  behind  an  arbitrary  charge  distribution  A(z).  The 
longitudinal  potential  for  the  distributed  charge  becomes 

C/*(zo)  =  C  A (z  -  z0)  Wz(z)  dz.  (2.31) 

2.  4.  2.  Transverse  Wake  Potential 

Consider  again  the  driving  charge  Q  and  a  test  charge  in  a  cavity  as 
shown  in  figure  2.3.  The  delta  function  transverse  wake  potential  Wx(z0)  is 
defined  as  the  transverse  momentum  kick  experienced  by  a  test  charge  fol¬ 
lowing  at  a  distance  zo  behind  a  driving  charge,  divided  by  Q,  or 

Wi.(zo)  (Ex  +  c(z  X  B)j.)  t_(l+ro)/c 

-<««.*  (2.32) 

where  the  transverse  gradient  operator,  Vx,  is  defined  as 

V^V-2^.  (2.33) 

Analogously  to  the  longitudinal  case,  the  transverse  wake  potential  Wx(z0) 
can  be  calculated  by  substituting  the  vector  and  scalar  potentials  which  are 
expanded  in  terms  of  the  eigenfunctions  into  eq.  (2.32).  But,  we  may  use 
a  rather  simple  formula  known  as  the  Panofsky-Wenzel  theorem30 


3^Wx(zo)  =  V±Wz(z0) 


(2.34) 


to  derive  the  transverse  wake  potential  fron  the  knowledge  of  the 
longitudinal  wake  potential.  Eq.  (2.34)  is  valid  for  the  selected  cavity 
geometries,37  e.g.,  for  cylindrically  symmetric  strut  tares  for  any  value  of  zo, 
for  closed  cavities  of  arbitrary  shape  for  zo  >  L,  and  for  structures  with 
end  plates  normal  to  the  path  of  particles  for  any  value  of  zo.  Thus,  we 
are  going  to  use  the  Panofsky-Wenzel  theorem  in  the  next  chapter  to  derive 
the  transverse  wake  potential  in  an  elliptical  pill-box  cavity.  Combining  eqs. 
(2.29)  and  (2.34),  the  delta  function  transverse  wake  potential  becomes 

(2.35) 

and  the  transverse  wake  potential  for  an  arbitrary  charge  distribution  A(z) 
can  be  obtained  from 


Ux(zo)  =  A(z  -  zo)  Wx(z)dz. 


(2.36) 


CHAPTER  m 


WAKE  POTENTIALS  IN  AN  ELLIPTICAL  CAVITY 

3.  1.  Eigenmodes  in  an  Elliptical  Pill-Box  Cavity 

In  chapter  2,  we  have  derived  the  longitudinal  and  transverse  wake 
potentials  solely  in  terms  of  the  z-component  of  the  vector  eigenfunctions, 
which  is  the  solution  of  the  homogeneous  Helmholtz  equation  in  a  source 
free  cavity  with  appropriate  boundary  conditions.  The  eigenfunctions  in  an 
elliptical  pill-box  cavity  are  derived  in  this  section,  and  the  wake  potentials 
are  formulated  in  the  following  sections. 

Consider  an  elliptical  pill-box  cavity  with  the  major  axis  2rb  and  the 
minor  axis  2y*  as  shown  in  figure  3.1.  For  a  cavity  of  elliptical  cross  sec¬ 
tion,  the  eigenfunctions  can  be  found  in  terms  of  known  functions  by  trans¬ 
forming  the  Cartesian  coordinates  to  the  confocal  elliptical  coordinates  as 
shown  in  figure  3.1.  In  these  coordinates,  the  boundary  conditions  on  the 
elliptical  cavity  walls  are  easily  satisfied. 

The  confocal  elliptical  coordinates  (£,  rj,z)  are  related  to  the  Cartesian 
coordinates  (x,y,z)  by  the  following  equations: 

x  =  h  cosh£  cos?;, 

y  =  h  sinh?  siiu/,  (3.1) 

where  h  is  the  semi-interfocal  distance  of  the  confocal  ellipses.  The 
Helmholtz  equation  is  then  transformed  to 


25 


/  |  =  const. 
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vV  +  = 


&+£>  +  £  +  *V-a 


/i2(sinh2£  +  sin2^) v  d£2  dr]2  dz 2 


(3.2) 


where  k  is  the  free  space  wave  number. 


Following  the  method  of  separation  of  variables,  we  seek  to  Find  solu¬ 
tions  of  the  form 


V=f(Z)80 7Mz). 


(33) 


Substituting  eq.  (3.3)  into  (3.2)  and  dividing  by  ip,  the  Helmholtz  equation 
is  split  into  three  ordinary  differential  equations: 


dt 


(a-2q  cosh2£)/(|)  -  0, 


+  (a  -  2q  cos2j/)£(i/)  =  0, 

^f-  +  k,2w(z)  =  0. 
dz 


(3.4) 

(3.5) 

(3.6) 


where  a  is  an  arbitrary  separation  constant.  The  parameter  q  is  defined  as 


2l2 


2q  = 


kch 


2  ’ 


(3.7) 


and  kt  and  kc  are  arbitrary  constants,  but  are  related  to  the  free  space 
wave  number  k  by 

k>  =  k2  -  kt2.  (3.8) 

The  above  equations  (3.4)  and  (3.5)  are  called  the  Mathieu  equations.31-52 
Solutions  to  these  equations  are  the  Mathieu  functions,  details  of  which  are 
given  in  Appendix  A. 

With  the  z-dependence  of  exp (ifc*z)  in  eq.  (3.6),  we  are  interested  in 
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the  solution  of  the  two  dimensional  Helmholtz  equation 

VxV  +  fcV  =  0  (3.9) 

in  an  elliptical  pill-box  cavity.  Here,  the  transverse  Laplacian  operator  V.2 
is  defined  by 


The  particular  modes  that  we  are  considering  are  the  modes  which 
have  longitudinal  acceleration  fields.  Therefore,  only  Transverse  Magnetic 
(TM)  modes  will  be  considered  for  the  calculation  of  the  wake  potentials. 
For  modes  TM  to  z,  we  may  express  the  fields  in  terms  of  the  vector 
potential  A,  or  equivalently  in  terms  of  the  vector  eigenfunction  a*,  having 
only  z-component.  Hence  tp  =  ait  in  eq.  (3.9).  Any  combinations  of  the 
Mathieu  functions  in  the  form  of  tp  =  /(£)  g(rj)  are  solutions  of  eq.  (3.9).  In 
addition  to  the  boundary  conditions  which  must  be  satisfied  on  the  walls 
of  the  cavity,  the  following  two  conditions  must  be  satisfied,  i.e., 

(i)  continuity  of  function  tp  on  the  interfocal  line: 

mv)=m-v)-  (3.ii) 

(ii)  continuity  of  gradient  of  function  tp  on  the  interfocal  line: 

^  (V>(£,  »7)){-o  =  -  ^  (V<£>  .  (3.12) 

Among  the  possible  combinations,  the  only  permissible  form  of  the  solution 
which  satisfies  above  two  conditions  are 

m  oo 

tp  =  2  SmSem^,q)sem(r},q)  +  X  CmCem($,  q)  cem(jh  q), 

ra  »  1  m  “  0 


(3.13) 
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with  the  z-dependent  factor  exp(iTcaz)  being  omitted.  Here,  Cm  and  Sm  are 
abitrary  constants.  The  function  cem(q,  q)  and  sem(q,q)  are,  respectively,  the 
even  and  odd  type  Mathieu  functions  of  the  first  kinds  of  integral  order, 
and  Cem($,  q)  and  Sem(£,  q)  are  the  modified  Mathieu  functions  of  the  first 
kinds  of  integral  order. 

We  first  notice  that  symmetry  of  the  function  ip  is  determined  by 
cem(jl,q)  and  sem(t],q).  The  terms  in  the  first  summation  are  odd  and  the 
terms  in  the  second  summation  are  even  with  respect  to  q.  From  Appendix 
A,  we  have 

sem(r],  q)=  -  sem(-rj,  q),  (3.14) 

cem(rj,  q)  =  cem(-rj,  q).  (3.15) 

Therefore,  the  first  summation  is  always  zero  on  the  median  plane  (y  = 
0  plane)  where  the  wake  potentials  are  to  be  evaluated,  and  as  a  result 
it  does  not  contribute  to  the  calculation  of  the  longitudinal  and  transverse 
wake  potentials.  For  this  reason,  we  are  now  considering  only  the  TM 
modes  of  even  type  solutions,  rpm  =  Cem(^,q)cem(rj,q),  for  the  wake  potential 
calculation. 

The  boundary  conditions  are 

ip  =  0  (3.16) 

at  the  elliptical  boundary  where  |  =  £a-  There  remains  the  condition  that 
no  tangential  component  of  a*  exists  at  the  end-plate  walls  at  z  =  0  and 
d,  which  is  satisfied  if  we  choose 

kt  —  ,  p  =  0, 1, 2, ... , 


(3.17) 
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where  d  is  the  gap  distance  of  the  pill-box  cavity  as  shown  in  figure  3.2. 
From  eq.  (3.16),  we  have 

Cem(£a,  q)  =  0,  (3.18) 

which  are  satisfied  by  particular  values  of  q.  Let  qmn  be  the  n,h  root  of 
eq.  (3.18)  for  a  given  mode  m.  Then,  we  can  calculate  resonant  frequency 
from  the  root  qmn  .  From  eqs.  (3.7)  and  (3.8),  the  resonant  frequency  is 

UJmnp 
C 

From  this  equation,  we  see  that  the  wavelength  of  the  dominant  mode 
becomes 


4tfn 


<>2 


(3.19) 


Aoio  =  (3.20) 

where  ec  is  the  eccentricity  of  the  boundary  ellipse,  ec  =  l/cosh|t.  Figure 
3.3  is  a  plot  of  qoi  as  a  function  of  the  eccentricity  ec.  The  ratio  Aou/jCb 
is  plotted  against  ee  in  figure  3.4.  In  this  figure  we  see  that  as  eQ  -*  0,  i.e., 
as  the  ellipse  tends  to  a  circle,  the  ratio  approaches  2.61,  which  is  the 
ratio  of  the  wavelength  to  the  radius,  A/ra  =  2jt/xoi,  for  the  TMoio  mode  of 
a  circular  pill-box  cavity,  where  *oi  is  the  first  zero  of  Bessel  function  Jo(x). 

The  transverse  components  of  the  vector  eigenfunction  can  be  easily 
found  from  the  z-component.  These  eigenmodes  are 


Ce m(£»  ^ma)  C£m(7*  ?mn)  CO 


(3.21) 


(£, ?mn) cem(jjf  qmn)  sin^7 ^  , 


(3.22) 


wavelength  of  the  fundamental  mode  as  a  function  of  eccentricity  ec 
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(3.23) 


where 


Dm  a  — 


_ 


cosh2|  -  cos2?;  ** 

2  J  ' 


(3.24) 


The  resonant  mode  patterns  for  some  of  the  lower-order  modes  are  shown 
in  figures  3.5  -  3.8.  Now,  we  proceed  to  calculate  the  wake  potentials  using 
eigenmodes  given  in  eqs.  (3  21)  -  (3.23). 


3.  2.  Wake  Potentials  in  an  Elliptical  Pill-Box  Cavity 
3.  2.  1.  Longitudinal  Wake  Potential 

The  longitudinal  and  transverse  wake  potentials  on  the  focus  axis  of 
an  elliptical  cavity  are  of  particular  interest  since  when  the  driving  charge 
Q  is  passing  through  one  of  the  foci  of  the  cavity  where  £  =  rj  =  0  (driving 
beam  path),  all  the  electromagnetic  fields  radiated  will  get  focused  again  at 
the  other  foci  where  $  =  0  and  rj  =  n  (accelerated  beam  path)  in  first  order 
approximation,  which  are  then  used  to  accelerate  a  trailing  charge. 

Consider  a  test  charge  which  is  traveling  through  the  cavity  where 
|  =  0  and  t]  =  n  (accelerated  beam  path)  and  trailing  the  driving  charge  Q 
by  a  fixed  distance  zo  in  the  z-direction.  The  delta  function  longitudinal 
wake  potential  Wu  or  the  energy  gained  by  the  unit  negative  test  charge 
is  given  by  eq.  (2.29): 

FPz(zo)  =  ^  cos(^)  forzo  >  0,  (2.29) 

where  Vx  is  the  voltage  induced  by  a  point  charge,  and  the  normalizing 


H 


H 


Figure  3.5.  Configuration  of  the  resonant  mode  a ^  in  an  elliptical  cavity  (m  =  0, 


«lz 


Figure  3.6.  Configuration  of  the  resonant  mode  a in  an  elliptical  cavity  (m  =  l,  n  =  l,  p  =  0) 
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EH 


Figure  3.7.  Configuration  of  the  resonant  mode  a ^  in  an  elliptical  cavity  (m  = 


•K  *B 
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1 


Figure  3.8.  Configuration  of  the  resonant  mode  a ^  in  an  elliptical  cavity  (m=3,  n=l,  p=0) 


factor  ux,  which  can  be  interpreted  as  the  stored  energy  in  the  mode,  is 
given  by  eq.  (2.17): 

^  fax-ax' dV=  uxdw .  (2.17) 

After  substituting  the  eigenfunction  in  eq.  (3.21)  into  eq,  (2.30),  the  voltage 
Vx  on  the  driving  beam  path  and  any  position  in  a  cavity  become 

nw  =  f0  dz  exp(^)  a»,(0, 0,  z) 

=  Cem(0,  qmn)  cem( 0,  qmtt)  dz  exp(ff  )  cos^z,  (3.25) 

Vx(r)  -  J^dz  exp(f^)  «,(£,  rj,  z) 

-  Cem($,  qmn)  cem(rj,  qmn)  dz  exp(~)  cos ^z,  (3.26) 

and  further 


2Cmn(—)2  ,,  , 

Vx  Vx  = - - - (1  -  (-l)p  cos^*% 

((f)2 -(^)2)2 


(3.27) 


where  the  coefficient  Cmn  is  given  by 


Cmn  ~  Ce m(£,  ^mn)  CCm(0,  ^mn)  qmn)  C6 m(0,  qm n). 


(3.28) 


The  stored  energy  ux  in  the  mode  becomes 


Wrnnp  —  ~  fax'  Hi'  dV 


(CO  mnp^2 

% '  C~2  -  /  Cem2(£,  ?mn)  cem\ij,  qmn)  cos2f~)z  dV. 


(3.29) 
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h2 

After  integration  of  eq.  (3.29)  using  dV  -  —  (cosh2£  -  coslrj)  d%  drj  dz,  it  is 
further  reduced  to 


Wmnp  ~ 


/Wmn?s2 


8 


4*  Cea\l  qmtl)  x  (c .sh2|  -  ©m)  d£, 


(3.30) 


where 

©m  =  ^  jf  cem\ij,  qmn)  cos2>7  drj 
A^A[m)  +2A^A&2  , 

_  r-0 

hA[m))2  +  2A$lxA<i% , 

z  r-0 

and  A$m)  in  eq.  (3.31)  are  the  coefficients  for  the  series  representation  of 
the  Mathieu  function  cem(rj,q).  Substituting  eqs.  (3.27)  -  (3.31)  into  (2.29) 
gives  the  delta  function  longitudinal  wake  potential  Wz. 


form  =  0, 2, 4, ... , 


form  =  1, 3, 5, ... , 


(3.31) 


On  the  accelerated  beam  path  (£  =  0,  rj  =  n ),  the  delta  function  lon¬ 
gitudinal  wake  potential  becomes 

,  -  -  .  erC..  (l-t-l/cos2^) 

WJzo)  =  -^2  2  2  - - - r, - C—L  cosC51^),  (3.32) 

,0  n- 1  p-0  V  C  *' 


where  eP  =  Vi  for  p  =  0  and  1  for  p  *  0,  and  Nmn  is  given  by 


^mn  ~  4  (It  ?mn)  X  (C0sh2£  ©m )  d£. 


(3.33) 


The  coefficient  Cmn  can  be  expressed  in  terms  of  A (rm)  as 

Cmn  =  ( C£m(0 ,  (Jmn))  C£m(0,  q mn)  C€ m(^t»  <?mn) 
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\2aW 

lf»0 
/ 


form  =  0, 2, 4, ... , 

4 


(3.34) 


1A&, 

V-° 


form  =  1, 3, 5, ... , 


When  the  boundary  ellipse  tends  to  a  circle,  the  confocal  hyperbolae 
in  figure  3.1  become  radii  of  the  circle  r,  and  the  confocal  ellipses  become 
concentric  circles  of  that  radius.  In  this  case,  we  get  the  longitudinal  wake 
potential  Wz  for  a  circular  pill-box  cavity,  which  is  expressed  analytically11,12 
in  the  form 


fVz(zo)  = 


£o nd  n=tl  p»0 


o  ^onp  d\ 

P  cos — ) 


X<*Jl  fron) 


cos(- 


^onp  ZO 


% 


(3.35) 


where  x<*  is  the  n,h  zero  of  Bessel  function  Jo(x).  We  can  easily  show  that 
eq.  (3.32)  is  exactly  reduced  to  eq.  (3.35)  by  using  the  limiting  properties 
of  the  Mathieu  functions.  The  functions  cem{rj,q)  and  Cem($,  q)  degenerate 
into  the  following  forms  as  the  semi-interfocal  distance  h  -*  0  and  £ 
while  keep  the  product  /icosh£  -*■  r.  From  eq.  (3.7),  we  see  that  q  -*•  0  as 
h  -*  0.  Thus,  as  q  -*  0,  we  have  from  Appendbc  A 


i  W  , form  =  0 

cem{rjy  q)  ■*  CQSm<p  f 6r  m  *  0 

> 

(3.36) 

Cem(M,  q)  Pm  Jm(kcT) 

(3.37) 

A -*■  0  (except  A^  1  for  m  *  0  andA^  -*•  Vv£) 

(3.38) 

where  pm  is  a  constant  multiplier  and  /m(r)  is  the  Bessel  function  of  the 
first  kind  of  order  m.  Using  above  equations,  Cmn  in  the  numorator  of  eq. 
(3.32)  becomes  simply 
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Cmn  —  (CCm(0,  ^mn))  C6m(0 ,  (Jmn)  C6m(?t,  Qma) 

■*  (pm  /«( 0)  cosnvp'f  =  0,  for  m  *  0 

-  (po  Jo(0)f/2  =  po2/2 ,  for  m  =  0.  (3.39) 

It  is  apparent  from  eq.  (3.39)  that  contributions  from  m  *  <)  modes  become 
zero  as  expected  as  an  ellipse  tends  to  a  circle.  The  denominator  in  eq. 
(3.32)  for  m  =  0  becomes 

9oo  4'  CeoU  ?«.)  cosh2|  f- jfykcr)  rdr 

(3.40) 

where  kc  =  ^on/ra  for  a  circular  pill-box  cavity  of  radius  ra.  By  substituting 
eqs.  (339)  and  (3.40)  into  (3.32),  we  can  get  exactly  the  same  expression 
for  the  wake  potential  in  a  circular  pill-box  cavity  as  in  references  11  and 
12.  No  closed  expression  is  known  for  the  infinite  sum  in  eqs.  (3.32)  and 
(3.35),  which  must  be  evaluated  numerically. 

The  wake  potential  for  the  distributed  charges  can  be  calculated  from 
eq.  (2.31).  For  a  Gaussian  bunch  with  rms  bunch  length  of  a, 

A(z)=dbexp($)'  (3-41) 


the  bunch  wake  potential  is 


Ufa)  - 


V2/;r 


-  -  ePCmB(l-(-l)pcos~1^) 

r  \  C 


2  22 


ttfido  m  _0  a-i  p-o 


<7mn  M 


mo 
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w  rl o  ✓  Z  x  ^mnp(^0  Z) 

x  e*P(^j)  cos -  c - 


(3.42) 


When  2o  >  >  cr,  the  bunch  wake  potential  becomes 

£pC„.(l-(-l)pcos^> 
x  - — - cos - K=-. 

?mn”mn  C 


(3.43) 


It  is  apparent  from  this  equation  that  contributions  from  the  modes 
whose  resonant  wavelengths  are  much  shorter  than  the  bunch  length  a  be¬ 
come  negligible.  For  the  dominant  mode,  eq.  (3.43)  can  be  written  as 


_  _  2nd  2g, 

*"  ea^declSNm  1  (1  tKX)“T’ 


(3.44) 


where  S  =  nxbyb  is  the  cross  sectional  area  of  the  cavity,  X  is  the 
wavelength  of  the  dominant  mode,  ec  is  the  eccentricity  of  the  boundary 
ellipse.  The  coefficients  Coi  and  Noi  are  found  from  eqs.  (3.33)  and  (3.34). 
From  this  equation,  we  see  that  the  wake  potential  scales  as  m~2  and 
S'1.  When  the  boundary  ellipse  is  very  close  to  a  circle,  i.e.,  ec  <  <  1, 
the  wake  potential  scales  as  ec~2. 

3.  2.  2.  Transverse  Wake  Potential 


As  discussed  in  Chapter  2,  we  can  derive  the  transverse  wake  potential 
in  a  pill-box  cavity  from  the  longitudinal  wake  potential  by  using  the 
Panofsky-Wenzel  theorem,  eq.  (2.35).  That  is 

,  x  vx  KiVj.  Vi  .  ,(t)izo x 

Wi(a)=?c'2^rs,n(T- >• 


(2.35) 
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On  the  accelerating  beam  path  where  £  =  0  and  rj  =  7t,  we  have 


V\  VxVx  =  Cem( 0,  O')  cem(Q,  qx)  cem(n,  qx)  VxCem(0,  qx) 


i^dztx p(~)  cos^zj 


(3,>5) 


Since  the  driving  and  test  charges  are  assumed  to  move  on  each  focus  ;xis 
which  are  parallel  to  the  z-axis  and  on  the  plane  of  symmetry  (y  =  0 
plane),  only  5-component  of  the  transverse  wake  potential  exists  at  the  foci 
of  the  elliptical  cavity.  Hence, 


•  VxCem(0,  qx)  =  I  lim  ^  Cem (£,  qx) 


f  5  {2r?AW ,  for  m  =  0,  2, 4, ... 

-  "r-0 

(2r+l)2^4^i ,  form  =  1,3, 5,  ... 

«r-o 


(3.46) 


Therefore,  the  delta  function  transverse  wake  potential  is  written  as 


cpCr(l -(-!)> cos^f) 

W±(zo)  =  -~2  - r; - —  sin(^o), 

eared  x  (OxqxNx  vc 


(3.47) 


where 


Ci  —  C£m( 0,  <7mn)  CCm(0>  ^mn)  CCm(^>  (^mn)  ^ xCCm (0,  ^mn) 


~  )3  (I  (^g0 ) ,  for  m  =  0, 2, 4, ... 

”  r-0  r-0 

1 -r1  ( 2 )3  ( 2  (>+  l)2/l?2i ) ,  for  m  =  1, 3, 5, ... 

n  r-0  r-0 


(3.48) 


For  the  driving  beam  with  a  Gaussian  charge  distribution,  the  transverse 
bunch  wake  potential  becomes 
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Ux(z0)  =  JZm  dz  X(z)  JV±(z0  -  z) 


/a  i  »  «  oo  £p  Cmn  (1  (  l)1*  COS  **  ) 

cV2/7r -2  E  E  —  c 


e&tdo 


m*0n*l  p»0 
,2 


ft>mnp  ^mn  Nmii 


(3.49) 


When  20  >  >  o',  the  transverse  bunch  wake  potential  becomes 


i  i  e*p(-c^)V/2) 

£o7tura=0n=l  p*0 


c 

to^mnpd 


ePCmn(l-(-l)pcos^) 
- — - sin - c-. 

Wmnp  ^mn  jVmn  C 


(3.50) 


3.  3.  Numerical  Examples 

The  longitudinal  and  transverse  wake  potentials  in  an  elliptical  pill-box 
cavity  are  calculated  by  using  eqs.  (3.42)  and  (3.49).  The  results  are  plotted 
in  figures  3.10  -  3.14.  In  these  calculations,  a  driving  beam  with  a  Gaussian 
charge  distribution  of  a  =  5  mm  is  assumed.  We  choose  the  same  cavity 
dimensions  bunch  length  as  in  ref.  22  to  compare  the  results  of  two 
methods.  The  cavity  dimensions  are  shown  in  figure  3.9. 

Figure  3.10  is  the  plot  of  the  longitudinal  wake  potential  on  the  ac¬ 
celerated  beam  path,  in  which  only  5  modes  are  included  in  the  calculation. 
The  modes  summation  converges  as  the  number  of  modes  are  increased. 
Figures  3.11  and  3.12  are  the  curves  of  the  longitudinal  wake  potentials 
respectively  on  the  accelerated  beam  path  and  on  the  driving  beam  path, 
in  which  different  number  of  modes  are  included  in  the  mode  summation 
(solid  lines  for  24  modes  and  broken  lines  for  12  modes).  From  figures 
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3.10  -  3.12,  we  readily  find  very  rapid  convergence  of  the  mode  summation, 
indicating  clearly  that  the  wake  potentials  are  predominated  by  a  few  lower- 
order  modes.  Even  5  modes  summation  (figure  3.10)  gives  a  well-defined 
curve. 

On  the  accelerated  beam  path,  we  obtain  about  125  MeV/m IfxC  ac¬ 

celeration  gradient,  while  about  110  MeV/m/^C  was  obtained  by  Y.  Chin22 
(figures  3.17  and  3.18). 

The  path  length  for  the  fields  radiated  from  one  focus  axis  to  reach 
the  other  after  one  reflection  at  the  elliptical  boundary  is  10  cm  for  this 
cavity.  The  peak  accelerating  potential  on  the  accelerated  beam  path  (point 
A  in  figure  3.11)  is  the  very  fields  radiated  from  the  driving  beam  path 
(point  B  in  figure  3.12)  and  focused  again  after  one  reflection  at  the  el¬ 
liptical  wall.  From  figures  3.11  and  3.12,  we  can  see  that  the  delay  distance 
between  the  point  A  and  point  B  is  about  10  cm. 

The  transverse  wake  potential  on  the  accelerated  beam  path  and  that 
on  the  driving  beam  path  are  shown  in  figures  3.13  and  3.14  respectively 

by  the  broken  lines.  The  corresponding  longitudinal  wake  potentials  are 
plotted  together  in  order  to  see  if  we  can  find  the  positions,  such  as  the 

point  C  in  figure  3.13,  at  which  the  accelerating  potential  is  large  while 
the  transverse  potential  is  small.  The  charge  to  be  accelerated  should  be 

positioned  at  such  a  point  in  order  to  avoid  large  transverse  deflection 
during  acceleration.  However,  we  also  found  that  at  the  point  of  peak  ac¬ 
celerating  potential,  the  transverse  wake  potential  is  not  zero.  The  ac¬ 
celerated  particles  positioned  at  this  peak  will  be  deflected  and  efficiency 
of  acceleration  will  drop  gradually  because  of  the  displacement  from  the 
focus  axis  and  proper  phase. 


47 


The  longitudinal  wake  potentials  for  different  size  cavities  are  plotted 
in  figures  3.15  -  3.17.  Figures  3.18  and  3.19  are  the  curves  for  the  lon¬ 
gitudinal  and  transverse  wake  potentials  calculated  by  pure  numerical 
method,  the  code  WELL,22  which  include  the  effects  of  the  beam  apertures 
of  !cm  diameter.  Comparing  figures  3.13  and  3.14  with  3.18  and  3.19,  very 
good  agreements  in  magnitudes  and  frequencies  are  observed. 


5  inodes  summation 
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Figure  3.10.  Longitudinal  wake  potential  on  the  accelerated  beam  path 
(Q  =  1  pC,  a  -  5  mm,  5  modes  are  included) 
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Figure  3.16. 


Figure  3.1 


CHAPTER  I\ 


WAKE  POTENTIAL  IN  A  DIELECTRIC-LOADED 
ELLIPTICAL  WAVEGUIDE 


4.  1.  Introduction 

As  discussed  in  Chapter  3,  disk-loaded  elliptical  cavity  is  capable  of 
producing  an  acceleration  gradient  greater  than  100  MeV/m/^C.  But  one  dis¬ 
advantage  of  using  this  cavity  as  a  wake-field  acceleration  device  is  its  un- 
cancellable  transverse  wake  fields  even  when  alternately  rotated  stages  are 
used  to  minimize  net  transverse  deflection.  The  transverse  wake  fields  are 
useful  as  an  electric  wiggler  fields  for  the  FEL  applications,19  but  are  the 
source  of  beam  instability  and  blow  up  for  the  acceleration  application.  One 
possible  solution  to  this  non-zero  transverse  wake  fields  is  to  use  a 
waveguide,  or  a  cavity,  which  are  partially  filled  with  dielectric  (known  al¬ 
ternatively  as  Dielectric  Wakefields  Acceleration  scheme).  In  this  scheme  an 
intense  driving  charge  excites  the  wake  fields  through  the  Cerenkov  radiation 
mechanism.  This  radiation  is  then  used  to  accelerate  a  second,  less  intense 
particles  to  higher  energies. 

This  method  has  a  particular  advantage  over  other  wake-field  accelera¬ 
tion  schemes  that  the  transverse  wake  fields  can  be  made  quite  small  for 
the  ultra-relativistic  driving  beam.  Recent  experimental16  and  theoretical17 
studies  have  shown  that  the  transverse  wake  potential  in  a  dielectric-loaded 
circular  waveguide  vanishes  in  the  limit  of  the  particle  velocity  v  -*  c  even 
when  the  driving  beam  is  off-center,  and  the  longitudinal  acceleration 
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gradient  is  of  the  same  order  of  magnitude  that  can  be  obtained  in  the 
disk-loaded  metallic  cavity. 

Being  a  co-linear  g^ametry  like  the  circular  waveguide  used  in  that 
experiment,  staging  is  not  quite  attractive.  Since  we  are  primarily  concerned 
with  the  geometry  which  can  permit  easy  staging  in  this  investigation,  the 
wake  fields  excited  in  the  dielectric-loaded  waveguide  of  an  elliptical  cross 
section  are  of  greater  interest.  In  this  chapter,  we  will  investigate  the  wake 
fields  excited  in  the  dielectric-loaded  elliptical  waveguide.  The  property  of 
vanishing  transverse  wake  fields  is  investigated  for  a  use  as  a  wake-field 
acceleration  device.  We  use  the  Fourier  transform  method  to  solve  this 
problem  for  convenience  in  calculations,  and  consequently  all  quantities  are 
treated  in  the  frequency  domain.  It  is  also  assumed  that  the  waveguide  wall 
is  perfectly  conducting  and  the  waveguide  is  filled  with  isotropic  dielectric. 

4.  2.  Formulation  of  Solution 

Consider  the  elliptical  waveguide  which  is  partially  filled  with  medium 
of  dielectric  constant  e  and  permeability  n  as  shown  in  figure  4.1.  Again 
the  problem  is  solved  in  the  confocal  elliptical  coordinates.  The  coordinate 
of  the  waveguide  wall  and  that  of  the  vacuum-dielectric  interface  are, 
respectively,  |  =  |,  and  |  =  lb.  A  point  charge  Q  is  assumed  to  move  in 
a  vacuum  with  velocity  v  along  the  arbitrary  line  (|o,  tjq)  which  is  parallel 
to  the  focus  axis  of  an  elliptical  cavity.  The  electromagnetic  fields  produced 
by  the  motion  of  a  charge  in  a  structure  is  given  again  by  the  Maxwell 
equations.  In  this  chapter,  we  use  cgs  units. 

The  Maxwell  equations  are 


Figure  4.1.  Cross  section  of  the  elliptical  waveguide 
partially  filled  with  dielectric 
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V 


x  E  =  - 


IdE 
c  a/’ 


(4.1) 


V 


x  H  = 


j.  ao 
c  dt  * 


VD  =  Anpt 
V-B  =  0, 

D  =  eE, 
B=//H. 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 


The  electric  field  E  and  the  magnetic  induction  B  can  be  written  in  terms 
of  a  vector  potential  A  and  a  scalar  potential  <p  as 


— & 

(4.7) 

B  =  V  x  A. 

(4.8) 

In  the  Lorentz  gage,  V-A  +  ^^ 
are  transformed  to 


2  _  a2A  _  _  4 Jtp 


VZA 


c2  a/2 


J, 


v> 


c2  a/2 


4;r 

e 


=  0,  the  Maxwell  equations  (4.2)  and  (4.3) 

(4.9) 

(4.10) 


The  charge  density  in  the  elliptical  coordinates  (£,  rjy  z)  is  given  by 
P  =  £?<5(r-  n>)<5(z  -  vt) 

«  £  ag  -  f 0)  <9(>?  -  *,)  <5(z  -  w),  (4.11) 


where  ro  is  the  transverse  coordinate  of  an  exciting  charge  Q,  and  A  is 
the  Jacobian  relation  given  by 
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A-  =  y(cosh2£  -  cos27),  (4.12) 

where  h  is  a  semi-interfocal  distance  as  shown  in  figure  4.1.  Through  the 
continuity  equation,  the  current  density  becomes 

J  -  zpv,  (4.13) 

where  i  is  the  unit  vector  in  the  2-direction.  Since  the  current  density  of 
the  particle  has  only  2-component,  and  the  metallic  waveguide  wall  and 
dielectric  surface  are  smooth,  the  vector  potential  A  has  only  2-component, 
i.e,  A  =  zAz.  As  a  result,  scalar  and  vector  potentials  are  proportional  to 
each  other.  In  fact, 

Az  =  —  v<p.  (4.14) 

c 

Therefore,  we  can  construct  a  solution  with  the  knowledge  of  a  scalar 
potential  <p  alone.  We  are  going  to  Fourier  transform  the  inhomogeneous 
wave  equation  (4.10)  and  solve  in  the  frequency  domain.  In  addition,  the 
dependence  on  z  and  t  are  such  that  all  quantities  depend  on  the  combined 
variable  (z  -  t it). 

First,  we  expand  the  delta  functions  in  eq.  (4.11).  The  well-known 
Fourier  expansion  of  <5 (z  -  vt)  is34 

d{z-vt)  =  ~  Re/0*e'^  *  ">  da,  (4.15) 

where  Re  stands  for  the  real  part  of  the  complex  quantity.  It  is  convenient 
to  expand  <5(7  -  70)  in  terms  of  the  othogonal  Mathieu  functions  of  the  first 
kind,  cem(7,  - q ),  as 
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l 

d(rj  -r}o)  = -2  cem(rj,  - q )  cem(rjo,  -q),  (4. 16) 

m“0 

where  q  is  given  by  the  relation 

=  y(£)2(l-W?!).  (4.17) 

where  p  =  v/c.  Here,  we  use  the  othogonal  property  of  the  Mathieu 
function  cem(j],  -q), 

jc  ft  ~4)  -4)  =  (4.18) 


to  derive  eq.  (4.16),  and  dmo  is  a  Kronecker  delta.  Substituting  eqs.  (4.15) 
and  (4.16)  into  (4.11),  the  charge  density  becomes 

p  =  2 cem(ri,  -q) cem(t] o,  -q)  X  Re ft e'^1 ' w) dco.  (4.19) 

Similary,  the  scalar  potential  <p  can  be  expanded  in  harmonics  of  the  form 

1  " 

<P(£,  V,  z  -  vt)  =  3  2  cem(rj,  -q)  cem(rjo,  -q) 

Ji  m  -0 

x  Re  ft  e'^z '  v°  a) )  d(o.  (4.20) 

From  the  definition  of  <p  given  by  eq.  (4.20),  the  Laplacian  operator  V2  in 
the  elliptical  coordinates  ($,  rj,z)  becomes 


V’-  V*  +  a\  +  -  l(*  +  *\_f “V 

32  32 

and  the  differential  operators  — j  and  — j  become 

dt  dz 


(4.21) 


(4.22) 
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(4.23) 


Substituting  eqs.  (4.19)  -  (4.23)  into  (4.10),  the  equation  for  Om(£,a>)  is 
reduced  to 


*-  -  A  (7>!  t1  -  £^!) 4>-  =  -  ' i(4  -  ?o).  (4.24) 


Since,  the  Mathieu  function  cem( 7,  -<j)  is  a  solution  of  the  Mathieu 
equation  with  negative  q ,  y"  +  (a  +  2q  cos2rj)y  =  0,  the  second  term  in  eq. 
(4.24)  can  be  expressed  in  terms  of  parameters  a,  q,  and  rj.  Finally,  we 
obtain  the  following  equation  for  <I>m(£,  a>): 


<o)-(a  +  2q  cosh2$)  <K„(|,  ^  <5({  -  !o).  (4.25) 


Eq.  (4.25)  together  with  (4.7),  (4.8),  (4.14)  and  (4.20)  completely  determines 
the  fields  produced  by  the  moving  charge,  from  which  we  can  determine 
the  longitudinal  and  transverse  wake  potentials. 

4.  3.  Wake  Fields  Calculation 

4.  3.  1.  Solution  to  the  Boundary  Value  Problem 

Two  independent  solutions  of  the  corresponding  homogeneous  equation 
of  (4.25)  are  the  modified  Mathieu  functions  of  the  first  kind  Cem(£,  -q) 
and  the  second  kind  FekmQ -q),  which  are  given  in  Appendix  A.  Here, 
we  denote  the  vacuum  region  as  region  1  and  the  dielectric  medium  as 
region  2.  Let  us  first  find  the  particular  solution  of  inhomogeneous  Mathieu 
equation  (4.25)  in  the  vacuum  region  where  £  < 


In  region  1,  we  have 
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e=H  =  \,  /S2sl, 

^  =  T(^)2(1_^)-°-  (4’26) 

Tne  particular  solution  of  eq.  (4.25)  can  be  found  by  using  the  Green 
ft  action  technique,  i.e.,  <I>m  is  continuous  while  the  derivative  of  <J>m  has  a 
jump  at  |  =  £0.  Furthermore,  if  we  demand  Om  be  finite  at  £  =  0  and 
vanish  as  £  -*  »  for  all  values  of  q,  then  the  solution  is 

On,,  inhomo  =  -77^  CV?m(£<  ,  -4)  Fe/:m(£>  ,  -q) ,  (4.27) 

UtoV 

where  subscript  <  (or  >)  is  the  smaller  (or  the  bigger)  of  £  and  £0,  and 
Dm  is  a  Wronskian,  which  is  independent  of  £  and  is  a  function  of  q  only. 
That  is 

Dm{i i)  =  Cem(£o,  ~q)  Fekm'(> £0,  ~q)  -  Cea'(£ 0,  -q)  Fekm(£ 0,  -?) 

=  Cem(0,  -q)  Fekm'iO,  -<?),  (4.28) 

where  f  =  <i//<i£,  and  Cem'( 0,  -<?)  =  0  is  used. 

In  order  to  form  a  complete  set  of  solutions  which  satisfy  all  the 
boundary  conditions,  we  have  to  have  a  homogeneous  solution  of  eq.  (4.25), 
which  is  given  by 

Om>homo  =  Cm  C£m(£o»  ~q)  CCm(£,  — <?)»  (4.29) 

where  Cm  is  an  arbitrary  constant.  The  solution  of  the  form  Fefcra(£,  -q)  is 
not  considered  in  eq.  (4.29)  for  the  same  reasons  discussed  in  chapter  3.1., 
i.e.,  the  product  -q)  Fekm(£,  -q)  must  satisfy  two  additional  conditions 
given  by  eqs.  (3.11)  and  (3.12).  We  see  that  it  is  continuous  but  its 
derivative  with  respect  to  £  is  not  continuous  on  the  interfocal  line. 
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Combining  eqs.  (4.27)  and  (4.29),  <bm  in  region  1  can  be  written  in  the 
form 

«Vl  “  -7 r^T  -?)  Fel:m(Z>,  -q) 

4“  Cm  C<?m(|o,  —  (j)  ~ <?)]•  (4.30) 

In  a  dielectric  medium  (|b  <  |  <  £«),  we  have 
q-i,  c/32  >  1, 

?  =  f(f)2(  l-c^2)=-i<0.  (4.31) 

The  solution  in  region  2  is  simply 

^m.2  =  Am  Feym(£,  s )  +  Bra  Cem(|,  s),  (4.32) 

where  Am  and  Bm  are  arbitrary  constant.  We  use  the  parameter  s  instead 
of  q  in  this  region.  Feym(!,$),  the  y-type  modified  Mathieu  function  of  the 
second  kind,  is  used  as  a  second  solution  in  (4.32)  for  convenience  in 
calculations. 

The  coefficients  Cm,  Am  and  J9m  are  determined  from  the  boundary 
conditions  at  the  metallic  wall  (|  =  |a)  and  at  the  vacuum-dielectric  interface 
(I  =  |b).  The  boundary  conditions  are: 

(i)  £t  =  0  at  |  =  I*. 

(ii)  Et  is  continuous  at  |  =  |b. 

(iii)  D(  is  continuous  at  |  =  |b. 

From  eqs.  (4.7),  (4.14)  and  (4.20),  we  see  that  Ez  and  D{  are  proportional 
to 
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4.  3.  2.  Longitudinal  and  Transverse  Wake  Fields 

The  longitudinal  electric  field  can  be  calculated  from  the  scalar  poten¬ 
tial  directly.  From  eqs.  (4.7),  (4.14)  and  (4.20),  the  electric  field  in  region 
1  becomes 
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£*(£.  ?,  z  - 1>/)  =  -^  (  1  -  ^2)  S  cem(?7  -  (?)  cem(f]o,  ~  q)  Re  S*da)  iai 

JtV  m  “0 

x  Cef:-q)  -«)  +  C.  C«„({.  -?))  «*'  *  (4.41) 

Dm\q) 

For  a  test  charge  at  position  t  =  (z  +  zo)/c,  or  trailing  the  driving  charge 
by  zo,  the  longitudinal  electric  field  experienced  by  a  test  charge  is  given 
by 


£*(£,  *I,zo)  ~  (1  _  P2)  £  cem(>7,  -  <?)  cem(r)Q,  -  q )  Re  /"dw  m 

JtV  m  «0 

X (*».(&  -9)  +  c. Cs.g,  -9)) 

=  2  £)m).  (4.42) 

ra-0 

This  is  the  wake  fields  left  by  the  source.  The  longitudinal  wake  potential 

Wz,  defined  as  the  energy  change  of  a  test  charge  per  unit  driving  charge, 

becomes  simply 

W,  =  ^-lL0E,dz.  (4.43) 

where  L  is  the  structure  length. 

The  wake  fields  excited  by  the  ultra-relativistic  particle  are  of  particular 
interest.  Consider  the  case  when  the  particle  velocity  v-*  c,  or  P~*  1.  From 
eq.  (4.26),  we  see  that  the  parameter  q  becomes  zero  as  p  -*  1.  In  this 
case,  we  can  use  the  following  Lmiting  properties  of  the  Mathieu  functions 
to  simplify  eq.  (4.42).  As  q  -*  0,  we  have 


where  ^0)  and  are  given  in  eqs.  (4.39)  and  (4.40).  Substituting  eq. 
(4.49)  into  (4.42)  then  gives  the  longitudinal  electric  field  as 
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£?’(£,  7.  *>)  = 


Re/r  dm  io) 


^0)  + 


s  sinh2|b 


(4.50) 


Here,  the  integral  involving  Feko  can  be  neglected  because  the  integrand 
involving  Feko  is  real  and  has  no  poles  when  /?2  -*  1.  Thus  it  does  not  make 
a  contribution  to  the  real  part  of  the  integral  as  fi2  -*■  1  and  it  is  reduced 


(1  -/92)  Re  /"  do)  io)  e'*7*  -*  0. 


(4.51) 


The  real  part  of  the  integral  in  eq.  (4.50)  comes  from  the  contr'butions 
from  the  poles  of  the  integrand.  The  integrand  is  purely  imaginary  and  the 
real  part  can  come  only  from  the  residues  at  the  poles.  From  complex 
analysis,  we  know  that 


I « 


(4-52) 


From  eq.  (4.31),  we  have  the  relation 

d  ds  d  2s  d 
d(o  da)  ds  ds  ’ 


(4.53) 


then  the  longitudinal  electric  Geld  of  the  monopole  mode  can  be  written 


as  an  infinite  sum  over  all  discrete  harmonics  as 


£<0)  _  12.  y _ _ |  c0™ 

eh 2  T  d  ^o)  ,  jsinh^b^  c 


m 

COS — Zo, 


(4.54) 


where  sx  are  the  roots  of 


vS»  +  h^°)  =  o. 


(4.55) 
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From  eq.  (4.31),  the  resonant  frequencies  a)x  are  related  to  sx  by 

“  =  T02(e_1)-  (456) 

Notice  that  there  is  no  £  or  rj  dependance  of  £l0)  in  eq.  (4.54).  It  is 
apparent  from  eq.  (4.54)  that  there  is  no  transverse  wake  fields  associated 
with  m  =  0  mode  in  the  limit  of  v  =  c,  since  the  transverse  wake  fields 
are  related  to  the  transverse  gradient  cf  Et  through  the  Panofsky-Wenzel 
theorem,  eq.  (2.34). 

Similarly,  Et  for  m  >  1  can  be  found  from  eq.  (4.42)  as 

£(*»)  =  (i  -  p2)  2  cem(ijt  - q )  cem(tjo,  - q ) 

JtV  m  "0 

x  Re  /;  dm  im  Ce^g~?)  (FeM£.  -q)  +  C„  Cemtf ,  -«))  e'?“ 

”  ^2  C1  -  A  £  ««(9.  -9)  «■(<*  -q) 

7CV  m-0 

x  Re  f’  dm  im  Ce„(|o,  -?)  Ce„(£,  -?)  (4.57) 

Again,  by  using  the  limiting  values  given  in  eqs.  (4.44)  -  (4.48),  we  get 


(1  _/)2)  ~q)  Ce"«.  -9) 


-(1  -)S!) 


coshm^o  cosh/ng  . 
m  sinhmgb  exp(mgb)^ 


and  consequently  the  longitudinal  electric  field  of  the  multipole  mode 
becomes 


^  .task  c°sm’'  c°shmi  x  (i  ~f)  Re;" d<u  ■“  ** 


0,  as/9  -*  1. 


(4.59) 


Here,  we  found  that  only  monopole  mode  (m  =  0)  is  excited  in  a 
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dielectric  loaded  elliptical  waveguide  in  the  limit  of  the  particle  velocity 
v  -*  c.  The  contributions  from  multipole  modes  become  negligible  for  the 
ultra-relaticistic  particles  because  Ez  tends  to  zero  as  y~2. 

The  transverse  wake  fields  can  be  calculated  from  the  longitudin*  l 
wake  field  by  using  the  Panofsky-Wenzel  theorem,  eq.  (2.34): 

Awx  =  VXWZ  =  VxeEz  (2.34) 


where  Vx  in  the  elliptical  coordinates  is  given  by 


Vx  = 


h  Vsinhz£  +  sin^ 


d  ,  .  a , 


(4.60) 


As  discussed  earlier,  there  is  no  transverse  wake  fields  associated  with 
m  =  0  mode  because  Ei0)  is  independent  of  the  transverse  coordinates  £ 
and  tf.  For  m  *  0  modes,  the  transverse  wake  fields  are  found  by  substitut¬ 
ing  eq.  (4.57)  into  (2.34).  Again,  using  the  limiting  properties  in  eq.  (4.64) 
-  (4.68),  we  have 

/■Wx  =  VxeEz 

oZq 

Of 

■*  constant  x  Vx(cosm7  coshml)  x  (1  -yS2)  Re  /“  dco  ico  e ,,«2# 

-*  0  as/9  -*  1.  (4.61) 

We  see  again  that  in  the  vacuum  region  the  transverse  wake  fields  become 
neglegible  for  the  ultra-relativistic  particles  and  tend  to  zero  as  y~2. 

For  the  driving  beam  with  an  arbitrary  charge  distribution  of  A(z),  the 
longitudinal  electric  field  becomes 


eh2  T  _d  ^o)  +  sinh2|b 


„/!.«)  cos^-zo).  (4.62) 


Let  us  consider  the  following  two  limiting  cases.  First,  if  the  waveguide 
is  completely  filled  with  a  medium  (|b  =  0),  then  eq.  (4.54)  is  reduced  to 


£?>  =  % 

eh 


s )  Feyo(0, s)  -  CeoCQ,  s)  Fey o(|„ s)  oa 

_ _ rfCafe.*)  "4  (463) 


Feyo'(Q,s) 


where  sx  are  roots  of  Ceo($i,s)  =  0.  Here,  Feyo'(Q,  s)  has  no  pole  since 


Feyo'(0,  j)  =  constant  x  ceo(~,  s)  >  0 


(4.64) 


for  any  s.  When  ^“0  but  finite,  we  see  that  the  fields  are  little  different 
from  those  in  a  completely  filled  waveguide  because 


dCe o(|b,  s) 


(4.65) 


for  |b  =  0.  The  important  fact  is  that  a  charge  moving  close  to  a  medium 
radiates  as  if  it  were  moving  through  a  medium,  and  the  fields  in  the 
vacuum  is  not  very  sensitive  to  the  vacuum-dielectric  boundary  conditions. 

When  lb-* |a,  we  see  that  both  ^ra)  andWjm)  in  eq.  (4.54)  become 
zero.  Therefore,  no  wake  fields  are  generated  in  the  perfectly  conducting 
smooth  waveguide,  as  discussed  previously  in  section  2.1. 


4.  3.  3.  Transition  to  the  Dielectric-Loaded  Circular  Waveguide 


When  the  boundary  ellipse  degenerates  into  a  circle,  we  get  the  equa¬ 
tion  for  the  wake  fields  excited  in  a  dielectric-loaded  circular  waveguide. 
We  can  show  this  again  by  using  the  limiting  properties  of  the  Mathieu 
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functions  given  in  Appendix  A: 

s)  Pm  Jm(kr), 

(4.66) 

Cem'(£, s)  -* pm  krjm'(kr). 

(4.67) 

Feym(%,  s) -*  pm  Nm(kr), 

(4.68) 

Feym'(£, s)  -* pmkr Na '(kr). 

(4.69) 

cem(tj,  -s)  -*  pm  cos m<p. 

(4.70) 

where  pm  is  a  constant  multilpier,  and  /m,  Nm  are  respectively  the  Bessel 
functions  of  the  first  and  second  kinds.  The  relation  between  5  and  k  is 
given  by 


(4.71) 


From  eq.  (4.71),  we  have  the  relation 


d  _  2  d 
ds  kh2dk‘ 


If  confocal  ellipses  with  §  =  and  §  =  £b  become  concentric  circles  of 

radius  r  =  a  and  r  =  b  respectively  as  shown  in  figure  4.2,  then 
substituting  eqs.  (4.66)  -  (4.72)  into  (4.54)  gives  the  longitudinal  electric 
field  £l0)  excited  in  the  dielectric-loaded  circular  waveguide: 


c<o)  =  2£  y 

^  e*  t 


cos^’ 


(4.73) 


where 


/f°  =  /o(faz)  -  7o(W>)  No(faz), 

/S0)  =  Joika)  No’ (kb)  -  Jo' (kb)  N0(ka), 


(4.74) 

(4.75) 
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which  is  exactly  the  same  equation  as  in  ref.  18.  Here,  Jo  and  No  are  the 
zero  order  Bessel  functions  of  the  first  kind  and  second  kind  respectively. 

4.  4.  Numerical  Example  and  Discussion 

As  an  example,  the  longitudinal  wake  potential  in  the  dielectric- maded 
elliptical  waveguide  is  calculated  for  the  limiting  case  of  v  =  c.  The  dimen¬ 
sions  of  the  waveguide  are  shown  in  figure  4.3,  and  the  driving  beam  with 
a  Gaussian  charge  distribution  of  a  =  5  mm  is  assumed.  Figure  4.4  is  a 
plot  of  the  longitudinal  wake  potential  in  the  vacuum  when  the  dielectric 
constant  is  £  =  2  and  vacuum  hole  is  very  narrow,  i.e.,  &>  555  0.  Here,  we 
obtain  about  40  MeV/m  acceleration  gradient  per  1  micro-Coulomb  driving 
charge,  while  about  120  MeV/m  was  obtained  in  chapter  3  for  the  disk- 
loaded  cavity  which  has  the  same  cross  sectional  area.  We  found  that  in 
a  dielectric-loaded  elliptical  waveguide  only  monopole  mode  (m  =  0)  is  ex¬ 
cited  and  the  transverse  wake  fields,  which  exist  in  the  disk-loaded  elliptical 
cavity,  become  zero  as  the  particle  velocity  fi -*  1.  This  result  is  very  im¬ 
portant  when  we  consider  using  the  dielectric-loaded  waveguide  as  a  wake- 
field  acceleration  device.  It  can  support  high  accelerating  gradient,  but  has 
no  deflecting  force  which  might  lead  to  beam  instabilities.  The  amplitude 
of  the  longitudinal  electric  field  is  a  function  of  the  transverse  position  of 
the  driving  beam.  Since  the  longitudinal  electric  field  is  uniform  within  the 
vacuum  region,  we  can  inject  the  driving  beam  and  accelerated  beam  at 
two  different  locations  in  the  vacuum  region  so  that  staging  becomes  prac¬ 
tical. 


conductor 


h  coshf. 


(Zjh  cosh|b  -*  b  *1  •*  <P 


Figure  4.2.  Transition  to  the  dielectric-loaded  circular  waveguide 


Igure  4A-  Longitudinal 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

Acceleration  of  charged  particles  by  the  wake  fields  produced  by  a 
bunch  of  relativistic  charged  particles  passing  through  the  elliptical  wake-field 
cavities  has  been  investigated.  The  expressions  for  the  wake  potentials  have 
been  derived  analytically  for  the  disk-loaded  metallic  cavity  and  dielectric- 
loaded  waveguide.  For  the  disk-loaded  elliptical  cavity,  modal  analysis 
method  is  used  under  the  assumption  of  no  aperture  in  a  cavity,  and  the 
wake  fields  and  wake  potentials  are  derived  in  terms  of  the  Mathieu  func¬ 
tions.  It  is  found  that  these  analytical  results  agree  well  with  the  numerical 
analysis  of  Chin,22  in  which  finite  difference  method  is  used  to  directly  solve 
the  Maxwell  equations  in  the  time  domain.  For  a  dielectric-loaded  elliptical 
cavity,  the  longitudinal  wake  potential  is  derived  analytically  using  Fourier 
transform  method.  It  is  shown  that  both  disk-loaded  and  dielectric-loaded 
cavities  can  support  an  acceleration  gradient  in  excess  of  a  few  hundred 
MeV/m.  It  is  also  found  that  strong  transverse  wake  fields,  which  could  be 
useful  for  the  FEL  application  as  an  electric  wiggler,  are  excited  in  the 
disk-loaded  cavity,  while  no  transverse  wake  fields  are  induced  in  a 
dielectric-loaded  cavity  for  the  ultra-relativistic  driving  beam. 

The  idea  of  wake-field  acceleration  by  means  of  the  elliptical  cavity 
offers  a  promising  way  for  reaching  acceleration  gradients  of  a  few  hundred 
MeV/m.  The  wake-field  cavities  discussed  in  the  previous  chapter,  do  not 
represent  the  maximum  acceleration  gradient  that  can  be  achieved.  'T'he  lon¬ 
gitudinal  wake  potential  or  acceleration  gradient  is  proportional  to  the  num- 
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ber  of  particles  in  the  driving  beam,  and  inversely  proportional  to  the  cross 
sectional  area  of  the  cavity.  Taking  the  parameters  used  in  example  calcula¬ 
tions,  an  elliptical  cavity  with  the  major  axis  of  5  cm  and  the  minor  axis 
of  3  cm  wil  yield  600  MeV/m/uC.  Although  the  transformer  ratio  of  this 
elliptical  cavity  is  not  as  high,  a  wake-field  accelerator  not  based  on  the 
impedance  transformation  principle  can  be  achieved  by  using  multi-stage 
scheme  in  which  short  stages  are  used  with  the  driving  beam  which  is 
replaced  with  new  one  or  whose  energy  is  replenished  aftei  each  stage.  For 
the  elliptical  cavity,  we  can  simplify  the  replacement  of  decelerated  driving 
beam  considerably.  The  need  for  a  fast  kicker,  or -bending  the  driver  beam 
can  be  eliminated.  Furthermore,  the  fact  that  the  wake-field  devices  are 
simply  shaped  metallic  pill-boxes  or  waveguides  partially  filled  with  dielectric 
without  any  frequency  tuning  or  phasing  devices  make  this  kind  of  ac¬ 
celerator  unexpensive  and  reliable. 

The  limitations  to  achieving  a  high  acceleration  gradient  should  be 
pointed  out.  The  principle  limitation  is  due  to  breakdown  and  electrical  dis¬ 
charges  which  may  damage  the  metallic  surface,  dielectric  breakdown,  and 
change  of  dielectric  properties.  We  do  not  concern  ourselves  with  this  prob¬ 
lem  about  which  little  is  known  yet.  However,  breakdown  limits  are  con¬ 
sidered  to  be  less  severe  than  any  other  near-field  accelerating  devices  since 
the  high  field  strength  exists  only  a  limited  area  and  only  for  a  short 
period  of  time.  Another  difficulty  is  the  production  of  very  short,  high  cur¬ 
rent  driving  beam.  A  charge  of  1  f iC ,  which  corresponds  to  •'bout  6  x  101* 
electrons  within  the  driving  beam  or  a  peak  current  of  nearly  60  KA,  could 
be  achieved  in  the  induction  linac.  The  laser-driven  photoinjector  may  allow 
a  simpler  solution  to  this  problem. 
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As  mentioned,  we  have  made  several  simplifying  assumptions  in  order 
to  obtain  the  analytical  solutions  of  the  wake  potentials  in  the  elliptical 
cavity.  Further  investigations  are  needed  on  the  following  issues:  (1)  Beam 
aperture  effects  are  not  considered  for  the  disk-loaded  cavity.  Although  we 

have  demonstrated  that  these  analytical  results  are  in  good  agreement  with 
the  results  of  the  numerical  code  for  the  cavity  with  small  apertures,  large 
aperture  will  affect  the  calculation  in  both  magnitudes  and  frequencies  of 

the  wake  potentials.  (2)  The  wake  fields  inside  a  driving  beam  can  not  be 

accurately  calculated  with  this  method  due  to  the  slow  convergence  of  the 

modes  summation,  therefore  we  cannot  calculate  the  energy  loss  of  the  driv¬ 
ing  beam  accurately.  (3)  Beam  dynamics  calculation  and  instability  study  due 
to  the  transverse  deflecting  modes  in  the  elliptical  pill-box  cavity.  (4)  Non¬ 
linear  effects  in  the  dielectric.  A  systematic  numerical  study  will  give  a  bet¬ 
ter  picture  of  these  wake  fields  problem. 


APPENDIX  A 

THE  MATHIEU  EQUATIONS  AND  MATHIEU  FUNCTIONS 
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A,  1.  Canonical  Form  of  the  Mathieu  Equation 

The  Mathieu  equation  is  the  particular  case  of  a  linear  type  of  fhe 
second  order  differential  equation  with  periodic  coefficients,  which  has  a 
form 

+  (a  -  2q  cos2 z)y  -  0,  (A.1) 

dz 

where  the  parameters  a  and  q  are  limited  to  real  number,  but  z  is  usually 
unrestricted.  For  the  present,  we  shall  confine  our  attention  to  solutions 
having  period  x  or  2n  in  z.  If  we  write  /£  for  z  in  eq.  (A.1),  it  becomes 

(a- 2q  cosh2£)y  =  0  (A.2) 

A.  2.  Periodic  Solution  of  y"  +  (a  -2q  cos2z)y  =  0 

We  have  two  types  of  solutions,  one  is  even  and  the  other  is  odd 
function  in  z.  These  solutions  may  be  expressed  in  series  forms  as: 

ce2n(z,  q)  =  'Z  A%n)  cos2 rz,  (A.3) 

r-o 

ce2n+i(z,  q)  =  2  cos(2r+ 1  )z,  (A.4) 

r-o 

se2n+2(z,  q)  =  £  M2r"+22)  sin(2r+2)z,  (A.5) 

r»o 

q)  =  2  11}  sin(2r+ 1  )z,  (A.6) 

r*o 

for  n  =  0,  1,  2,  ...  In  these  series,  the  coefficients  A and  Brm)  are 
functions  of  q.  The  function  cem(z,q )  and  sem(z,q )  are  respectively  the  even 
and  odd  type  Mathieu  functions  of  the  first  kinds  of  integral  order  m.  A 


85 


few  of  the  lowest-order  Mathieu  functions  are  plotted  in  figures  A.l.  and 
A.2. 

A.  3.  Solution  of  y"  —  (a  -  2^cosh2|)y  =  0 

The  first  solutions  of  eq.  (A.2)  are  derived  by  substituting  /£  for  z  in 
eqs.  (A.3)  -  (A.6).  Thus,  we  have 


Ce 2o(i,  q )  =  cernig,  q)  =  2  A%n)  cosh2r£,  (A.7) 

r»o 

00 

Ce*+i@,  q)  =  «!„♦,(;?.  q)  =  cosh(2r+l)f,  (A.8) 

r»o 

S<WS,  q)  =  -l  «&«(<£.  <?)  =  £  sinh(2r+2)|,  (A.9) 

r-o 

Se*+i(lq)  *  -iseia+ii&q)  =  £  B&t1*  sinh(2r+ 1)|.  (A.  10) 

r*o 


These  are  defined  to  be  modified  Mathieu  functions  of  the  first  kind  of 
integral  order,  and  the  capital  letters  are  used  to  denote  the  modified 
functions. 

A.  4.  Othogonality  and  Normalization  of  Mathieu  Functions 

The  Mathieu  functions  of  the  first  kinds  have  the  property  that  for  a 
given  q , 

/J*  cem(z,  q)  cep(z,  q)dz  =  S*  sem (z,  q)  sep(z,  q)  dz  =  0  for  m*p.  (A.  1 1) 

With  this  othogonal  property,  various  normalizations  are  possible.  Here, 
cem(z,q )  and  sem(z,q )  are  normalized  according  to 


J  /f  «*2(z,  9)  J  ?)<*  =  ! 


(A.  12) 


for  all  real  values  of  q.  The  normalization  in  (A.11)  entails  a  mean  squre 
value  of  1/2  over  the  interval  (0,  2jt),  which  is  the  same  as  that  for  the 
circular  functions.  Under  this  normalization,  the  coefficients  A and  f$m) are 
given  by  the  relatior 


2  [A^f  +  2  \a%a)]2  =  1, 

r-l 

2  [A&V']2  =  i  [B&Xl)]2  =  2  [Mrn+22)]2  =  1. 
r»0  r-0  r-0 


(A.  13) 
(A.14) 


A.  5.  Degenerate  Forms  of  Mathieu  Functions 


A.  5.  1.  Mathieu  Functions  of  the  First  Kind 


When  the  fundamental  ellipse  tends  to  a  circle  by  letting  h  -*  0,  hence 
q  —  0,  Aim)  and  Mm)  in  the  series  tend  to  zero  except  that  ylSl)  -*  1  and 
-*  1,  and  the  confocal  hyperbolae  become  radii  of  the  circle  with 
rj  =<p.  So,  we  have 


cem(T},q) 


cosmrj  =  cos m<p  for  m  S  1 
A^  =  1  /VI  for  m  =  0 


sem(ri,  q)  —  sin mr\  =  sin/n#  for  m  2:  1. 


A.  5.  2.  Modified  Mathieu  Functions 


(A.15) 
(A.  16) 


When  a  confocal  ellipse  of  semi-major  axis  r  tends  to  a  circle  with 
this  radius  by  letting  h  -*  0,  £  -*  »  while  keeping  the  product  h  cosh£  -*  r, 
then  a  -*  m2  for  a  function  of  integral  order  m.  Consequently,  eq.  (A.2) 
degenerate  to  the  standard  Bessel  equation: 


(A.  17) 
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where  2 q  =  h2  k\/2.  Thus  we  have  the  following  degenerate  cases: 


Ce m(£,  q)  Pm  Jmikir), 

(A.  18) 

Seai%,  q)  -  5m  Jmikir), 

(A.19) 

•JZ  Cem(£,  q)  ■* Pm  kirJm'ikir), 

(A.20) 

Semi£,  q)  -*  Sm  kxr  Jmikir), 

(A.21) 

Fey  mi!;,  q)  -Pm  Ymikir), 

(A.22) 

-jgFeymi!;,  q)  -*Pm  kxr  Ym'ik\r), 

(A.23) 

Ce2«iS,  ~q)  pin  hnik\r), 

(A.24) 

Ce2n+ii£,  —q)  -*  szn+i  I2o+iik\r), 

(A.25) 

Se 2n+l(£,  -q)  Pln+l  Jin+iikxr), 

(A.26) 

Se 2n+2(^,  -q)  s2n+2 12t,-nik\r), 

(A.27) 

ji  Fek2ai£,  -q)  -  p2n  K^ihr), 

(A.28) 

n  Fek2n+ii%,  —q)  -*  s2n+i  K2n+\ik\r), 

(A.29) 

where  pm ,  sm  are  constant  multipliers,  and  fix)  =  df{x)/dx.  Feym{%,  q)  is  a 
linearly  independent  second  solution  for  positive  q.  For  negative  q ,  the 
solution  of  the  form  Fekmi%,  -q)  is  usually  used. 

A.  6.  Limiting  Forms  of  Feym(l q)  and  Fekm{!-,  -q)  as  q  -*  0 

Feyo(l  q)  ■*  ^  (y  -  ln2  +  $  +  V2  \nq)  (1  -  l/2  q  cosh2£)  (A.30) 

2Feko(S,  ~q) 

■*  -  l(y  -  ln2  + 1  +  V2  In q)  (1  +  V2q  cosh2|)  -  V2  sinh2|]  (A.31) 

Feym(S,  q)  ■*  -  22m~l  Tim)  T(m+ 1)  jT1  q~a  e~mi  for  m  >  1  (A.32) 

Fekmfa  -q)  -  Y2Feym($,q)  for  m  S  1  (A.33) 
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Figure  A.l.  Mathieu  function  of  the  first  kind  of  even  integral  order 

cem(z,  q)  and  sem(z,  q) 
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GLOSSARY 


a  Separation  constant 

a*(x)  ector  eigenfunction 

A(x,  t)  ector  potential 

Am  Arbitrary  constant 

A$m)  Coefficients  for  the  series  representation  of  cem(rj,q )  and  Cem(£,  q) 

B(x,  t)  Magnetic  flux  density 

Bm  Arbitrary  constant 

l£m)  Coefficients  for  the  series  representation  of  q)  and  Sem(£,  q) 

c  Speed  of  light 

cem(r}t  q)  Even  type  Mathieu  function  of  the  first  kind  of  integral  order  m 
Cem(£,  q)  Modified  Mathieu  function  of  the  first  kind  of  integral  order  m 
Cm  Arbitrary  constant 

Cmo  Coefficient  in  the  longitudinal  wake  potential 

Cmn  Coefficient  in  the  transverse  wake  potential 

d  Gap  distance  of  an  elliptical  pillbox  cavity 

D(x,  t)  Electric  flux  density 

Dm(q)  Wronskian  relation 

An(<7)  =  Cem( 0,  -q)Fekm'(0,  -q)  -  Cea\ 0,  -q)Fekm{ 0,  -q) 

ec  Eccentricity  of  the  boundary  ellipse 
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E(x,  t )  electric  field 

Feym(£,  q)  Y-type  Mathieu  function  of  the  second  kind 

Fekm(£,  -q)  K-type  Mathieu  function  of  the  second  kind 
h  Semi-interfocal  distance  of  an  ellipse 

H(x,  t)  Magnetic  field 

i  r 


J 

Jm(x ) 

• 

k 

ke 

kz 

L 

min(/,  Lie) 

n 

Nm(x) 

Pm 

®(0 

?ran 

Q 

r. 

n> 

n(t) 


Current  density 

Bessel  function  of  the  first  kind  of  integral  order  m 
Wave  vector 

Transverse  component  of  wave  vector  k 
Longitudinal  component  of  wave  vector  k 
length  of  a  cavity 
Minimum  of  t  and  Lie 
Unit  normal  vector 

Bessel  function  of  the  second  kind  of  integral  order  m 
Constant  multiplier 

Time  dependent  coefficients  of  the  vector  potential  A(x,  t) 

Parameters  in  the  Mathieu  functions 

Charge  of  a  point  source 

Outer  radius  of  dielectric-loaded  circular  cavity 

Radius  of  dielectric-vacuum  interface  of  dielectric-loaded  circular  cavity 

Time  dependent  coefficient  of  the  scalar  potential  <p(\t) 


to 

Ro 
Re  2 
s 

sem{rj,q) 

Sem(£t  q) 

Sm 

S 

t 

Tx 

ux 

Ux(zo) 

Uz(z0) 

V 

V 

Vx{r) 

Wx(z0) 

Wz(zo) 

X 

(*>  y,  z) 

Xb 


Transverse  coordinate  of  the  driving  charge 
Radius  of  a  circular  pillbox  cavity 
Real  part  of  complex  number  z 

Parameter  related  to  the  transverse  wave  number  in  dielectric 
Odd  type  Mathieu  function  of  the  first  kind  of  integral  order  m 
Modified  Mathieu  function  of  the  first  kind  of  integral  order  m 
Constant  multiplier 

Cross  sectional  area  of  an  elliptic  pillbox  cavity 
Time 

Normalizing  factor  of  the  scalar  eigenfunction  <px 
Normalizing  factor  of  the  vector  eigenfunction  ai 
Transverse  wake  potential  of  distributed  charges 
Longitudinal  wake  potential  of  distributed  charges 
elocity  of  a  particle 
olume  of  a  cavity 

oltage  induced  by  the  driving  charge  at  the  transverse  position 
Delta  function  transverse  wake  potential 
Delta  function  longitudinal  wake  potential 
Spatial  coordinates 
Cartesian  Coordinates  ' 

Semi-major  axis  of  the  boundary  ellipse 


yb 
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Xoo 
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Bo 
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^(X.0 
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m 

M 

Mo 

p 
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Semi-minor  axis  of  the  boundary  ellipse 
Unit  vector  in  z-direction 

Delay  distance  between  the  driving  and  test  charge 
Normalized  velocity  /9  =  v/c 
nth  root  of  Bessel  function  Jo(x) 

Kronecker  delta 
Jacobian  A  = 


.  l) 


Permittivity  of  dielectric 

Permittivity  of  vacuum 

1/2  for  p  -  0  and  1  for  p  *  0 

Scalar  eigenfunctions 

Scalar  potential  in  the  frequency  domain 

Scalar  potential 

Relativistic  factor  y  =  (1  -  {Pyxn 
Line  charge  density  of  the  driving  bunch 
Permeability  of  dielectric 
Permeability  of  vacuum 
Charge  density 

rms  length  of  a  Gaussian  bunch 
Resonant  frequencies 
Confocal  elliptical  coordinates 


(l  *l>z) 
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lb 

(lo»  rjo) 


£  at  the  boundaiy  ellipse 
£  at  the  dielectric-vacuum  interface 
Transverse  coordinates  of  the  driving  charge 


